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The auxiliary field method is a new technique to obtain closed formulae for the solutions of 
eigenequations in quantum mechanics. The idea is to replace a Hamiltonian H for which analytical 
solutions are not known by another one H, including one or more auxiliary fields, for which they are 
known. For instance, a potential V(r) not solvable is replaced by another one P{r) more familiar, or a 
semirelativistic kinetic part is replaced by an equivalent nonrelativistic one. If the auxiliary fields are 
eliminated by an extremization procedure, the Hamiltonian H reduces to Hamiltonian H. The approx- 
imation comes from the replacement of these fields by pure real constants. The approximant solutions 
for H, eigenvalues and eigenfunctions, are then obtained by the solutions of H in which the auxiliary 
parameters are eliminated by an extremization procedure for the eigenenergies, which takes the form 
of a transcendental equation to solve. If H = T(p) + V(r) and if P{r) is a power law, the approximate 
eigenvalues can be written T(p{)) + V(ro), where the mean impulsion po is a function of the mean dis- 
tance ro and where is determined by an equation which is linked to the generalized virial theorem. 
The general properties of the method are studied and the connections with the envelope theory pre- 
sented. Its mean field and (anti)variational characters are also discussed. This method is first applied to 
nonrelativistic and semirelativistic two-body systems, with a great variety of potentials (sum of power 
laws, logarithm, exponential, square root). Closed formulae are produced for energies, eigenstates, var- 
ious observables and critical constants (when it is relevant), with sometimes a very good accuracy. The 
method is then used to solve nonrelativistic and semirelativistic many-body systems with one-body and 
two-body interactions. For such cases, analytical solutions can only be obtained for systems of identical 
particles, but several systems of interest for atomic and hadronic physics are studied. General results 
concerning the many-body critical constants are presented, as well as duality relations existing between 
approximate and exact eigenvalues. 
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§1. Introduction: Analytical methods in quantum mechanics 

The aim of this work is to present the auxiliary field method (AFM), which is a new and 
remarkably simple method to find analytical approximate solutions of eigenequations corre- 
sponding to Hamiltonians admitting bound states. The quest for exactly solvable problems 
in quantum physics is actually as old as quantum mechanics itself, and it is worth making 
general comments about this topic before focusing on the auxiliary field method. 

The most famous systems for which the Schrodinger equation admits analytical solu- 
tions aie probably the harmonic oscillator and the hydrogen atom, i.e. Hamiltonians of the 
form + and — l/r respectively, where p and r are conjugate variables. The spectrum 
of the one-dimensional harmonic oscillator is actually known since 1925 thanks to Heisen- 
berg's pioneering work,^ ' while the spectrum of the nonrelativistic hydrogen atom was found 
in Schrodinger's papei^' as a first and crucial test of his celebrated equation. Notice that the 
same problem was solved the same year by 

Since those early results, there has been a considerable amount of works devoted to 
the computation of analytical solutions of the Schrodinger equation, especially in bound 
state problems. The Schrodinger equation has been found to be exactly solvable in many 
one-dimensional cases: Dirac comb, exponential and linear potentials, . . . However, only a 
few analytical three-dimensional solutions are known for any value of the angular momen- 
tum (the ^-wave channel is very similar to a simpler one-dimensional equation). Besides 
the aforementioned harmonic oscillator and Coulomb cases, one can mention the Kratzer 
potential, the particle in a box and the symmetric top, which are of quite limited use in mod- 
eling realistic quantum systems. An exhaustive discussion of quantum problems admitting 
analytical solutions would be out of the scope of the present introduction; therefore we refer 
the interested reader to the textbooks,'^''^''^''^''^''2J in which a wide range of problems is 
covered. 

When exact analytical solutions cannot be found, physicists can either resort to nu- 
merical computations or to methods leading to approximate analytical solutions. Even if 
the computational power that we have at our disposal nowadays allows to solve accurately 
nearly every eigenequation for few body systems, finding approximate analytical formulae 
is always useful in physics, not only because of the intrinsic mathematical interest of such 
a task, but also in view of obtaining informations about the dependence of the observables 
on the parameters of a model and on the quantum numbers of the states. Calculating an 
analytical expression being always much less time consuming than solving numerically the 
corresponding eigenequation, the use of closed formulae can be a great advantage when one 
tries to fit the parameters of a given model to some experimental data, or even to interpret 
new data. We find worth mentioning here three approximation schemes that are all "textbook 
material" and able to lead to analytic results: the WKB or semiclassical approximation, the 
variational method, and the perturbation theory. 

The principle of the WKB approximation has been found by Wentzel, Kramers and 
Brillouin around 1926,1^ but it is also known that Jeffreys brought important contributions 
to the field in 1924.^ It basically consists in a Taylor expansion of the Hamiltonian and the 
wavefunction in powers of h. At dominant order for a one-dimensional problem, such an 
expansion leads for example to the Bohr-Sommerfeld quantization rule Jl^l^} p{E ,x)dx = 
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7r/j(?i + 1/2), where p is the classical momentum expressed as a function of the energy and 
of the position x, and where a{E) and b{E) we. the turning points of the classical trajectory. 
The energy spectrum can then be extracted from this integral and is generally accurate, 
but the explicit expression for E{n) is not always possible to find. One has instead often 
to numerically compute the energy from the Bohr-Sommerfeld quantization rule. Further 
informations about the semiclassical approximation can be found in e.g. Refs. |4]|,|7]l. 

While the WKB method is particularly efficient for highly excited states, where % 
clearly becomes a small parameter, the variational method due to Ritz^^' allows to gain 
relevant informations about the ground state (for a modern presentation see Ref. 13)). In- 
deed, let us consider a Hamiltonian H and a normalized arbitrary function i/a belonging to 
space. Then, the Ritz theorem states that Eq < where Eq is the exact ground-state 

energy. Provided that i/a is chosen such that is analytical, one can then obtain an 

analytical upper bound on the ground-state energy. Typically, one chooses a gaussian- or 
exponential-type wavefunction; many examples of analytical upper bounds can be found in 
Ref. |9]l. Notice also the MacDonald's theorem stating that, given a set of n orthonormal- 
ized functions i//,- belonging to space, the eigenvalues of the matrix y/y) provide 

upper bounds on the n first eigenenergies of Z/.^*^' However, analytical results can hardly be 
obtained from such a generalization due to the complexity of the underlying calculations. 

A last technique to be mentioned is the perturbation theory.'^' Let us consider that the 
Schrodinger equation with a potential V is exactly solvable. Then the Schrodinger equation 
with potential V + can be solved thanks to a Taylor expansion in provided that this 
parameter is small. At first order in particular, the computation of {w) with the lowest- 
order states provides the dominant correction to the whole energy spectrum. Again, a wide 
range of quantum problems can be studied by using perturbation theory, and we refer to 
g gl3j9j explicit examples in which (w) is analytical. 

Those three methods being presented, it appears that a major challenge remains: To 
find approximate analytical solutions of eigenequations for the complete spectrum (not only 
the ground state as with the variational method) without using any Taylor expansion as done 
with the WKB method or the perturbation theory. The AFM, that will be presented in detail 
in the rest of this work, is an attempt to address that issue. 

§2. Generalities on the auxiliary field method 

2.1. Historical aspects 

Auxiliary fields, also known as einbein fields, are used in various domains of theo- 
retical physics. Historically, they have been introduced in order to get rid of the square 
roots typically appearing in relativistic Lagrangians.^^'^ The most obvious example is the 
free relativistic particle, described by the Lagrangian ^ = —m V^, where x is the world- 
velocity. An auxiliary field /i can be introduced in this expression so that one is led to the 
more convenient form ^{pi) = x^/ (2/x) + m^jx /2. This last expression is formally simpler 
than Lagrangian , but it is equivalent if the equations of motion of the auxiliary field are 
considered. The auxiliary field is indeed given on shell by }Xq = —yf^jm, and it is readily 
checked that ^(/^) = ^ . That feature can be seen as a definition of an auxiliary field: It is 
a field whose equation of motion is not dynamical but leads to algebraic relations that allow 
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to express it in terms of the other degrees of freedom of the problem. 

Besides pointhke particles, auxiliary fields have soon become a cornerstone of string 
theory because of the systematic replacement of the Nambu-Goto Lagrangian =S^g = 

—a f dQ^ (ww')2 — w^w'^, in which an unconvenient square root appears, by the Polyakov 

Lagrangian = —{a/2)J_dd V— det h hf"^ dbW^ dcW ^, where the induced metric /j'"' is 
made of auxiliary fields.'^^''^''^ In connection with string theory, supersymmetric field 
theories, linking bosons and fermions through supersymmetric transformations, also gener- 
ally demand the introduction of extra fields - auxiliary fields - in order to close the super- 
symmetric algebra. More information concerning the building of e.g. consistent supergrav- 
ity theories can be found for example in Refs. I22l).l23l).l24l).l25l). 

More generally, auxiliary fields have become widely used in field theory. Indeed, they 
allow to replace kinetic terms of the form \/^ by a formally nonrelativistic term such as 
y? / (2jU). Computations in the path integral formulation of field theories are then simplified 
because one is then formally led to Gaussian path integrals, about which analytical results 
are known.'^''^ In quantum chromodynamics (QCD) for example, the introduction of aux- 
iliary fields has allowed to gain relevant informations about confinement in hadrons within 
the framework of potential models^ '^ and in particular to give support to the phenomeno- 
logical QCD string model.^"''^' It is worth mentioning that renormalization problems^ 
as well as numerical computations in A/^— body systemd^J^D can also be addressed using 
auxiliary fields as a computational tool. 

A particular interest for our purpose is to come back to effective models of QCD. It is 
known that, in the strong coupling limit, the Wilson loop formulation of QCD supports a 
linear confinement between a heavy quark and antiquark.l^ That result can be generalized 
to light quarks also; it then appears that the effective confining interaction in mesonic sys- 
tems can be modeled by a Nambu-Goto string linking the quark with the antiquark,^^ ' whose 
dominant contribution is indeed a linear potential of type a r. An idea that is already present 
in Refs. I28l).|29l) is to introduce an additional auxiliary field, v, in a way similar to what is 
done for relativistic kinetic terms: ar can be replaced by a^r^ / (2v) + v/2, formally reduc- 
ing the linear potential to a harmonic oscillator, for which analytical solutions are known. 
Of course the remaining auxiliary field has eventually to be eliminated, but analytical ap- 
proximate mass formulae for mesons and baryons can be found following this philosophy. 
The question underlying the present report is nothing else but an attempt to generalize such 
a procedure, namely: Could any arbitrary potential appearing in an eigenequation be ap- 
proximated by an expression involving an auxiliary field and another potential for which an 
analytical solution is known? As we shall see in the following, an affirmative answer can 
be given to that question, leading to a method allowing to compute approximate analytical 
energy formulae of various eigenequations. Moreover the method can be generalized to treat 
problems involving more than 1 or 2 particles. 

2.2. The auxiliary field method 

The most general form that we consider for the eigenvalue equation of a one- or two- 
body problem in first quantized quantum mechanics is 



H\yif) = [T{p) + V{r)]W)=EW). 



(2-1) 
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The potential V in the Hamiltonian H depends on the variable r = | r| which is the distance 
from the center of forces in a one-body problem and the relative distance between the two 
particles in a two-body problem. The kinetic part T depends on the momentum operator 
p which is the conjugate variable of r. Practically, we will consider a nonrelativistic form 
p^/ (2m) (m is the particle mass for a one-body problem or the reduced mass of a two-body 
system) as in the Schrodinger equation and a semirelativistic one a\/p^ + (a = 1 for 
one-body problems and a = 2 for two identical particles) as in a spinless Salpeter equation. 

An exact analytical expression for all the eigenvalues is known explicitly only for 
very specific potentials P{r) with a nonrelativistic kinetic parts: the quadratic interaction 
P(r) = (harmonic oscillator) and the Coulomb potential P{r) = —l/r (hydrogen-like 
system) are the most familiar for a practical use. An exact expression is known for a number 
of other potentials, but only for ^-waves (see Ref. "9^ for a more detailed discussion) or for 
one-dimensional problems. 

Our goal is the search for approximate analytical solutions for the Hamiltonian of type 
(12-11) . relying on Hamiltonians for which solutions are well-known. In other words, we 
assume that we are able to obtain an analytical solution for the equation 

h{a)\a) = [T{p) + aP{r)] \a) = e{a) \a) , (2-2) 

in which, at this stage, a is a real parameter. This method was formulated and applied for 
the first time in Ref. [38l) . 

We summarize here the principle of the method introducing an auxiliary field, which is 
a priori an operator. It consists in four steps: 

1. We calculate the function K{r) (the prime denotes the derivative with respect to the 
argument of the function) 

V'(r) 

2. We denote by 7 = K^^ the inverse function of K. Thus, one has 

K{r) = z and J{z) = r. (2-4) 

Since both V{r) and P(r) do exhibit an analytical form, the same property holds for 
K{r). But it is by no means sure that J{z) can be expressed analytically. We will see in 
the following that an explicit analytical expression for J is not necessary to write down 
the basic equations. 

3. We introduce the "bridge function" B defined by 

B{y)=V{J{y))-yP{J{y)). (2-5) 

The name of this function comes from the fact that it makes a bridge between the 
potential P{r) for which an analytical expression for the energies is known and the 
potential V (r) for which the corresponding analytical expression is a priori not known. 
Lastly, it is interesting to define the AFM potential V through 

V{r,y)=yP{r)+B{y). (2-6) 

It depends not only on the position operator as the original potential V{r), but also 
on the y variable, which is undetermined at this stage. The AFM Hamiltonian is now 
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defined as 

H{y) = T{p)+V{r,y) = h{y) +B{y). {1-1) 

The construction of H seems rather artificial and one can wonder what could be the 
justification of such a procedure. In fact, there is no mystery. 

Let us suppose that the y variable is an operator denoted v. It plays the role of an 
auxiliary field, not present in the original Hamiltonian H, but which is an essential 
ingredient of the AFM Hamiltonian H. Among the infinite number of possibilities for 
the auxiliary field v, let us choose a very specific one Vq defined by 

%=K{r). (2-8) 

A more explicit notation would be Vo(r) indicating clearly that it is an operator de- 
pending on the position only. The very important property is that Vq, coming from 
(I2-8I) . is an extremum the AFM Hamiltonian, i.e. 



5H{v) 



5v 



= 0. (2-9) 

V=Vo 



Moreover, one has the additional property that the value of the AFM Hamiltonian taken 
for this operator coincides with the original Hamiltonian H 

H{vo)=H. (2-10) 



These properties are easy to show and follow from the definitions (I2-3I ). (l2-4b and (12-81) . 
Thus, considering the auxiliary field as an operator and affecting it the value given by 
(12-81) is just an alternative method to solve the original problem. 
4. The philosophy of the AFM is to consider the y variable, no longer as an operator v, 
but as a pure real number v. In this case the AFM potential V{r,v) is of the form 
vP(r) +B{v) where v and B{v) are no longer operators but must be considered as 
arbitrary constants. Taking into account (l2-2b . the eigenvalues of ^(v) are 

E{v)=e{v)+B{v), (2-11) 

where e{v) are the eigenvalues of h{v) which are supposed to be known whatever the 
radial and orbital quantum numbers n, 1. Provided that the J function is calculable, the 
same property is true for E{v). Then, we determine the value Vq that extremizes ^(v): 



dE( 



dv 



= 0. (2-12) 

V=Vo 



We propose to consider ^(Vo) as the approximate form of the exact energy E of the 
Hamiltonian H 

E^Eafm = E{vo). (2-13) 

In order to obtain an analytical expression for the eigenvalues E^fm, we must then 
fulfill a second necessary condition: to be able to determine Vq and, then, ^(Vo) in an 
analytical way. Denoting |v) the eigenstate of ^(v) corresponding to the eigenvalue 
^(v), we consider |vo) as an approximation of the corresponding genuine eigenstate 
ofH. 
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Presented as such, this procedure appears to be an empirical recipe. Arguments given in 
the third item show that this prescription makes sense and that the AFM must be considered 
in essence as a mean field approximation. This point will be developed more deeply below. 

Another interesting property concerning the AFM potential (12-61 ) is the following. Let 
us denote by ro the value of the radius defined by 

Vo=i^(ro) and ro=7(vo). (2-14) 

It is just a matter of simple calculation to verify that the AFM potential taken for the special 
value y = Vo can be recast as 

V(r,Vo) = Vo[P(r)-P(ro)]+V(ro). (2-15) 

Owing to the definition of the K function (I2-3I ). one can easily show that the AFM potential 
coincides both to the exact potential and its derivative at the point r = ro, namely 

V(ro,Vo) = y(ro) and t/'(ro, Vq) = V'(ro). (2-16) 

In other words, at this particular point ro, the AFM potential and the exact potential are 
tangent curves. This property will be exploited in Sect. I2.4l in connection with the envelope 
theory (see Appendix lAl). 

Let us note by |v) an eigenstate of H{v). The Hellmann-Feynman theoreni^Sj states 

that 

dE{v) _ 
dv 

Using this relation, one can show thaP^ 

(vo|P(r)|Vo)=P(ro). (2-18) 

This means that ro is a kind of "average point" for the potential P{r). That is why it will be 
often called "the mean radius" in the following. Using this last relation with the definitions 
above, we get 

(vo|Z(vo)|vo)=Z(vo) with Z{x)=P{J{x)). (2-19) 

So, our method can actually be considered as a "mean field approximation" with respect 
to a particular auxiliary field which is introduced to simplify the calculations: Vo is the 
mean value of the operator Vq = K{r) through a function Z which can be quite simple. For 
example, Z{x) =xifV{x) = P{xY /2 + Vb where Vb is a constant. 

It is in the step of passing from an operator to a constant interpretation for the auxil- 
iary field that lies the approximation. It is clear that the quality of the approximate results 
strongly depends on the choice of the function P{r). The cleverness of the physicist relies 
in his abiUty to guess a form for P{r) as close as possible of V(r) while leading to a man- 
ageable form of the corresponding eigenvalues e(v). A big bulk of this report is devoted to 
the discussion of the quality of this approximation. 

This method is completely general and a priori valid for any potential V{r). Of course, 
all the four steps mentioned above can be done numerically so that ^afm can be computed 
numerically. But it seems that the corresponding numerical treatment is as heavy as solving 



dH{v) 
dv 



(2-17) 
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directly the eigenvalue equation to get lastly a poorer result. This is true if we start with an 
arbitrary function P{r). However we have the complete freedom for this choice. 

Let us assume that, for a kinetic part T, we do start with a potential P allowing to obtain 
analytical results for its eigenvalues e{v) ie{v) is a shorthand notation for the most correct 
expression e(v; «,/)); the same is true of course for its derivative with respect to v : e'{v). 
Using the definition of the K and J functions defined previously, it is easy to show that 

E'{v)=e'{v)-P{J{v)). (2-20) 

The determination of Vq results from the condition ^'(Vo) = 0, that is to say 

e'ivo)=P{J{vo)). (2-21) 

It appears that obtaining the Vq value requires only to solve a transcendental equation. Once 
this value is obtained, the AFM energy is very easy to calculate 

Eafm = E{Vo) = e{vo)+B{vo). (2-22) 

Thus, provided that we have an analytical expression for e{v), the AFM results only need 
solving a transcendental equation, a much easier procedure than solving exactly the genuine 
eigenvalue equation. 

The practical usefulness of the AFM method is to get a final solution which is com- 
pletely analytical. In order to do that, we must fulfill, as we saw, three conditions: 

1. to choose a P function leading to analytical expression for e{v); 

2. to be able to invert relation (12-81) in order to have access to the function J defined by 
(l2-4b (it appears into the B function); 

3. to be able to determine Vq (through (12-211) ) and to calculate the corresponding value 
^(Vo) in an analytical way. 

In fact, the second condition is formally not necessary. Indeed, let us introduce the mean 
radius ro = /(Vq). Then Vq = K{ro) so that the bridge function is expressed as 

B{vo) = Viro)-K{ro)Piro) = B{ro), (2-23) 

where, for simplicity, we still maintain the notation B: B = V — KP. The transcendental 
equation (12-211) is transformed into 

e'iKiro))=Piro), (2-24) 

while the energy is given by 

EAFM = Eivo)=eiK{ro)) + B{ro). (2-25) 

Using the expression (12-231) for S(ro) and the definition of e{v), this last formula can be 
recast under the form 

EAFM = To + V{ro), (2-26) 

where, according to (12-181) . 



To = e{K{ro)) - K{ro)P{ro) = (vo| T{p) | Vq) 



(2-27) 
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is an average kinetic energy. 

In this new formulation, any reference to the J function has disappeared so that it is 
not necessary to try to get it. However, solving analytically the transcendental equation 
(12-241) is more or less of the same difficulty than getting an analytical expression for the / 
function. Anyhow, this new formulation leads very often to simpler practical calculations 
and is preferred most of the time. 

It is important to stress the following point: if P{r) is chosen to be V{r), it is trivial to 
check that K{r) = 1, B{r) = 0, ro is meaningless since obviously P{r) and V{r) are tangent 
curves everywhere, and Safm = ^(1) which is precisely the exact value E. In this particular 
case, AFM recovers the exact result. This property is sometimes used to check special 
results. 

2.3. Scaling laws 

Scahng laws represent an important property for quantum mechanical systems. They 
allow to give the expression for the eigenenergies (and wavefunctions) of the most gen- 
eral equation in terms of the corresponding eigenenergies (and wavefunctions) of a reduced 
equation with less parameters. In fact, the scaling laws are nothing else than a direct conse- 
quence of dimensional analysis applied to the various dimensioned parameters entering the 
problem. Starting from a general Hamiltonian //(«!,...,«„) depending on n dimensioned 
parameters a,, it is generally possible to write 

//(«!,...,«„) = 7/2(i8i,...,A>.), (2-28) 

where y has the dimension of an energy and where /j is a Hamiltonian expressed in terms 
of dimensionless conjugate variables and depending on m < « dimensionless parameters j8,-. 
The scaling properties of the nonrelativistic Schrodinger equation have been studied and 
used for a long time, but spinless Salpeter equations also benefit of scaling properties. They 
will be explicitly described below. In the following, in order to lighten the notations, the 
AFM will be generally applied to dimensionless Hamiltonians as in Refs. 38,).,40j.,4l]). So, 
the exact and the AFM solutions automatically share the same scaling properties. 

2.4. Eigenstates and upper/lower bounds 

If £'(Vo), which is an eigenvalue of ^(vq), is an approximation of the exact energy E, 
the corresponding eigenstate |vo) is an approximation of a genuine eigenstate of H. The 
shape of the corresponding wavefunction (r|vo) depends on the quantum numbers via the 
parameter Vq or the mean radius tq. Practically, only nonrelativistic harmonic oscillator or 
nonrelativistic Coulomb wavefunctions can be used. Explicit examples will be presented 
below. 

Since the potential V{r, Vq) is tangent to the potential V{r) at r = ro, and since H and 
^(vo) have the same kinetic part, the comparison theorem^'ES''^^ (for both nonrelativistic 
and relativistic equations) implies that the approximation ^(Vo) is an upper (lower) bound 
on the exact energy if ^(r, Vq) > V{r) {V{r,Vo) < V{r)) for all values of r. Equivalently, a 
function g{x) can be defined by 

Vix)=g{P{x)). (2-29) 

It can then be shown that, if g{x) is a concave (convex) function, that is if g"{x) < {g"{x) > 
0) \/x, the approximation ^(vo) is an upper (lower) bound on the exact energy. This property 
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has been demonstrated in the framework of the envelope theory, ' but can be apphed as well 
to the AFM.**^' Several examples will be presented below. 

The knowledge of lower and upper bounds on an eigenstate is a first technique to esti- 
mate the accuracy of the AFM. It has also been shown thaP^ 

E{vo) - {vo\H\vo) = V{ro) - (Vo|V(r)| Vq). (2-30) 

The right-hand side of this equation is the difference between the value of potential V com- 
puted at the average point ro and the average of this potential for the AFM state | Vq) consid- 
ered here as trial state. In some favorable cases (the trial state is a ground state for instance), 
E < {vo\H\vo) and a bound on the error can be computed by 

E{Vo)-E > Viro) - (vo|V(r)|vo). (2-31) 

If the mean value of V can be computed analytically, this constitutes a second procedure 
to estimate the accuracy of the AFM. Several examples of this calculation are presented in 
Ref. I38l) . At last, the eigenstates of a Hamiltonian of type (12-11) can be solved numerically 
with an arbitrary precision. So, as a third possibility, a direct comparison with the AFM 
results is always possible. So, one can wonder why to use the AFM? Let us recall that the 
interest of this method is mainly to obtain analytical information about the whole spectra 
(dependence of eigenenergies on the parameters of the Hamiltonian and on the quantum 
numbers), without necessarily searching a very high accuracy. Moreover, the AFM approx- 
imation can be extended to A/^-body problems for which exact eigenenergies are not easily 
reachable, even numerically. 

In some cases, it is possible to compute analytically the mean value £'*(Vo) = {vo\H\vo), 
considering the AFM solution | Vq) as a trial state (see Sect. l3.3l for two examples). We have 
then 

E*{vo)-E{vo) = {Vo\H-H{vo)\vo) = {Vo\V{r) -Vir,Vo)\vo) = AV. (2-32) 

If g{x) is concave, E{vo) > E and V{r, Vq) > V{r). In this case, AV <0 and E*{vo) < 
E{vo)- By the Ritz theorem, we know that £'*(vo) > E for the ground state. For this state, it 
is interesting to compute the mean value £'*(v*) = {v*\H\v*), where |v*) is an eigenstate of 
^(v*) with V* determined in order to minimize the energy E* (v*). We have then £'*(v*) < 
£'*(vo) since the parameter Vq is fixed by the AFM computation. 

If g{x) is convex, E{vo) < E and V(r, Vq) < V(r). Thus, AV >0 and E*{vo) > E{vo). 
In this case, we have also £'*(vo) > E and £'*(v*) < £'*(vo) for the ground state. 

We can gather these results to obtain: 

• If g{x) is concave or equivalently V{r, Vq) >V{r): 

- For the ground state, E{vo) > E*{vo) > E*{v*) > E; 

- For the other states, E{vo) > E and £(vo) > E*{vo). 

• If g{x) is convex or equivalently V{r, Vq) <V{r): 

- For the ground state, E*{vo) >E*{v*)>E> E{vo); 

- For the other states, E > E{vo) and E*{vo) > E{vo). 

The results mentioned above are directly applicable if H and ^(vq) have the same 
kinetic part. This is always the case for nonrelativistic Hamiltonians but not for the spinless 
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Salpeter Hamiltonians, whose treatment requires generally the replacement of the square 
root operator by a nonrelativistic operator. The bounds on these kinds of Hamiltonians are 
specifically studied in Sect. \5\ 

2.4. 1 . Extension to the form aP{r) + V (r) 

If the AFM is applicable for some potentials V\ (r) and V2 (r) independently, there is 
no certainty that it still applies for a potential which is their sum: V{r) = Vi{r) +V2{r). 
However, the method can be used in the case where one of the potentials V, can be identified 
with the basic P potential. Thus, in this section we consider a Hamiltonian of type 

Ha = T{p) + aP{r)+V{r), (2-33) 

whose exact eigenvalues are denoted Ea- The extension of AFM to potentials of this type was 
first presented in Ref. 40). One introduces an auxiliary field v as before, forgetting about the 
aP{r) contribution. The first 3 steps of the algorithm remain unchanged. Thus the Vq field 
is the same, as is the same the function B{v). The only difference arises in the expression 
(12-71) of H and h where vP{r) has to be replaced by {a + v)P(r). As a consequence, the 
corresponding energy (12-111) has to be replaced by 

Ea{v) = e{a + v)+B{v). (2-34) 

^^(v) is an eigenvalue of Hamiltonian 

Ha{v)=h{a + v)+B{v), (2-35) 

where h is defined by (12-21) and where e{a + v) is an eigenvalue of Hamiltonian h{a + v). 
An eigenstate of Hamiltonians Ha and h{a + v) is denoted |a + v), and we have e{a + v) = 
{a + v\h{a + v)\a + v). If Vq is the value of V which extremizes (12-341 ). then we could expect 
that 

Ea{vo) = e{a + Vo)+B{vo) (2-36) 

is a good approximation of Ea, an exact eigenvalue of Hamiltonian (12-331 ). It seems that 
(12-341 ) is very similar to (12-1 II) . the only difference being the replacement of v by a + v in 
the argument of the function e(v). Nevertheless, this small difference is important, because, 
even if the determination of Vq from (12-1 II) is technically easy and analytical, it may happen 
that its determination from (12-341 ) could be much more involved and very often not analyt- 
ical. The alternative formulation in terms of ro is also slightly modified in this case. The 
transcendental equation now writes 

e'ia + K{ro))=Piro) (2-37) 

and the AFM energy is given by 

EAFM = e{a + K{ro))+B{ro). (2-38) 
Using again the Hellmann-Feynman theorem,'22j' it can be shown that 



{a + Vo\P{r)\a + Vo) = P{J{vo)) = P{ro). 



(2-39) 
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So, /(Vo) = ro is also in this case a kind of "average point" for the potential P(r). Recalling 
that Z{x) =P{J{x)), we have 

(a + Vo|Z(v)|fl + Vo)=Z(vo). (2-40) 

We confirm also in this case the "mean field approximation" with respect to a particular 
auxiliary field which is introduced to simplify the calculations. 

2.5. Recursion procedure 

An idea for obtaining analytical expressions of the eigenenergies for an arbitrary poten- 
tial is the following (see also Ref. I40l)). We start with a potential P{r) = /'[*'l(r) for which 
the energies of the corresponding Hamiltonian are exactly known. We then proceed as 
above to find approximate solutions for the eigenenergies of a Hamiltonian //['^ in which the 
potential is at present V{r) = p['l(r). In general, a large class of potentials can be treated in 
that way. Moreover, as we will see below, by comparison with accurate numerical results, 
we can even refine the expressions so that analytical forms for the energies are very close to 
the exact solutions. 

Considering now these approximate expressions as the exact ones, we apply once more 
the AFM with P{r) = P^^\r) to obtain approximate solutions for the eigenenergies of a 
Hamiltonian //PI in which the potential is at present V{r) = pPl(r). Even if analytical 
solutions for Hamiltonian //^l were not attainable directly with P{r) = pl^'l (r), it may occur 
that they indeed are with P(r) = p[^l(r). Pursuing recursively such a procedure, one can 
imagine to get analytical solutions for increasingly complicated potentials. Presumably, the 
quality of the analytical expressions deteriorates with the order of the recursion. 

§3. Schrddinger equation with a power-law potential 

In this section, we discuss in detail the case of the Schrodinger equation with a power- 
law potential for two reasons. First, it is a simple case to illustrate the AFM. Second, as we 
will see below, this kind of potentials is in the thick of our method. The Hamiltonian can be 
written ^ 

H=^ + ^gn{X)ar^, (3-1) 
2m 

where a is a strictly positive constant and where the sign function is defined by sgn(A) = 
A /I A I with A 7^ 0. The physical values of A must be such that A > —2, otherwise the 
wave equation leads to a collapse. The two unavoidable starting potentials P{r) = r^ and 
P{r) = —\/r are indeed two particular cases of power-law potentials. In fact, the only 
interesting values studied in this paper are those comprised between these extreme values: 
— 1 < A < 2. They represent most of the physical applications. For example, the linear 
potential A = 1 is the traditional form of the confining potential in hadronic physics,'23J32)J3D 
the value A = 2/3 was shown to give the good slope for Regge trajectories in nonrelativistic 
treatmenP^ and A = 0.1 was considered by Martin.!^ 

The energy E of (13-11 ) depends on the physical parameters m, a, A but also on the 
quantum numbers n (radial) and / (orbital). It is easy to show that the scaling laws allow to 
write 

E{m,a,X;n,l) = 2~-a7^nC~ e{X;n,l), (3-2) 
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where e{X;n,l) is the eigenvalue of equation 

//=^ + sgn(A)r\ (3-3) 

This form is chosen in order to match the conventions of Ref. [38]l. The scaling laws are 
unable to say anything about the dependence of e on the power A and on the quantum 
numbers {nj). It is the virtue of the AFM to shed some light on these very important 
questions. 

An another interesting potential is the logarithmic potential, discussed for instance in 
Ref.EB 

V{r)=a\n{br). (3-4) 

It is in strong connection with the power-law interaction since it can be rewritten into a 
similar form 

1 5 

ln;t;= lim-(x^-l). (3-5) 
Using the the scaling law, we can write 

2Z?^ tncL 
E{m,a,b) = e{l5) with j8 = — 2, (3-6) 

where e(j3) is the eigenenergy of the reduced Schrodinger Hamiltonian 

h = ^+l5lnr. (3-7) 

We will apply the AFM for these two interactions with the two starting potentials P(r) = 
and P{r) = — l/r. Let us note that a number of results presented in the following sections 
have been already obtained in the framework of the envelope theory.133 Nevertheless we 
remind them in order to keep some consistency in our presentation. 

3.1. Energies with P{r) = 

The harmonic oscillator (HO) Hamiltonian h{v) = /A + vr^ has the following eigen- 
values 

e{v) = ^/vQHo with Qho =^n + l + 3/1. (3-8) 

In the following, a combination of quantum numbers as in Quo (Qho is a simplified notation 
for QHo{n,l)) will be denoted "principal quantum number". 

Let us apply the AFM for the power-law potential V{r) = sgn(A)r'*'. We have then 
K{r) = |A|/-V2, e{K{r)) = QhoVW^^^^'^^'^ ^ ^'i^ir)) = Qnor^^-^^^V and 
B{r) = ((2 — X)/2)V {r). The transcendental equation (12-241) for the extremization of the 
energy leads to the value of the mean radius 



ro 



HO 



l/(A+2) 

(3-9) 



2\X\ 

Inserting this value in E{ro) (12-251 ) gives the AFM approximation of the eigenenergy 



£AFM(A;«,/) = ^|A|5Tx2-2fx2^. (3-10) 
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It is important to emphasize that this value is an upper bound on the exact energy for the 
domain of A we aie interested in. Moreover, putting A = 2 in the previous expression, one 
recovers the exact result eAFM(2;«,/) = e{2;n,l) = Qho< as expected. 

Let us consider now the logarithmic potential V{r) = j3lnr. We have then K{r) = 
j8/(2r2), e{K{r)) = QnoVm/r, e'{K{r)) = rQHo/V2p and B{r) = V{r)-^/l. The 
transcendental equation (12-241) for the extremization of the energy leads to the value of the 
mean radius 

Qho ^^\ 
r() = —^- (3-11) 



2^ 

Inserting this value in E{rQ) (12-251) gives the AFM approximation of the eigenenergy 



£AFM(i8;n,/) =jS In 



'^Qho 



(3-12) 



It is well known that the eigenenergies E{m, a, b; n, I) resulting from a Schrodinger equation 
with the potential (13-41 ) satisfy the propertj^ ' 



E{am,a,b;n,l) = E{m,a,b;n,l) — — Inot. 



(3-13) 



The immediate consequence is that the corresponding spectrum is independent of the mass 
of the particle. It is remarkable that the basic property (13-131 ) still holds for our approximate 
expression (13-121) . 

It is worth mentioning that the energy formula (13-121) can be understood as a particular 
limit case of formula (I3-2I ). as it was suggested in Ref. 47). Guided by (13-51) . let us consider 
the following Hamiltonian 



(3-14) 



1 



It reduces to the Hamiltonian (13-71) in the limit A — 0. A simple rewriting of formula (13-21 ) 
for a — j3/| A I and m = 2 gives the eigenenergies of Hamiltonian (13-141 ). We have thus 



£(A) 



2 + Ai8iT2 2^ 
2A ' ^ 



X 

I'M. 



no 



A' 



and as expected 



lim£'(A) = j3ln 



(3-15) 



(3-16) 



that is precisely formula (13- 121 ). directly obtained from the logarithmic potential. This con- 
firms the idea that the logarithmic potential can be seen as the limit of a power-law potential 
when A goes to zero. The same conclusion was obtained with the envelope theory (see 
Ref.|52l)). 

3.2. Energies with P{r) = —l/r 

Now we consider the Coulomb (C) Hamiltonian h{v) = — v/r whose eigenvalues 
are given by {Qc is a simphfied notation for Qdn, /)) 

e(v) = --j with Qc = n + l + \. (3-17) 
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Let us apply the AFM for the power-law potential V{r) = sgn(A)r'*'. We have then 
K{r) = \X\r^+\e{K{r)) = -r^^+^{\X\yQl), e'{K{r)) = -r^+' / QD andB(r) = (1 + 
X)V{r). The transcendental equation (12-241 ) for the extremization of the energy leads to the 
value of the mean radius 

r g2 1 1/(1+2) 

Inserting this value in E{rQ) (12-251) gives the AFM approximation of the eigenenergy 

2+A 2 A ^ 

£AFM(A;«,Z) = ^|A|5Tx2-2ne2+\ (3-19) 

It is important to emphasize that this value is a lower bound on the exact energy. Putting A = 

— 1 in the previous expression, one recovers the exact result £afm(— 1;'^>0 = £{—^\n,l) = 

— l/2c' expected. It is worth noting that formulae (13-101 ) and ( |3-19b exhibit the same 
formal expression with just the exchange Qc o Qho- This important point will be discussed 
inSectgll 

The resolution of the Schrodinger equation with the logarithmic potential in (13-71 ) is 
even simpler with the present choice of P{r). Indeed, one finds in this case K{r) = j8r, 
e{K{r)) = -{rp/Qcf, e'{K{r)) = -2^r^/Ql and B{r) = V{r) + jS. The transcendental 
equation (12-241 ) leads to the value of the mean radius 

.0 = ^. (3-20) 
Inserting this value in £'(ro) (12-251) gives the AFM approximation of the eigenenergy 

eAFM(j8;«,/) = i3 



(3-21) 

Again, formulae (13-121 ) and (13-211 ) exhibit the same formal expression with just the exchange 
Qc ^ Qho- Consequently, the properties mentioned above concerning the spectrum and the 
limit of a power-law potential when A — )■ still hold in this case. 

3.3. Wavefunctions and observables 

In this section, the quality of the wavefunction Vq) given by the AFM is tested for the 
Schrodinger Hamiltonian with the linear potential (results for logarithmic and exponential 
potentials are presented in Ref. 53)). This Hamiltonian is chosen because it can be solved 
analytically for a vanishing angular momentum / (see Appendix |G|. When / 7^ 0, numer- 
ical results have been obtained with two different methods.'^''^ In order to simplify the 
notation, we will denote \n,l) an eigenstate of H and \n) = |?i,0), with {r\n) given by (IG-31 ). 
The corresponding energies will be denoted £"„ / = {n,l\H\n,l) and = E„ o = {n\H\n). 
As we need a Hamiltonian H with a central potential which is completely solvable to apply 
the AFM, we can only use in practice a hydrogen-like system iP{r) = — 1/r) or a harmonic 
oscillator iP{r) = r^). As a linear potential seems closer to than —1/r, we can expect that 
the use of a harmonic oscillator to start the AFM will give better results. Using the scaling 
properties, we can consider the following simple Hamiltonian 

H = p^ + r, (3-22) 
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in order to match the conventions of Ref . 

As it is shown in the previous section, the approximate AFM energies, denoted here by 
are given by 

end = :^Q^'\ (3-23) 

with Q = 2n + J + 2>/2 for P{r) = and Q = n + l + \ for P(r) = — 1/r. Exact energies, 
given by (IG-21 ). reduces to £"„ = — a„, where a„ is the {n + l)^ zero of the Airy functions Ai. 
A simple approximation of E„ can be obtained using the expansion (IG-41 ) at the first order 

En = -an^[ — \ U + « 2.811 U + . (3-24) 

It is worth noting that the linear potential is not only a toy model to test the AFM method. 
Effective theories of QCD have proved that it is a good interaction to take into account the 
confinement of quarks or gluons in potential models of hadronic physics.'^''^''^ 

3.3.1. AFM with P(r) = -1/r 

One can ask whether it is possible to obtain good approximations for the solutions of 
a Schrodinger equation with a linear potential by means of hydrogen-like eigenf unctions. 
This will be examined in this section. Using the previous results for f (r) = — 1/r, we find 
Vo = ''o = 2^/^(?i + / + 1)^/^. Exact eigenstates are then approximated by AFM eigenstates 
which are, in this case, hydrogen-like states (IF-31 ) with 

' 2 2^/3 ^ ' 

Such states are denoted |Hy;?i,/) and |Hy;?i) = |Hy;«,0). Using (IF-61 ) and results above, 
it can be shown that (12-181 ) is satisfied, with P(ro) = -1/ro = -l^^l^in + / + 1)^2/3 

also worth noting that (12- 191 ) gives (l/^fv^ = We denote and e^^ = e^Q the 

approximated energies which are given by (13-231 ) with Q = n + l+l. Since V (r, Vq ) — V (r) = 
— (r — ro)^/r < 0, are lower bounds on the exact energies. This can also be determined 
with the function g defined by ( 12-291 ). In this case, g{y) = — l/y with j < 0. The function 
g"{y) = — 2/y^ being positive, g is convex as expected. 

The quantum number dependence of the scaling parameter T] corrects partly the differ- 
ence between the shapes of {r\n,l) and (r|Hy;?i,/). Consequently, (Hy;?i,/|Hy;?i,/') = 5/// 
because of the orthogonality of the spherical harmonics, but (Hy;?i,/|Hy;«',/) / ^m'- Using 
the definition (IH-ll) . we find 

(Hy;«,/|Hy;.',/) =F„^, ((^TTTl)'^') ' ^^'"^^ 

with given by (IH-31) . Table U gives some values of |(IIy;n,Z|Hy;n',Z)p. We can see 
that the overlap is not negligible for n close to n' , but it decreases rapidly with \n — n'\. The 
situation improves when I increases: |(Hy;0,Z|Hy; 1,/)|2 = 0.43, 0.29, 0.22, 0.18, 0.15, 0.13 
for / = ^ 5. 
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Table I. Results for P„,„/,/ = |(Hy;n,/|Hy;«',0|^. Values for / = (/ = 1) are given in the lower-left (upper- 
right) triangle of the Table. P„^n',l = Pn'.n.l ™d ^«.n,/ = 1- 





n = 


1 


2 


3 


n' = 


1 


0.29 


0.028 


0.0039 


1 


0.43 


1 


0.36 


0.036 


2 


0.055 


0.43 


1 


0.39 


3 


0.0097 


0.049 


0.43 


1 



It is interesting to compare £„ ^ with (13-241 ) 

= ^(« + 1)'^' « 1.890(^+1)2/3. (3-27) 

The ratio e^^ /E„ is respectively equal to 0.808, 0.734, 0.7 12 for « = 0, 1 , 2 and tends rapidly 
toward the asymptotic value 3^^^ /k^^^ 0.672. As expected, these ratios are smaller than 
1 since e^^ are lower bounds. Two wavefunctions are given in Fig. [T] We can see that 
the differences between exact {r\n) and AFM (r|Hy;?i) wavefunctions can be large. The 
overlap |(?i|Hy;?i)p between these wavefunctions can be computed numerically with a high 
accuracy. We find respectively the values 0.934, 0.664, 0.298 for « = 0, 1, 2, showing a rapid 
decrease of the overlap. It is worth noting that, asymptotically, {r\n) oc exp (— | r3/2) while 
(r|Hy;?i) is characterized by an exponential tail. Nevertheless, if an observable is not too 
sensitive to the large r behavior, this discrepancy will not spoil its mean value. 




Fig. 1. Normalized wavefunctions, as a function of the variable r in ( I3-22I >. for n = (left) and 1 (right): (r|n) 
(solid line), (r|Hy;«) (dotted-dashed line) and (r|HO;n) (dashed line). 



The various observables |v/„,o(0)P> ('^) and (/j^) computed with the AFM states can 
be obtained using formulae (IF-41) . (IF-61) and (IF-71) for / = with the parameter rj given by 
(13-251 ). Results are summed up in Table [III A direct comparison between the structure of 
exact and AFM observables can be obtained if we remind that the exact ones depend on | a„ | 
which can be well approximated by j8„ (see (IG-41 )). Let us look in detail only at the mean 
value (r). For the exact and AFM solutions, we have respectively 

{n\r\n) = ^ ^ (^) + ^) « 1 .874 + ^) , (3-28) 
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{Uy;n\r\Uy;n) = _(„ + 1)2/3 ~ 1. 890 («+ 1)^/3. (3-29) 

Some observables like |v^n,o(0)P and {p^) are very badly reproduced. Others can be ob- 
tained with a quite reasonable accuracy. Despite the fact that exact and AFM wavefunctions 
differ strongly when n increases, their sizes stay similar. 

Table II. Various observables A computed with the AFM and P{r) = — \/r. 

(Hy;«|A|Hy;M) (Hy;«|A|Hy;«)/(«|A|w) 
(A) (2m = a = 1) n = n = l n = 2 « > 1 

k„,o(0)P 2 4 6 

{r) 1.212 1.101 1.068 1.009 + ^ + O (^^j) 

{r^) im^TW ^-^^^ ^-^^^ ^-^^^ 0.942 + 



{P^) ^^^$7^ 0.808 0.734 0.712 0.672 + ^ + 

5)(n- "' 

24/3 



(/) {^n+5Kn+i)''' ^_g^5 3_ggQ 2.009« + 1.926 + 



Since it is possible to compute analytically the mean value {¥ly;n\p'^ + r\Yly;n) for an 
arbitrary value of the scale factor rj, the relations between the bounds given in Sect. l2.4l can 
be checked. We can verify that E = -a„ > E{vo) = + 1)^^^ and £'*(vo) = + 
1)2/3 > ^(vo) = + 1)^/3. Moreover, for the ground state, we have E*{vo) = ^ > 

E*{v*) = 1^ > £ = -ao > E{vo) = ^. 

The behavior of observables computed with the AFM is similar for values of / = 0, 1 , 2. 
We do not expect strong deviations for larger values of /. This is illustrated with some 
typical results gathered in Table Hill Observables are generally not very well reproduced, but 
this is expected since eigenstates for a linear potential are very different from eigenstates 
for a Coulomb potential. Actually, the agreement is not catastrophic, except for I'/'fi.olO)^ 
and (/j^) as mentioned above. It is even surprising that the AFM with P{r) = — l/r could 
give energies and some observables for a linear potential with a quite reasonable accuracy. 
Moreover, lower bounds on the energies are obtained. 

Table III. Ratios between the AFM results (energies and (r)) with P{r) = — l/r and the exact results, for 
several quantum number sets {n,l). 

I n=0 n=\ n=2 n=3 n=4 n=5 

0.808 0.734 0.712 0.702 0.696 0.692 

1 0.893 0.805 0.767 0.746 0.733 0.724 

2 0.925 0.846 0.805 0.779 0.762 0.749 

(Hy;?i|r|Hy;?i) / {n\r\n) 

1.212 1.101 1.068 1.053 1.043 1.037 

1 1.116 1.118 1.103 1.089 1.079 1.071 

2 1.080 1.110 1.110 1.104 1.096 1.089 
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3.3.2. AFM with P{r) = 

Since the quadratic potential is closer to the linear potential than a Coulomb one, one 
can expect better results with harmonic oscillator wavef unctions. This will be examined in 
this section. Using the previous results for P{r) = r^, we find Vq = 1 / (2ro) = 2^^/^ {2n + 1 + 
3/2)^^/^. Exact eigenstates are then approximated by AFM eigenstates which are, in this 
case, harmonic oscillator states (IE-31) with 

X = vl'^ = 2-^l\2n + l + 3/l)-^l\ (3-30) 

Such states are denoted |H0;«,/) and |HO;?i) = |HO;?i,0). Using (IE-61 ) and results above, 
it can be shown that (12-181 ) is satisfied, with P{n)) = = 2^/^(2?! + / + 3/2)"*/^ It is also 
worth noting that (12- 191 ) gives (l/v^) = l/Vg. This is in agreement with (20) in Ref. 41]) 
(an auxiliary field = 1/v is used in this last reference). We denote e^f and e^*-* = the 
approximated energies which are given by (13-231 ) with Q = 2n + l + 3/2. Since V{r, Vo) — 
V{r) = {r — r^Yl (2ro) > 0, e^f are upper bounds on the exact energies. In this case, g{y) = 

y(y with y > 0. The function g"{y) = — 1/ (4^^^^) being negative, g is concave as expected. 

Since the scaling parameter A depends on the quantum numbers, (HO;?i,Z|HO;«,Z') = 
Sir because of the orthogonality of the spherical harmonics, but (HO;?i,Z|HO;?i',Z) / ^nn'- 
Using the definition (IH-ll ). we find 

(HO;„,/|Ha„',,)=F™((^±|±i)"'), (Ml) 

with F^*^ given by (IH-41 ). Table IIVI gives some values of |(HO;?i,Z|HO;?i',Z)P- We can see 
that the overlap is always small and decreases rapidly with \n — n'\. The situation is even 
better when Z increases: for Z = ^ 5, |(HO;0,Z|HO;1,Z)|2 = 0.029, 0.023, 0.019, 0.017, 
0.014, 0.013. 



Table IV. Results for Pn^j = |(HO;;j,/|HO;;j',/)|1 Values for Z = (/= 1) are given in the lower-left (upper- 
right) triangle of the Table. -P„.„'./ = Pn'.n.l and P„,„j — 1. 

n=0 1 2 3 

n' = i 0023 0.0026 0.00039 

1 0.029 1 0.026 0.0027 

2 0.0036 0.027 1 0.026 

3 0.00064 0.0031 0.027 1 



Let us look at the energies £„ 



HO 



3\^/^ 

e„"° = 3 « + - . (3-32) 



By comparing with (13-241) . we can see immediately that the situation is more favorable 
than in the previous case. The ratio /En is respectively equal to 1.059, 1.066, 1.067, 
for « = 0,1,2. The asymptotic value 3^l^2^l^ /ii^l^ ^ 1.067 is rapidly approached. As 
expected, these ratios are greater than 1 since are upper bounds. Two wavefunctions 
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are given in Fig.[T] We can see that the differences between exact {r\n) and AFM (r|HO;?i) 
wavef unctions ai^e quite small. The overlap |(?i|HO;?i)p between these wavef unctions can 
be computed numerically with a high accuracy. We find respectively the values 0.997, 0.979, 
0.951 for « = 0, 1,2. A value below 0.75 is reached for « = 6 and below 0.25 for « = 14. A 
wavefunction (r|HO;«) is characterized by an Gaussian tail, while we have asymptotically 
{r\n) oc exp (— f?"^^^)- Again, if an observable is not too sensitive to the large r behavior, 
this discrepancy will not spoil its mean value. 

The various observables |i//„o(0)p, {r^) and {p'^) computed with the AFM states can 
be obtained using formulae (IE-41) . (IE-61 ) and (IE-71 ) for / = with the parameter A given by 
(13-301 ). Results are summed up in Table IVl Again, a direct comparison between the structure 
of exact and AFM observables can be obtained by using j8„ (see (IG-41 )) instead of | a„ | . Let 
us look in detail only at the mean value (r^). For the exact and AFM solutions, we have 
respectively 

{n\r^\n) 
(HO;?i|r^|HO;«) 

In contrast with the previous case, all observables are very well reproduced. This is also 
the case for \, while the overlap |(?i|HO;?i)p tends towards very small values in this 
limit. This is due to the fact that some observables are not very sensitive to the details of the 
wavefunctions and that the sizes of exact and AFM states stay similar. 



Table V. Various observables A computed with the AFM and P{r) = r^. 
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4 n + T 



n + 



4/3 



: 4.214 n + 



4/3 



4/3 



, (3-33) 
(3-34) 



The mean value {¥lO;n\p^ + r\¥lO;n) being analytically computable for an arbitrary 
value of the scale factor A, the relations between the bounds can also be checked. We 
can verify that ^(vo) = ^(2?i + 3/2)^/3 >E = -a„ and E{vo) = ^(2« + 3/2)^/3 > 
E* (vo) = {n + 3/4)2/3 _^ ^ _^ 6y/^r{n + 3/2). Moreover, for the ground state, we have 

£(vo) = IS > £*(vo) = |g + ^ > (V*) = > £ = -ao. 

As in the previous case, the behavior of observables computed with the AFM is similar 
for values of Z = 0, 1,2. We do not expect strong deviations for larger values of /. Some 
typical results are presented in Table |Vll Agreement between AFM and exact results are 
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very good, much better than for the previous case. This is expected, since eigenstates for 
a quadratic potential are closer to eigenstates for a linear potential than eigenstates for a 
Coulomb potential. The quantum number dependence of the scaling parameter X corrects, 
much better than ^ , the difference between the shapes of AFM and exact eigenstates. 

Table VI. Ratios between the AFM results (energies and (r)) with P[r) — and the exact results, for several 
quantum number sets in, I). 

I n=0 n=\ n=1 n=3 n=4 n=5 

1.059 1.066 1.067 1.067 1.067 1.067 

1 1.036 1.055 1.060 1.063 1.064 1.065 

2 1.026 1.046 1.054 1.058 1.061 1.062 

{llO;n\r\UO;n) / {n\r\n) 

0.976 0.964 0.962 0.962 0.961 0.961 

1 0.985 0.972 0.968 0.965 0.964 0.963 

2 0.990 0.978 0.972 0.969 0.967 0.966 



3.3.3. General considerations 

Within the auxiliary field method, if the problem studied is analytically manageable, an 
eigenstate can then be determined with the same computational effort for any set of quan- 
tum numbers («,/). So this method is very useful if one is interested in obtaining analytical 
information about the whole spectra, wavefunctions, and observables of a Hamiltonian with- 
out necessarily searching a very high accuracy. The selection of the potential P(r) seems 
crucial to obtain good results for the eigenstates. For a hnear potential, P{r) = is clearly 
the best choice. In Ref . |53] ). it is shown that P(r) = —l/r provides much better results 
for an exponential potential, while the choice of P(r) is not so crucial for the logarithmic 
interaction. 

An eigenvalue equation can also be solved within the variational method by expanding 
trial states in terms of special basis states. The correct asymptotic tail can be well repro- 
duced if the basis states are well chosen. With this method, a matrix representation of the 
Hamiltonian is obtained and the solutions are computed by diagonalizing this matrix: M up- 
per bounds on the energies are determined with the corresponding M states, where M is the 
order of the matrix. A very good accuracy is possible if M is large enough. Even for M = 1 , 
the accuracy can be better than the one provided by the AFM (see Sect. 16.4.3] ). However, if 
one is interested in closed-form results, the matrix elements must have an analytical expres- 
sion and the number of computed states M must be limited to 4. But even for M = 2, the 
eigenvalues can have very complicated expression, not usable in practice. So, the variational 
method can only provide, at best, a very limited number of eigenvalues and eigenstates with 
an analytical form. 

The WKB method is also a popular method to solve eigenvalue equations.'^'DSJ'in) in 
principle, it is only valid for high values of the radial quantum number n, but it can some- 
times yield very good results for low-lying states.^^' Wavefunctions are not necessary to 
compute some observable^^ but they can be determined for arbitrary value of n with this 
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method. An advantage is that their asymptotic behavior can be correct but, unfortunately, 
the WKB method is mainly manageable for S-states. Indeed, for / 7^ 0, the interaction V{r) 
must be supplemented by the centrifugal potential, which complicates greatly the integrals 
to compute. Moreover, these wavefunctions are piecewise-defined whose different parts 
must be connected properly at the turning points. So they are not very practical to use. 

3.4. Improving the eigenenergies 

3.4.1. Modifications of the principal quantum number 

We have a priori no idea of the dependence of the exact result e{X;n,l) on the param- 
eter A and on the quantum numbers («,/). The AFM gives approximate answers: (13101) or 
( I319I ). Moreover, a simple glance at both formulae convinces us that they have the same 
dependence in A, but they differ in the dependence on («, Z) only through a different expres- 
sion for the principal quantum number Q which is a remnant of the function P(r) used in 
the AFM. 

Owing to the fact that AFM gives the exact result for A = — 1 and A = 2, a continuity 
argument suggests that the exact result should not be too different from the previous formu- 
lae, but with a modified expression for the principal quantum number. Thus we propose for 
the eigenvalues of the Schrodinger equation with power-law potential (13-31 ) the following 
prescription 

£AFM(A;n,0 = ^A\x\2TI2-^Q{X■n,l)i^. (3-35) 

The choice of the function Q{X;n,l) is just a matter of guess but we require that it is a 
continuous function of A and that it coincides with Qc{n,l) for A = —1 and with QnoinJ) 
for A = 2. Besides, there is some freedom. Asymptotically for large I, both Qc and Qho 
are proportional to / with a slope equal to unity. Owing to the fact that Qc and Qho lead to 
lower and upper bounds, a unity slope is maintained asymptotically whatever the value of A. 
Asymptotically for large n, both Qc and Qho are proportional to n and we maintain also this 
characteristic behavior for other values of A. In this work we propose two forms of Q{X;n,l) 
which improve largely the results as compared to the Qc and Qho values. Although in a 
context slightly different, R.L. Hall has considered an improvement of the results based on 
a function /'„,/(A) which was shown as monotone increasing. His prescription to improve 
the results (see Ref . 1571 ) was different from ours, but the spirit is the same. 
The first prescription depends on two free functions b, g and is chosen as 

Qs{X\n,l)=b{X)n + l + g{X). (3-36) 

The preceding conditions impose b{2) = 2, g{2) =3/2 and ^(—1) = 1 = g(— 1)- The simple 
form Qs is linear (it is called sometimes the "linear approximation") both in n and / and, 
thus, looks as simple as Qc and Qho but allows much better results, as will be shown below. 
The reader may use it if he needs a simple formulation giving an average relative precision 
of order of 10^^. 

The second prescription depends on five free functions a, c, d, e, f and reads 

a{Xy + f + c{X)nl + d{X)n + e{X)l + f{X) 
Q,{^;n,l) = ^^y-^ . (3-37) 
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The boundary conditions are in this case {a,c,d,e,f){?i = 2) = (2,3,7/2,5/2,3/2) and 
{a,c,d,e,f){X = —1) = (1,2,2,2,1). The form Qq is a rational fraction with a quadratic 
numerator (it is called sometimes the "quadratic approximation"). This form is more so- 
phisticated that Qs but it gives results with a very high accuracy of order of 10"^. Some 
examples are given in the next section. Other forms have been tested but formulae (13-361 ) 
and (13-371 ) appeared very convenient. Contrary to the case Q = Qho^ for which £afm is a 
upper bound and the case Q = Qc for which it is a lower bound, we have no certainty con- 
cerning the position of the AFM energy as compared to the exact energy in the cases Q = Qs 
Q = Qq \ but in any case, we will show that the approximate value is always very close to 
the exact value. 

Except the very special case of the linear potential, which is studied in detail above, 
no analytical expression for the exact energy e{X;n,l) exists (even in S-state) for potentials 
with A / —1 or 2. Therefore to compare our AFM results to the exact ones, one needs 
very accurate eigenvalues £num('^;'^,0 of Hamiltonian (13-31 ). For that purpose, we rely on 
a very efficient method, called the Lagrange mesh method."^ Without difficulty, one can 
get an accuracy of 10^^ much higher than what could be expected from AFM. Thus, in the 
following, we identify fully the exact value e and the numerical value £num- 

The ultimate aim is to determine the approximate AFM energies £afm that stick as close 
as possible to the exact values. Two kinds of approximations, denoted generically {a), are 
considered; in both of them £afm given by (13-351 ) but in the linear approximation ((a) = 
s) the principal quantum number is given by (13-361 ) while in the quadratic approximation 
((a) = q) this number results from (13-371 ). The free parameters of each approximation are 
denoted cy(A). There are two of them a = {b,g) for the linear approximation and five of 
them a = {a,c,d,e,f) for the quadratic approximation. 

3.4.2. Comparisons to numerical results 

We now search for an analytical formula giving CJ(A). In order to do that, we proceed 
in the following way. In a first step, we choose a sample {Ai = — l,A2,...,Ap = 2} of 
p values of A in the considered domain — 1 < A < 2. We found that = 10 is a good 
compromise between the researched quality and the numerical effort. For each value A, 
belonging to the sample, we calculate, using the Lagrange mesh method,^"*' an array of 
"exact" eigenvalues £(A,;«,/) for «max + 1 values of the radial quantum number n (n = 
0, . . • ,«max) and for Z^ax + 1 values of the orbital quantum number / (/ = 0, . . . , /max)- For the 
precise study that we want, it is enough to choose «max = 8 and /max = 5. The AFM results 
depend on the free parameters a (A,) = a,. The game is to find the values of a, which make 
the AFM energies the closest to the exact ones. 

A tool to do that consists in the minimization of the chi-square quantity 

^(^') = 7 -TTT^— — rTl£(£(A,-;«,/)-£AFM(a,;«,/))^ (3-38) 

C^max I ) ('max i '-) g 1=0 

Other choices are possible but we find this one very convenient. In some very specific 
situations, this function may lead to bad results and we find convenient to introduce an 



25 



alternative minimization based on the difference quantity 



"max 'mux 



^{(^i) = 7 —TTfj -TT £ £ |£(A,-;?i,Z) -£AFM(a;;?l,/)| . (3-39) 

The minimization of one of these functions provides the values a, which will be denoted 
anum(^) because they rely on a numerical prescription. Generally the chi-square and the 
difference procedures give very similar results, but sometimes a glance at the results may 
lead us to favor one against the other. Of course for A = — 1 and A = 2, we know that the 
values Onumi^i) given by the values discussed earlier lead to the exact results so that the 
minimization procedure is useless and both xi'^i) ^^id ^{^i) vanish. Over the full range of 
A sample (10 values), we found that the maximal value of the chi-square is 4.4 10^'* and 
2.2 10^^ for the linear and quadratic approximations respectively, while the accuracy given 
by A is 1.33% and 0.07% respectively. Thus, one sees that the very simple prescription 
(a) = s gives already a relative accuracy of the order of 10^^ over the whole spectrum 
(54 values), while the more sophisticated prescription (a) = q gives practically the exact 
results. Several formulae exist in literature for approximate values of e{X;n,l) resulting 
from various techniques. Among them, one of the most precise is reported in Ref. [5TI) with 
the use of a WKB approach. It must be stressed that our method, even in its crudest version, 
looks much simpler and gives much better results. Thus, we are confident that our formulae 
are really an improvement compared to the previous existing ones. 

The second step of our study is now the guess of continuous functions CJ(A) which 
stick as much as possible to the values anum(A) for the particular values A = A,-. The choice 
is huge because we have freedom on the form of the function, and then on the parameters 
entering it. Since we research simplicity above all, we find that an hyperbola form for all 
the parameters is very well suited. The fit of the 3 parameters entering the definition of an 
hyperbola is done using again a chi-square function x' based on the numerical value a, and 
the proposed value a (A;). 

We summarize below the definite choice for the various parameters. For the linear 
approximation (13-361) . one gets 

,,,, 41A + 86 5A + 17 

For the quadratic approximation (I3-37I ). the corresponding functions look like 
,., 43A-f82 ... 171A+675 136A+330 

^^^ = 1^1^' ^W=mTI^' 'W = ^5aTT2^' ^'-''^ 

109A + 517 ^ 225A+729 

30A+234 ' 94A + 598 ■ 

One verifies that for A = —1, 2 the exact values are recovered. We choose to use integer 
numbers below 1000 in order to match at best the real numbers yielded by the minimization 
procedure. 

As a matter of test, we present in detail the special case of the linear potential (A = 1) 
studied in Sect. 13.31 £(1;«,0) is then just proportional to the opposite of the {n + l)th 
zero of the Airy function. An asymptotic solution of Hamiltonian (13-31 ) valid for large n 
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is given by (|3^ . Equating (|3^ and (l336l) . we find = 7r/\/3 1.814 and g{\) = 
\/3n/4 1.360. These values are the exact ones to describe the spectrum for / = and 
large n, but not necessarily for / / or low n values. However the values coming from our 
fit (13401) are ^^(l) = 127/71 ^ 1.789 and g{\) = 22/16 = 1.375, very close to these new 
values. Assuming that this procedure can be extended for any values of /, the hyperbola can 
be determined completely by using the exact values ^(—1) = 1, g(— 1) = l^ b{2) = 2, and 
g{2) = 3/2. The coefficients are then given by 



V371, 

(4z-18)A + (18-2z) 
(3z-15)A + (21-3z)' 



(7z-36)A + (36-5z) 
^^^)=(6z-32)A + (40-6z)- ^'"''^ 

These relations could look quite complicated but it is worth noting the symmetries existing 
between the absolute values of the numbers present in the formula. When a number appears 
twice in a coefficient of the numerator (denominator), it is the arithmetic mean of two other 
corresponding numbers in the denominator (numerator). A number which appears twice in 
g{X) is the double than the corresponding number mb{X). 




Fig. 2. Coefficients b{X) and g(A) from J3 36b as a function of A, obtained from various approximations: 
numerical fit minimizing the function x (solid line), hyperbolae from | |3-42I ( (short-dashed line), hyperbolae 
from ( l3-40b (dashed-dotted line). These last curves are nearly indistinguishable from the solid lines. 

In Fig. 121 we give an idea of how good is the fit of CJ(A) with hyperbolae. We show 
only the plots for the functions b{X) and g{X) for the linear approximation. The results look 
similar for the functions a(A) of the quadratic approximation. The solid lines represent the 
numerical values anum('^i) obtained by a chi-square procedure; in a sense they are the best 
that we can find. The short-dashed line represent an hyperbola determined analytically; the 
corresponding functions are given by (l3-42! ). The dashed-dotted line represents the hyper- 
bola which sticks as close as possible to the numerical curve; the corresponding functions 
are given by (l3-40! ). This last prescription appears to be very good. 

In Table IVIII we present the lower part of the spectrum (?imax = ^max = 3) obtained 
for different types of approximations and compare the various results to the exact ones. It 



27 



Table VII. Values of £(A = l;n,l) for various approximations. For each set (n.l), the first line (bold case) 
is the exact result obtained by numerical integration, the other lines are obtained with the AFM approxi- 
mation ( l3-35t . For the principal quantum number 2(1;"! 0' 'he second line uses the prescription i3-42t . 
the third line the prescription l l3-36t with l l3-4Qt (linear approximation), and the fourth line the prescrip- 
tion i3-37\ with i3-41t (quadratic approximation). 
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is clear that all the considered approximations give very good results. As expected, the 
approximation with (13-421 ) is very good for / = and large n values, but it becomes worse 
and worse (although not so bad) as / increases. The linear approximation (13-361 ) is very 
good everywhere and should be considered as a good compromise between simplicity and 
accuracy while the quadratic approximation is really excellent everywhere and must be used 
if we want analytical accurate expression of the exact energies. 

Other approximate forms for the eigenvalues of Hamiltonian (13-31 ) exist and give good 
results,E2}'EH) but they are not considered here because they do not give the exact result for 
the cases A = — 1 and A = 2. 

Obviously the functions o(X) depend on the points used for the fit but also on the par- 
ticular choice of the functions % and A . Other definitions-relative error instead of absolute 
error or different summations on quantum numbers-would have given other numbers. Nev- 
ertheless, the quality of the results tends to prove that they are very close to the best possible 
ones. 

§4. Schrodinger equation with an arbitrary potential 

4.1. Scaling laws 

Let E{m,G,a) be the eigenvalues of a Schrodinger equation corresponding to a system 
of reduced mass m subject to a potential of intensity G and characteristic inverse length a. 
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The scaling law gives the relationship between E{m,G,a) and E{m' ,G' ,a'). Let us start 
from the coiTesponding Schrodinger equations 



~Ar + GV{ar)-E{m,G,a) 
2m 



2m 



-Ar + GV{ar)-E{m',G,a) 



•f^'(r) =0. 



(4-1) 
(4-2) 



The important point is that it is the same function V{x) which appears in both equations. In 
(14-21) . let us make the change of variables r = ax and multiply it by %■ Now, we choose the 
arbitrary parameters a and % in order to fulfill the conditions xjijda}) = \/m and aa' = a. 
In other words, we impose the following values 



a 



X 



m / a 



a m \a 

With these values, (14-21 ) can be recast into the form 



1 ,m' f a\'^ , , m' f , , , „ 

A, + G'-(-) V{ax)- — [-] E{m',G',a'] 
2m m \a J m \a J 



(4-3) 



0. 



(4-4) 



Equation (14-11 ) can be recovered, provided one makes the identification G 
G' {m' /m){a/a'Y and a similar relation for the energies. 

The scaling law is thus expressed in its most general form as 



Eim,G,a) = — [- ] E \ m',G' = G— { - 
m \a' J \ m' \ a 



(4-5) 



In fact, it is always possible to define the function V{x) so that a' = \ (for example in 
defining a new position operator as a'r). In what follows, and without loss of generahty, we 
will apply the scaling law for energies under the form 



E{m,G,a) 



m'a^ ^ ( . . mG , 
E[ m',G' = ^^A 



m 



ma'^ 



(4-6) 



This equality is very powerful since it is valid for the eigenvalues of a nonrelativistic 
Schrodinger equation with an arbitrary central potential. It allows to express the energy 
in terms of a dimensionless quantity and some dimensioned factors, as we will see below. 
Moreover it is possible to give to the mass m' any particular convenient value (for example 
m' = 1) so that E{m,G,a) is expressed through the energies of a very simple reduced equa- 
tion depending on a single parameter G'. Instead, one can choose to impose a given value 
of the intensity G' (for example G' = 1) so that E{m,G,a) is expressed through the energies 
of a very simple reduced equation depending on a single parameter m'. 

To summarize the previous discussion, the scaling laws are very powerful since they 
allow to express the most general form of the energy E{m,G,a) in terms of the eigenvalues 
of a reduced equation in which two parameters among the three available can be set at a 
given fixed value. The chosen parameters as well as their values are determined in the most 
convenient way. 
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4.2. General formulae for energies and mean radius 

In Sect. [3l we remarked that, using either P{r) = r^ or P{r) = — 1/r, we obtain, for 
power-law potentials, AFM energies which have exactly the same form but differ only by the 
expression of the principal quantum number Q. In other words, this means that ei^YM{^\n,l) 
appears always in the form F{X;Q{n, I)) with a universal form of the function F. The only 
difference is that, in F, we must use QnoinJ) for P{r) = r^ and Qc{n,l) for P(r) = —1/r. 
This is obvious from the comparison of (13- 101 ) and (13- 191 ) (or even from (13-121 ) and (13-211 )). 

This property allowed us to modify the expression of this number Q in order to im- 
prove drastically the quality of the results. This universal property seems also valid for the 
expression of the mean radius, as you can see from (13-91 ) and (13-181 ) (or even for (13-111 ) and 
(13^ ). 

We will show in this part that these invariance properties are even more general and are 
valid not only for a power-law potential but whatever the potential V{r) to be used.^^' This 
property is even stronger since it persists for any power-law type for the starting function 
P{r) and not only for a quadratic or Coulomb form. Let us demonstrate this fundamental 
property. 

For beginning, let us recall that the AFM energies for the Hamiltonian (13-11 ) are given 
by (see (EH) and (13351) ) 



A 



£AFM(m,a,A;n,/) = ^A{a\l\)jh ('^(^i^') . (4-7) 
2A \ m J 

If we ask that EAFM{tn,a,X;n,l) = E{m,a,X;n,l), we can consider that this formula gives 
the definition of the principal quantum number Q{X;n,l) for all power-law potentials. The 
exact form is known only in a limited number of cases. Obviously, Q{—l;n,l) = n + l +1 
and Q{2;n,l) = 2« + Z + 3/2. Using the notations of Appendix iGl we can write 

Q{l-n,0)=2[-^f"^j=n + ^7r. (4-8) 



In the other cases, we showed that (14-71 ) was able to give the results with a relative accuracy 
better than 10^^ if a good choice is made for the value of Q. Following the philosophy of 
Sect. 12.51 we consider that expression (14-71 ) represents the exact expression of our problem 
and pursue the general AFM method with the starting function P{r) = P^^\r) = sgn(A)r'''. 

In this case, K{r) = V'{r)r^^^ /\X \ and the extremization condition (12-241 ) leads to the 
transcendental equation 

m\X\ ■ ^^"^^ 

It seems that ro is really A -dependent, but replacing K{r) by its expression in terms of the 
potential, this equation turns out to be 

rlv\ro) = ^^^. (4-10) 
m 

For a given potential, the inverse function &{x), 

&{x^V'{x))=x or &{xYv'{&{x)) =x, (4-11) 
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is obviously universal and one has ro = &{Q^ /m). 

The value of the AFM energy follows from (12-251) and reads 



?APM = ^^ + nro). (4-12) 



Thus the function 

2&{x 



^{x) = 7^7;^+y{&{x)) (4-13) 



is also universal and one has 

£afm = =^(GV'«)- (4-14) 

Since it is always possible to choose the same mass for the genuine Schrodinger equa- 
tion with V{r) and the Schrodinger equation with the AFM potential P^^\r), the previous 
demonstration shows that both the mean radius and the eigenenergy depends on A through 
the principal quantum number Q{X;n,l) only. 

We end up with the very important conclusion that can be stated as a theorem (expressed 
below for the energy but also valid for the mean radius) 

If, in the expression E{Q{X;n,l)) of the approximate energies resulting from the 
AFM with P^'^\r), one makes the substitution Q{X;n,l) — )• Q{ri;n,l) (so that 
E{Q{X;n,l)) — ?• E[Q[r\;n,l)) with the same functional form for E), one obtains 
the approximate eigenenergies resulting from the AFM with P^^\r). 
In a sense, as long as we use a power-law potential P^^\r) as starting potential, there is 
a universality of the approximate AFM expression of the eigenvalue, depending only on 
the potential V{r). The only remainder of the particular chosen potential p(^)(r) is the 
expression of Q(X;n,l), as given by (13-361 ) for instance. This result holds whatever the form 
chosen for the potential V{r), even if we are unable to obtain analytical expressions for the 
AFM approximation. 

This property was first remarked for the particular case of a power-law potential V{r) = 
r^ switching from the harmonic oscillator (A = 2) to the Coulomb potential (A = —1). We 
proved here that it is in fact totally general. It is probably related to a well-known property 
in classical mechanics: one can pass from the motion of a harmonic oscillator to the Kepler 
motion by a canonical transformation. This universality property will be used extensively 
in the following applications; it allows us to choose the form P^'^^ir) which is the most 
convenient for our particular situation. Nevertheless, the existence of possible lower or 
upper bounds can only be guaranteed for A = —1 and 2, since formula (14-71 ) is exact only in 
these cases. 

By rewriting (14-101 ) and recalling that the kinetic energy is denoted by P(jp) with Fix) = 

2 

the AFM formula for the energies can be written into the form 

£AFM = r(/?o)+V(ro), (4-15) 

PQ = -. (4-16) 

Po7^'(A))=roV'(ro), (4-17) 

where Q stands for 2(A ; «, Z) to lighten the notations and to stress that we have some freedom 
in the choice of this quantity. Equation (14-161 ) defines a "mean impulsion" p^ from the 
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mean radius ro resulting from (14- 101) . Looking at equation (I4-17I ). one can recognize the 
general form of the virial theorem.^^ ' Finally, (14-151 ) gives the energy as a sum of the kinetic 
energy evaluated at the mean impulsion po and the potential energy evaluated at the mean 
radius rg. This makes clearly appear the physical content of the AFM approximation. Let 
us point out that po is not a linear function of Q, as suggested by (14-161 ). since ro has a 
complicated dependence on Q through (I4-17I ). The approximate eigenstate |vo) is a solution 
of h{vo) = I (2/m) + sgn(A)Vor^ and has a characteristic size depending on mVo = mK{ro) 
(see (IE-31 ). (IF-31 ) and (IG-31 )). This quantity can be easily computed once (14-171 ) is solved. 
The virial theorem applied to Hamiltonian ^(Vq), given by (13-11 ) with a — )• Vo, implies 

that 



Vo 





Vo^ = (vo 


Vo|A|/ 


m 







Vo 



(4-18) 



with Vo = ^(''o) given by (14-91) . The use of definition (14-161 ) with property (12-181 ) gives 



Vo 



pI 



Vo 



Vo 



Vo 



(4-19) 



The last equality is another way to write the fundamental property (12-181 ). These equations, 
as well as the boundary character of the solution, are only applicable when the exact form 
Q{X;n,l) is used. In practice, this occurs when A = —1 or 2, or A = 1 for S-states only. 
Nevertheless, a better accuracy can be obtained by an appropriate choice of Q as shown in 
Sect. 1X41 

In the rest of this section, several interactions will be studied with the AFM. Approxi- 
mations for all states in the spectrum will be given. When the potential allows only a finite 
number of bound states, formulae for corresponding critical constants will be also presented. 

4.3. Square root potentials 

The square root potential that we study in this section is fundamental for the under- 
standing of hybrid mesons, which are exotic mesons actively researched nowadays. Indeed, 
there are two possible descriptions of hybrid mesons within constituent models: 

• First, a genuine three-body object made of a quark, an antiquark and a constituent 
gluon. 

• Second, a two-body object made of a quark and an antiquark in the potential due to the 
gluon field in an excited state. 

It has been shown that these two pictures of the same object are, to a large extent, equiva- 
lent.'^''^''^ In this part, we focus only on the second aspect of the description. 

In general, the string energy, and, therefore, the potential energy between the static 
quark and antiquark in the excited gluon field is given b}'^-'^ 



V{r) 



(4-20) 



where a is the usual string tension while = InaK + C is a term exhibiting the string 
excitation number K and a constant C. These values depend on the model adopted: a pure 
string theory^-'' or a more phenomenological approach.^^'-^^' For the study of heavy hybrid 
mesons, it is therefore very interesting to calculate the eigenenergies of the Schrodinger 
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equation governed by the Hamiltonian 

„2 

// = f- + Va^r2 + Z72, (4-21) 
2m 

where m is the reduced mass and where the parameters a and b depend on the model adopted, 
which must be confronted to experimental results ultimately. 

In this section we give an analytical approximate expression for the eigenvalues of 
Hamiltonian (I4-21I) . In Ref. |65]) . those results have been exploited to study the mass of a 
heavy hybrid meson in the picture of an excited color field. With this mass formula, it was 
possible to predict the general behavior of the mass spectrum as a function of the quantum 
numbers of the system and search for possible towers of states. These results can then be 
used as a guide for experimentalists. 

4.3. 1 . Expression of the energies 

Using scaling laws (I4-6I ). dimensionless variables e and j3 can be defined so that the 
general energy E{m,a,b) is easily calculated 

/2a^ \ o / m \2/3 

E{m,a,b) = ^—j e{P), with P=b\^) . (4-22) 

e (j8 ) is an eigenvalue of the dimensionless reduced Hamiltonian 

r,2 

h = ^ + \/r2 + j8. (4-23) 

Let us apply the AFM in order to find approximate expressions for the eigenvalues of 
the Hamiltonian (14-231) . Equation (14-101) gives immediately 

4rl-Q\l-pQ' = 0. (4-24) 

A simple change of variable x = 2r^ / Q^^^ plus the definition of the parameter 

16)3 



32^/3 



(4-25) 



allows to transform the transcendental equation in the fourth order reduced equation 4x^ — 
8x — 37 = whose solution is given in terms of the G function discussed in Sect. IB. 21 This 
allows the determination of Cafm by (14-121) . 



eAFM(j8;?i,/) = 2y^ 



Gl{Y) + 



G-{Y) 



(4-26) 



with Y = Y{p;n,l) given by (14-251 ). The problem is entirely solved. 

Equation (14-261 ) is complicated but quite accurate. In order to get a better insight into 
this formula, it is interesting to calculate the limits when Y ^ I and Y <^ I because they 
allow an easier comparison with experimental data. This technical point was treated in 
Ref.f 
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4.3.2. Comparisons to numerical results 

As we discussed in Sect. 14.21 it is possible to improve drastically the results given 
by (14-261) if, instead of using the natural forms Q = Quo or 2 = Gc> we switch to a more 
convenient expression. The only relevant dynamical parameter in this case being j8, it is 
natural to propose an expression depending on this parameter: 2(j3 ;«,/). For power-law 
potentials we proposed a "linear" (in terms of n,l) approximation Qg and a "quadratic" 
approximation Qq. For our special potential, one can adopt the same prescriptions, or choose 
a new better one, if one wishes. Since our aim is not to obtain the most accurate possible 
analytic expressions at the price of complicated formulae, but rather to show that the AFM 
is able to give already very good results even with very simple expressions, we adopt here a 
linear expression, so that we choose Q(P;n,l) under the form 



Q{l5;n,l)=A{p)n + l + C{p). (4-27) 

We allow the j3 value to vary between j3 = for which the potential is purely linear and 
j8 = oo for which it is harmonic. We are again in a situation for which Q = Quo gives an 
upper bound while Q = Qc provides a lower bound. 

The procedure used is similar to the one presented in Sect. 13. 41 and is described in detail 
in Ref. I65l). In order to obtain functions which are as simple as possible, continuous in j8, 
and which reproduce at best the exact values, we choose hyperbolic forms for A(j3) and 
C(j3). Explicitly, we find 

8^+102 30^+53 

The integers appearing in A and C are rounded numbers whose magnitude is chosen in 
order to not exceed too much 100. The A and C functions have been constrained to exhibit 
the right behavior A — > 2 and C — )• 3/2 for very large values of j3. formulae (14-281 ) give 
A(0) = 102/57 « 1.789 and C(0) = 53/39 ^ 1.359. These values are such that A(0) ^ 
7r/\/3 f« 1.814 and C(0) \/3n/4 « 1.360, as expected from the results of a nonrelativistic 
linear potential. 

Our results are exact for j3 — )• oo, and the error is maximal for small values of j8. But, 
over the whole range of j3 values, the results given by our analytical expression can be 
considered as excellent. Just to exhibit a quantitative comparison, we report in Table lVIlTl the 
exact e{p;n,l) and approximate £AFM(i3;«,/) values obtained for j8 = 1, a value for which 
the corresponding potential is neither well approximated by a linear one nor a harmonic 
one. As can be seen, our approximate expressions are better than 1 % for any value of n and 
/ quantum numbers. Such a good description is general and valid whatever the parameter j8 
chosen. 

The upper bounds obtained with P{r) = are far better than the lower bounds com- 
puted with P{r) = —\/r. This is expected since the potential \/ cP-r^ + is closer to a 
harmonic interaction than to a Coulomb one. Better lower bounds could be obtained with 
P{r) = r. But, the exact form of Q is not known for this potential, except for Z = for which 
Q is given by (|4^ . With the approximate form Q = {n/\/?))n + 1 + ^/Zn/A^^ we have 
checked that results obtained are good but the variational character cannot be guaranteed. 
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Table VIII. Comparison between the exact values £{fi\n,l) (2nd line in bold) and analytical approximate ex- 
pressions £f^fyi{fi;n, I) for the eigenvalues of Hamiltonian l l4-23b with /3 = 1. For each set the exact 
result is obtained by numerical integration. 3rd line: approximate results are given by l l4-26t with i4-25t , 
| |4-27I ( and ( l4-28t : 1st line: upper bounds obtained with Q = Qho\ 4th line: lower bounds obtained with 

Q = Qc- 
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2.99541 


3.90193 


4.72059 


5.47723 


6 18692 




2.94841 


3.77899 


4.53310 


5.23246 


5.88996 




2.95032 


3.1161^ 


4.52783 


5.22459 


5.87970 




2.75000 


3.23240 


3.68492 


4.11355 


4.52250 


3 


3.46197 


4.32027 


5.10556 


5.83725 


6.52732 




3.41419 


4.19405 


4.91307 


5.58628 


6.22329 




3.41969 


4.20097 


4.91998 


5.59242 


6.22821 




3.23240 


3.68492 


4.11355 


4.52250 


4.91485 


4 


3.90193 


4.72059 


5.47723 


6.18692 


6.85935 




3.85430 


4.59335 


5.28251 


5.93264 


6.55111 




3.86189 


4.60620 


5.29790 


5.94903 


6.56756 




3.68492 


4.11355 


4.52250 


4.91485 


5.29295 



4.4. Exponential-type potentials 

In this section, we apply the AFM to find approximate closed analytical formulae for 
central potentials of exponential form, that is — ae^^'^'')'' where a, j8, T] are positive real 
number. We thus start with the following Hamiltonian 

H = ^-ac-''P'-^\ (4-29) 
2m 

Two particular cases are of current use in physics: 

• The first one is the Gaussian potential (rj = 2) which was intensively used in molecular 
and nuclear physics because it allows very often an exact analytical expression of 
various matrix elements. 

• The second one is the pure exponential (t] = 1), used also in nuclear physics and for 
which exact analytical expression for the energies in S-state are known. 
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4.4. 1. Expression of the energies 

Using the scaling laws, the eigenenergies E{m,a,l5 ,ri;n,l) can be written 



E{m,a,p,r];n,l) = —e{ri,g;n,l), 
where £{ri,g;n,l) is en eigenvalue of the dimensionless Hamiltonian 



with 



h = p -ge 
2ma 



(4-30) 



(4-31) 



(4-32) 



An interesting feature of exponential-type potentials is that they all admit a. finite number of 
bound states that depends on the dimensionless parameter g, ruling the potential depth, and 
defined by (14-321) . There exists thus "critical constants": Potential depths beyond which new 
bound states appear. We refer the reader to Ref. \66h for detailed explanations about how 
to compute critical constants in a given potential. The definition and the properties of the 
Lambert function, that will frequently appear in our calculations, are given in Appendix |Cl 
Some of these potentials have been studied with the AFM in Ref.lFTb. 
The application of (14-101 ) gives 



n+2 / Tlx 

exp(-r,;) = — . 



(4-33) 



Introducing the parameter 



ri+2\7]g 



17/(1+2) 



(4-34) 



the solution of (14-331 ) is obtainable using the properties of the Lambert function 



ri+2 



1/f? 



(4-35) 



Lastly the AFM energy is given by (14-121 ). A straightforward calculation gives the final 
result 



£AFMiTl,g;n,l) 



_^exp( ^Wb(-F„; 



^7 



^7 + 2 , 



n±2 
I) 



l + R±lwo{-Y^] 



(4-36) 



with Yri given by (14-341 ). The energy £afm must be negative, which implies that W{—Yrj) > 
—2/(77 +2)- This is only possible for the branch Wq with a negative argument, which is 
the case since — F,j < 0. The quality of this formula for the Gaussian potential (t] = 2) is 
discussed in Ref.|68l) and for the pure exponential (rj = 1) in Refs. I53l).l67b. 
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Let us just mention that eigenenergies of Hamiltonian (14-311 ) for the pure exponential 
(t] = 1) can be analytically computed for Z = only. But even in this case, the expression of 
£ {g) is not very tractable since it is formally defined by the relatioiJ^ 

^2V^(2V?)=0' (4-37) 

where Jp is a Bessel function of the first kind. Consequently, the AFM result is interesting 
since it yields an analytical formula for the energy levels of the pure exponential potential 
that is of simpler use than (14-371 ) for / = and, above all, that remains valid for arbitrary n 
and / quantum numbers. 

4.4.2. Critical constants 

By definition, the critical constants grj-„i of the potential are such that SAFuiV^S = 
gr];ni',n,l) = 0, with £ given by (14-361 ). Of course this is the AFM approximation of these 
critical constants, but according to this last equation, one finds that the energy vanishes for 
y,j = 2/{r\ +2)e^^/(''+^), a value which is lower than 1/e as required by (14-351 ). Equiva- 
lently (equating this particular value of with the general value (14- 341 )). one can say that 
the critical constants are given by 

gnM=[-Y) QM ■ (4-38) 

They are such that, if ^ > gnoio, the potential admits a bound state with the given quantum 
numbers hq, Iq. It is remarkable that our approximation scheme, based on the potential 
P{x) = for which an infinite number of bound states is present, is able to predict that only 
a finite number of bound states will be present in exponential potentials. This is another 
test of the ability of the AFM to analytically reproduce the qualitative features of a given 
eigenvalue problem. Let us note that the AFM does not give any information about the 
optimal power-law potential to determine the form of the number Q{n,l). For that, we can 
rely on direct comparison with numerical solution. It is also very remarkable that whatever 
the power Tj, the critical constants are always proportional to the square of the principal 
quantum number Q. 

The critical constants for the Gaussian potential are discussed in Ref. "68^ and for the 
pure exponential in Ref. 67). Let us just discuss a little bit this last case. Equation (14-371 ) 
leads to a determination of the exact critical constants g„o (t] = 1 is here always assumed 
and is no longer indicated); let us denote them as g*^. Indeed these critical constants are 
such that they correspond to a energy level £(^*o) = 0. We are thus led to the equation 

Jo{2y^=0^g:o = ^-f, (4-39) 

where jn is the (n + 1)* zero of the Bessel function Jq. At large n, these are given by 
jn ^ K{n + 3/4),^ leading to the asymptotic expression 




(4-40) 



This result is qualitatively similar to (14-381) . stating that g„Q <x Q[n,0)'^. 
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We can try to use formula (14-401) to improve the result (I4-38I ). Assuming that 2,,/ = 
bn + l + c with no constraint on parameters b and c, it is easy to see that formulae (14-381 ) 
and (14-401 ) coincide for g„o with b = n/e and c = 37i/ (4e). We can then try a new relation 
to compute the critical constants 

Se.i = [-^n + -l + -j . (4-41) 

This formula is in good agreement with exact results: For n G [0,5] and I £ [0,5], the mini- 
mal, maximal and mean relative errors are respectively 0.03%, 8.7% and 3.5%.'^ 

4.5. Yukawa potential 

Among all the central interactions, the Yukawa potential is widely used in atomic 
physics (effective interaction), nuclear physics (long range behavior of the nucleon-nucleon 
interaction, due to one pion exchange), hadronic physics (screened Coulomb force). It is 
thus interesting to apply the AFM to this very important case. We start with the following 
Hamiltonian 

H = ^-a . (4-42) 

2m r 

4.5.1. Expression of the energies 

Using the scaling laws, the eigenenergies E{m,a,li;n,l) can be written 

E{m,a,p;n,l) = ^e{g;n,l), (4-43) 
where e{g;n,l) is en eigenvalue of the dimensionless Hamiltonian 



h = p^-g—, (4-44) 
r 

with 

S=-p-- (4-45) 

This potentials admits also a finite number of bound states that depends on the dimensionless 
parameter g, ruling the potential depth, and defined by (14-451 ). 

The application of the AFM to this case leads to the transcendental equation (see (14-101 )) 

e-'-Vo(l + ro) = — = r. (4-46) 

8 

As far as we know, there does not exist an analytical expression giving ro in terms of T. Let 
us introduce formally the inverse function Q = t^^ of the function t{u) = e"u{u — 1). By 
definition, we have the properties 

Q[e"u{u-l)] = u = e^''"^Q{u){Q{u)-\). (4-47) 

Thus one has ro = —Q{T). A calculation giving the AFM energy through (14-121 ) leads to 
the expression 

.*PMfe;«,/) = f (4-48) 



e 
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Q. (m) is a multi-valued function composed of 3 monotonic branches (see Fig. |3]l. However, 
we have the constraints that ro > and Eafm < 0. This imposes the condition — \<Q.{T)< 
0, with r > 0. This condition automatically selects the branches, passing through the origin, 
defined in [— exp(l/0)/0^,exp(— 0)0^] and whose image is in [—0,1/0], where is the 
golden ratio. So, there is no ambiguity. It is possible to obtain an approximate AFM energy 
formula in terms of usual simple functions, but it is then necessary to use an approximate 
form for the solution of the transcendental equation (14461 ).^^' Since Q.{u) — m for m ^ 1, 
we find 



limeAFMU;«,0 



422 



(4-49) 



With the choice Q = Qc, this corresponds to the Coulomb energy. It is expected since g 
is equivalent to j8 — ^ 0. 






m 


2/e 




-1 















Fig. 3. Left: the three branches of the function £2. Right: in bold, the function £2{T) relevant to compute 
formula J4-48l >. (j) is the golden ratio. 



Upper bounds can be obtained by choosing Q = Qc = n + l + l in (14-481) . This seems 
quite natural since a Yukawa potential is similar in form to a Coulomb one. It is possible 
to improve the accuracy of the AFM formula by using an appropriate form Q = A{g)n + 
I + C(gj^ but the variational character of the approximation cannot then guaranteed. An 
accurate energy formula has been found in Ref. 70) from a fit of the numerically computed 
energy levels of the Yukawa potential. In our notations, it reads 

s-2A'{Qc + a)^ + 2B'Ql 



with 



-jp^ig-gG-nl)- 



8G;nl 



g-gG;nl+2B'Ql 



Zi + — 

s, 



(4-50) 



(4-51) 



where A' = 1.9875, B' = 1.2464 and a = 0.003951, and where 

Z, =Zo(l+p/ + T/2), Si=So{l + rl + 5f), 
Zo = 0.839908, p= 2.7359, t= 1.6242, 
5o = 1.1335, 7 = 0.019102, 5 = -0.001684. 



(4-52) 



This formula is rather different from ours. However, they coincide at the limit g ^ °° when 
Q = Qc (see ( |4-49l >). Equation (14-511 ) gives the critical constants of the Yukawa potential 
with a relative accuracy around 0.4%. 
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Table IX. Eigenvalues £{g; n, I) for a Yukawa potential with g = 30 as a function of (n, /) sets. First line: exact 
value; Second line: formula l l4-48t with Q = Qc = n + l+l; Third line: formula l l4-50t . A * indicates a non 
real or a non negative value. 



/ 













-196.44 


-31.51 


-5.47 


-0.22 




-195.98 


-29.97 


-2.92 






-196.36 


-31.62 


-5.67 


-0.30 


1 


-30.74 


-4.94 


-0.029 






-29.97 


-2.92 


* 






-30.52 


-4.85 


-0.019 




2 


-3.81 










-2.92 










-3.70 









The quality of the two approximations presented above can be appraised in Tables |IX] 
for a particular value of g. Not all bound states can be found with (I4-48I ). Let us note that 
formula (14-501 ) gives generally the correct number of bound states, but not in all cases.^^ 
Formula (I4-48I ) with Q = Qc gives quite good results for the lowest eigenvalues. The quality 
of the fit for formula (14-501) is better but comparable to the quality of AFM formula (14-481 ) 
with a more sophisticated parameterization of Nevertheless, formula (14-501) is em- 
pirical, while AFM result (14-481 ) is obtained from an explicit analytical resolution of the 
Schrodinger equation. 

4.5.2. Critical constants 

Let us discuss below more deeply the critical constants relative to the Yukawa potential. 
It is easy to find from (14-481) that the T value which cancels the AFM energy is equal to 
r = 2/e. Using the definition (14-461) . one obtains the AFM critical constants for the Yukawa 
potential 

gyut,nl=eQ^- (4-53) 

It is again remarkable that the critical constants for the Yukawa potential are also propor- 
tional to the square of the principal quantum number, a property that seems universal for 
exponential-like potentials. 

Due to the smallness of /and 5 in (14-511 ). we can replace Si by in (I4-51I ). It clearly ap- 
pears that gG-„i oc ^So^/ZqP I + asymptotically. This is precisely the quadratic behavior 
predicted by our analytical results. It is important to stress that formula (14-511) was specially 
designed for fitting the Yukawa critical constants and that it needs several free parameters, 
while our AFM result (14-531 ) follows a general algorithm valid in any circumstance. The 
fact that both give the same asymptotic behavior is again a very strong argument in favor 
of the relevance of the AFM. Introducing a constant which allows to obtain the exact re- 
sult for n = Z = and considering the asymptotic expansion for large values of /, (14-511 ) is 
approximately given by 



gG;ni ~ (1.248?i+ 1.652/ + 1.296)^ 



(4-54) 
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This is not too different from our result (14-531 ) with Q = Qc = n + l+l 

§yuk,nl (Gc) = (V^« + Ve/ + ^/i)', (4-55) 

since ^/E 1.649. This formula is unfortunately not very accurate: For n G [0,4] and / G 
[0,4], the minimal, maximal and mean relative errors are respectively 6%, 42% and 26%. 
The minimal, maximal and mean relative errors are respectively 0.01%, 0.9% and 0.3% with 
formula (14-511 ). The very good quality of these last results is due to the use of a complicated 
formula fitted on exact results. Our formula is simpler and the general behavior is predicted 
by the AFM. Let us note that it is possible to improve the values of AFM critical constants 
for the Yukawa potential using an appropriate form Q =An + Bl + C, with A, B and C fixed, 
and keeping the simplicity of the original formula (l4-53b .l^ 

4.6. Sum of power-law potentials 

In Sect. I2.4.T] we suggested that potentials of type aP{r) + V{r) could be treated with 
the AFM in some occasions. Since a power-law potential is very often a good starting point, 
it is interesting to see whether an analytical AFM solution exists for a second contribution 
which is itself a power-law potential. The resulting total potential is thus the sum of two 
power-law potentials. Those types of potential are frequently used in various domains of 
quantum mechanics. This is the subject of this section, whose calculations are based on 
(1^371) and (l238l) . 

4.6. 1. Expression of the energies 

It is very instructive to start with the potential P(r) = sgn(T])r'' = P^''^\r). In this case, 
the extremization condition (12-371 ) is written (see also (14-91 )) 

a\ri\r^,^^ + riv\ro) = ^, (4-56) 
m 

where Qri is a priori the optimal principal quantum number for P^^\r). The AFM energy 
(12-381 ) is given by 

EAFM = ^ + ^gn{ri)ar^+V{ro). (4-57) 

With this in mind, let us introduce the potential V{r) in the form of a power-law potential 
sgn{X)br^ so that the total potential for our system is given by 

W(r) =sgn(77)ar'l +sgn(A)Z7/. (4-58) 

The above considerations show that the extremization condition and the AFM energy look 
like 

a\7l\r^,+^ + b\Mr^+^ = ^ (4-59) 
m 

Eafm = ^ + sgn(T] )ar^ + +sgn{X)br^ . (4-60) 

Should have we taken instead P{r) = P^^\r) and V{r) = P^'^\r), we would have re- 
covered the same physical problem with the potential W{r) given by (14-581 ). Due to the form 



41 



of basic equations (12-371 ) and (12-381) the only change would have been Q-q ^ Qx- Because 
of the symmetry of the problem there is no reason to prefer the choice Qj^ , as in the previous 
equations, instead of the choice ; moreover in any application, we already saw that it is 
convenient to change the value of the principal quantum number. In consequence, in the 
following, we will use the more neutral notation Q. This result can also be directly obtained 
by using (I4-15I )- (I4-17I ). but it is interesting to illustrate that the AFM depends only on the 
potential P{r) through the principal quantum number Q. 

In general, one does not have an analytical expression for the root of the extremization 
equation (14- 591 ). This is possible only for very specific values of the powers rj and A. 

In this section, we will study such a problem for the most favorable cases, where Tj is 
chosen to give an exact expression for the eigenvalues (in practice rj =2 and T] = — 1) and 
where A is chosen in order to have an analytical result. In the following, it is assumed that 
a > 0, b > and -2 < X <2. 

4.6.2. Solvable potentials 

Let us examine first the case P(r) = r^. The extremization equation is a transcendental 
equation for which an analytical solution does not exist automatically. The only cases for 
which we are sure that an analytical solution exists is when it can be transformed into a 
polynomial of degree less than or equal to 4. In order to investigate this condition, let us 
put {X + 2)/A = p/q {p and q are relatively prime integers), and define the new variable 
X = Tq^^. The corresponding extremization condition is transformed into 

2aX'^ + b\X\XP = —. (4-61) 
m 

All the solvable potentials should verify the conditions < p < 4, 1<^<4. An 
exhaustive research of all the solvable potentials in our domain for X leads to the following, 
non trivial, values of the power X : 

A = -2,-1, -1^1. (4-62) 

We will study in detail the cases: 

• X = —2 because it corresponds to a centrifugal term with a real parameter instead of 
the usual /(/ + 1) term; 

• A = — 1 because it corresponds to a simplified potential for hadronic systems with a 
short range Coulomb potential and a quadratic confinement; 

• A = 1 because this anharmonic potential is sometimes used in molecular physics. 
Now we investigate the case P{r) = — l/r. The same type of argument with X+2 = p/q 

and X = leads to the following polynomial equation 

aX'^ + b\X\XP = —. (4-63) 
m 

An exhaustive Ust of the non trivial solvable potentials is given below 

7 5 3 4 5 2 1 

'■—'•-4--r-r-r-4--r-r'-'- 

Among them we will study: 
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• A = 2 because it corresponds to a potential —a/r-\r br^ which is already studied with 
P(r) = r^. This is the only potential that can be described with either P(r) = or 
P{r) = — 1/r and this property allows very fruitful comparisons. 

• A = 1 because it corresponds to the funnel potential (Coulomb + linear) which is 
widely used in hadron spectroscopy.^^' Finding approximate analytical values for the 
energies corresponding to this potential is thus a very interesting question. To our 
knowledge, such formulae have not been proposed in the literature. 

It could be also interesting to introduce potentials with A < — 2 but with the restriction 
that it is repulsive at the origin (for instance, van der Waals forces or Lennard- Jones type of 
potentials). In this case, the transcendental equation is modified a little bit since A + 2 < 
implying a negative fraction p/q. An exhaustive research of the different values of A gives 
the following results: 

A = -5,-4, -3, -^,-^,-2. (4-65) 

In this list, we will just consider the case: 

• A = — 2 because the corresponding potential, known as the Kratzer potential, exhibits 
its spectrum under an analytical form for all values of radial quantum number n and 
orbital quantum number /. 

4.6.3. Kratzer potential 

The Kratzer potential^ ' is defined in its simpler form as 

2a 

V{r) = ---. (4-66) 

It can mimic the interaction of two atoms in a diatomic molecule. Moreover, it presents 
some interest as a benchmark since it is one of the rare potentials for which ones knows an 
exact analytical expression of the energies valid for any n and / quantum numbers. In fact, 
one can define the Kratzer potential in a more general formulation (see Ref . iTTi l) as 

y(r)=G(l-{). (4-67) 

One recovers the simplest form putting G = a^ and f = 2/a. 

It can be shown (see for example the hint given in Ref. |9l)) that the exact eigenenergies 

read 



E{n,l) 



1/2 + ^(/ + 1/2)2 + 2mG 



2 • 



(4-68) 



Applying the previous AFM expressions to this potential leads to the following mean 
radius (t] = -I, X = -1, a = Gf and b = -G in (14-581 )) 

ro = ^ (4-69) 

mGj 



which, inserted in (14-601) . gives the desired result 

'2[2mG + Q'^y 



^^Hn,l) = -^,;^$l^. (4-70) 
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It is natural to set Q = Qc = n + 1 -\- 1, considering that P{r) = —l/r. We thus remark that 
the approximate value E^^^ presents the correct asymptotic behavior for large n and for large 
/. Just to have an idea of the quality of this approximation, let us calculate the difference 
d between the terms in brackets appearing in the denominators of E and E^^\ It is just a 
matter of simple algebra to find 



5 = (2?i + l)(Z + l/2) 



'1 + 



2mG 



(/ + 1/2)2 



1 



(4-71) 



One always has 5 > so that E{n,l) > E^^\n,l) meaning that the AFM approximation 
is a lower bound. Moreover, for small intensity and/or mass, mG ^ 1, or for large angular 
momentum, Z » 1, the approximate value tends to the exact one and we have more explicitly 

mG 



{2n+\) 



(/ + 1/2) 



< 1. 



(4-72) 



For these limit conditions, the AFM approximation tends towards the exact result. This be- 
havior is easily understandable because, under those conditions, the contribution due to 1 /r 
is predominant as compared to the contribution of l/r^, and both expressions tend towards 
the same exact Coulomb result. 

4.6.4. Quadratic + centrifugal potential 

We consider now the potential (for an attractive centrifugal potential, not all values of 
b are relevanP^)) 



(4-73) 



The incorporation of the term ^hjr^ into the /(/+ l)/r2 term already present in allows 
to get the exact eigenvalue using the same kind of arguments than those developed in the 
harmonic oscillator case.'^ Explicitly, we obtain 




2(2?i + l) + A/(2/ + l)2±8mZ? 



(4-74) 



Using the AFM with P{r) = p- and setting Y = Imb/Q^, the value that extremizes 
this energy comes from a first degree equation and reads 



b{\±Y] 
aY 



1/4 



(4-75) 



Substituting this value into the energy (I4-60I ). one obtains a very simple expression for the 



approximate energy E = 2y/ab{l ±Y)/Y or, alternatively 

E^'^''\m,a,b;n,l) =2^J^yQ^±2mb. 



(4-76) 



This quantity and the corresponding exact one depend on three parameters m, a, b but 
we know that the general scaling law properties allow us to write them in a more pleasant 
form 



E{m,a,b;n,l) 



' — e(B = 2mb;n,l), 
2m 



(4-77) 



44 



where e{^;nj) is an eigenvalue of the reduced Schrodinger equation for which the Hamil- 
tonian depends now on a single dimensionless parameter jS 

h = p^ + r^±^. (4-78) 
The exact eigenvalues of this Hamiltonian are given by 



£ (i3 ; /) = 2(2?i + 1 ) + (2Z + 1 )2 ± 4j8 . (4-79) 

The approximate values immediately come from (14-761) 

e(^^'(j8;n,Z)=27e2^. (4-80) 

It seems natural to choose Q = Qho^ considering that P{r) = r^. In this case, one can check 
that the relative error between e and g^'''^' decreases as l^^ for a fixed value of n and large /; 
it decreases as n^^ for a fixed value of I and large n. This behavior is easily understandable 
because, for large values of the quantum numbers, the contribution due to is predominant 
as compared to the contribution of 1 /r^ and both expressions tends towards the same exact 
harmonic oscillator result. 

Let us assume that j3 ^ 1 and let us choose Q = Qho\ ^ Taylor expansion truncated to 
first order leads to 

£(q^)(i8;n,/)«2G//o±7r-- (4-81) 
Quo 

It is easy to check that this expression can also be obtained by perturbation theory. In 
particular for j8 = 0, one recovers the exact value 2Qho^ as expected. 

4.6.5. Anharmonic potential 

The potential under consideration reads 

V{r)=ar^ + 2br. (4-82) 



The extremization condition writes 



Let us introduce the parameter 



2ari + 2brl = —. (4-83) 
m 



Y = -a[^] (4-84) 



3 \ mb'^ 



1/3 



and the new variable x = ^^^j The previous equation (14-831 ) is put in the simplest 

form Ax^ — 8a; — 37 = 0. This is the reduced quartic equation presented in Sect. IB. 2l whose 
solution is given by x{Y) = (see (IB-71 )). Substituting this value into the energy (14-601 ) 

and making a little algebra leads to the desired approximate energy 

E^^\m,a,b;n,l) = ('^-(^) + ^^F) ) ' ^^'^^^ 
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with Y given by (14- 841 ) . It is natural to use the value Q = Qho^ considering that P(r) = r^. 

As in the previous case, this quantity and the corresponding exact one depend on three 
parameters m, a, b but the general scaling law allows us to write them in terms of a reduced 
quantity depending on a single parameter /3 




'ib^ [3m 



Eim,a,b;n,l) = ^—e\^p = _y j , (4-86) 
where e{^;n,l) is an eigenvalue of the reduced Schrodinger equation for the Hamiltonian 

h = ?^ + 3? + %^r. (4-87) 
The approximate value corresponding to this reduced equation follows from (14-851 ) 



e(-->il5-n,l) = 2l5Y(Gi{Y) + -^-^j, Y = j . (4-88) 

The parameter j8 could also be associated with the quadratic potential, but this less interest- 
ing case is not considered here. 

Let us assume that j8 ^ 1 and let us choose Q = Qho\ ^ Taylor expansion truncated to 
first order leads to 

£(-)(i8;n,Z) « V3Qho + 4^^^^- (4-89) 

In particular for j8 = 0, one recovers the exact value V^Qho^ as it should be. This result 
comes also from the perturbation theory. 

The limit j8 — )• oo is not physically relevant, but it is interesting to consider it in order to 
check the formula. In this limit, we find 

£("")(j8;«,/) = 3{4pQy/^ + O ■ (4-90) 



The dominant term is the result expected for a pure linear potential, as given by (I4-7I ). 
4.6.6. Quadratic + Coulomb potential 

We study now the quadratic + Coulomb potential defined as 

V{r) = ar^--. (4-91) 
r 

The equation giving the mean radius looks like 

2a4 + bro = —. (4-92) 
m 

Let us introduce the parameter 

8^2 /4a \ 1/3 
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1 /3 

and the new variable a; = (^) tq. The equation that leads to the extremization of the 
energy is then 4x^ + 8;c — 3F = 0. This reduced quartic equation is studied in Sect. IB.2I Its 
solution is given by x{Y) = G+{Y) (see (IB-71) ). Substituting this value into the energy (l4-60b . 
one is led, after some manipulations, to the desired approximate energy 



E(^^\m,aMn,l) = l(^ 



1/3 



GliY) G+(y) 



(4-94) 



with Y given by (I4-93I ). As usual with the AFM, the quality of the approximation depends 
on the choice of Q. For instance, the use of 2 = 2c and Q = Qho in the expression for Y 
gives respectively a lower bound and an upper bound on the true energy. 

The quantity (14-941) and the corresponding exact one depend again on three parameters 
m, a, b but the general scaling law allows us to write them in terms of a reduced quantity 
depending on a single dimensionless parameter j8 . One can imagine two formulations de- 
pending on whether j8 is part of the quadratic contribution or of the Coulomb contribution: 

1/6 \ 



2a 



E{m,a,b;n,l) =4^/— £ ( j3 = - 



1 /SAm^f^ 



E{m,a,b;n,l) 



3m 
3mb^ 



16 



T7 /3' = 4 



54m3M 



;n,/ 



1/6 



;n,l 



(4-95) 



(4-96) 



The £ and rj energies are the eigenvalues of the reduced Schrodinger equations for the 
respective Hamiltonians et /j,, : 



h 



16 4 r ' 
16 r 



(4-97) 



(4-98) 



The approximate values corresponding to these reduced Hamiltonians follow from ( 14-941) : 

2 



4 



Gl{Y) G+(F) 



Y=\ ^ 



T](^^)(i8';n,/) 



3j8 



/2 



Y' = {QP') 



(4-99) 



(4-100) 



Let us assume that j3 ^ 1 or j8' ^ 1 (it is equivalent to the formulation V{r) <^ P{r)): 
• A Taylor expansion truncated to first order for the formulation based on the £ form 
leads to 

V3. 



,(qC) 



-Q 



lip 



(4-101) 



In particular for j8 = 0, one recovers the exact value ^/3Q/A ii Q = Qho- 
The Taylor expansion for the formulation based on the r] form gives 



T](^^)(i8';«,/) 



8 9p"'Q 



6^4 



■ + 



3^2 ' 128 

In particular for j8' = 0, one recovers the exact value -^/{3Q^) if 2 = Qc- 



(4-102) 



47 



In all these cases, the approximate formulae resulting from AFM agree with the result of 
perturbation theory. 

The limits j3 — oo and j8' — oo are interesting to consider in order to check the formulae: 

• The £-based formulation gives 

e(^c){l5-n,l) = -^ + 0{l5-^). (4-103) 

If 2 = 2c> the dominant term is the exact result for a Coulomb potential. 

• The Tj -based formulation gives 

r7(^^)(/3';«,/) = ^p'^Q + O^P''^') . (4-104) 

If 2 = Qho^ the dominant term is the exact result for a quadratic potential. 
Let us emphasize the point that both (14-971 ) and (14-981 ) can be related with the scaling 
laws. As a consequence, it can be shown that both the exact eigenvalues and the AFM 
approximate ones fulfill the relation 

It is worth mentioning that analytical solutions of the Schrodinger equation with 
quadratic plus Coulomb potentials are presented in Ref. 73), within the framework of a 
two-dimensional system of two-interacting electrons (— 1/r) in a confining magnetic field 
(r^). Closed-form solutions are found for particular values of magnetic field and spatial con- 
finement length. It is argued that a generalization to a three-dimensional space is possible. 
But, as the corresponding formulae are not given explicitly, a comparison with our results 
is not available. More generally, the Schrodinger equation with a. —a/r + br + cr^ potential 
(particular cases are studied in Sect. 14.6.61 and 14.6.71) is directly linked with the biconfluent 
Heun's equation.l^ 

4.6.7. Funnel potential 

In this section, we study the funnel potential defined as 

b 

V{r)=ar--. (4-106) 
r 

This potential is particularly important and its interest is discussed in Sect. 14.6.21 The ex- 
tremization condition writes in this case 

(92 

brQ + arl = —. (4-107) 
m 



Let us introduce the parameter 




(4-108) 



and the new variable x = a/%''o- The extremization condition for the energy is the solution 



of the equation x +?>x — 2Y = 0. This reduced cubic equation is studied in Sect. IB. II Its 
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solution is given by x{Y) = F^{Y). Inserting this value into the energy (14- 601 ) . one finds the 
expression of the approximate AFM energy 



E^^\m,a,b;n,l) = ^/Zah 



(4-109) 



F+(y)2 f+(F)_ 

with Y given by (14-1081) . Let us mention that another form of this equation can be found 
thanks to the following relation 

Y 2 . . _ 1 



sinh 



4sinh0' 



(4-110) 



with the change of variables Y = sinh (30). 

The quantity (14-1091) and the corresponding exact one depend again on three parameters 
m, a, b but the general scaling law allows us to write them in terms of a reduced quantity 
depending on a single dimensionless parameter j8. In hadronic physics, the dominant inter- 
action between a quark and an antiquark is a confining linear potential.'^ So, although the 
linear potential has no analytical exact solution for all values of the quantum numbers, it is 
also interesting to consider two formulations depending on whether j8 is part of the linear 
contribution or of the Coulomb contribution: 



E(m,a,b:n,l) = 3 ( — 
^ ^ '2m 



1/3 



( 4m^b 



3 \ 1/4 



V 27a 



;n,l 



E(m,a,b;n,l) 



2mb 



f 11a \ 



1/4 



3V 



(4-111) 



(4-112) 



The £ and T] energies are the eigenvalues of the reduced Schrodinger equations for the 
respective Hamiltonians he et h^^ : 



h 



^4/3 



3 ' 3 r 

2 3I/3 



(4-113) 



+p r. 

3 r 



— - — '• (4-114) 
The approximate values corresponding to these reduced Hamiltonian follow from ( 14-1091 ): 



2/3 



77W(j8';n,/) = 32/3/3'' 



Y' : 



F+{Y')^ F+{Y') 



, r = {Qi5') 



A 2 



(4-115) 
(4-116) 



Let us assume that j8 <C 1 or j8' <C 1: 

• The Taylor expansion for the formulation based on the £ form gives 



22/3 



2Q^ 



1/3 



(4-117) 



In particular for j8 = 0, one recovers the value expected for a pure linear potential, as 
given by (14-71 ). 
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• The Taylor expansion for the formulation based on the Tj form gives 

^'^^(/3';«>0«-^ + ^. (4-118) 

In particular for j3' = 0, one recovers the exact value — 3^^^/ {4Q^) for the Coulomb 
potential if 2 = Qc- 

In all these cases, the approximate formulae resulting from AFM agree with the result of 
perturbation theory. 

The limits j3 — )• oo and j3' — )• oo are interesting to consider in order to check the formulae: 

• The e -based formulation gives 



= _^ + o[r'") ■ (4-119) 

The dominant term is the exact result for a Coulomb potential if 2 = Qc- 
The T] -based formulation gives 

r]«(^';n,/) = ^lp''Qy\o[l5''^') . (4-120) 



The dominant term is the result expected for a pure linear potential, as given by (I4-7I ). 
Here again, both (14-1131 ) and (14- 114b can be related with the scaling laws. As a conse- 
quence, it can be shown that both exact eigenvalues and AFM approximate ones fulfill the 
relation 

eW(j8;«,/) = ^T]W(l/j8;«,/). (4-121) 

4.6.8. Comparison with numerical results 

The potentials studied in this section are rather sophisticated and, presumably, the AFM 
should lead to less good results than power-law potentials. Thus, it is interesting to test 
this method in these less favorable cases. A detailed discussion can be found in Ref. l40l) . 
Here we focus our analysis on the funnel potential only, because of its physical importance, 
and rely on the reduced Hamiltonian (14-1 131) . The AFM energies e^^\l5;n,l) are given 
by (14-1151) . For most of the hadronic systems, the physical values of jS vary from to about 
1.5. 

As usual, the simplest prescription to improve harmonic oscillator or Coulomb-like 
results is to adopt a principal quantum number of the form 

Q{p)=b{p)n + l + c{p). (4-122) 

The protocol to determine the best values for the b and c coefficients follows essentially the 
same steps than those explained in Sect. 13.4.21 and is not repeated here. The only difference 
is the form of the functions b{j5) and c(j8) which are better fitted with Gaussian functions 
than hyperbolae. Explicitly, the parameterization retained is 

b{l5) = l+p,exp{-pll5^), c{l5) = \+q,exp{-qll5^). (4-123) 

With this choice, ^(oo) = 1 = c(oo) as it should be since in this case we have a pure coulomb 
potential. For j8 = we are in presence of a pure linear potential and for / = we do know 
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the exact result. With this in mind, we choose b(0) = 7r/\/3 and c(0) = The only 

free parameter is p2 for b{[5) and q2 for c(j8). The parameters retained for our study are: 

P' = 73-^'P^ = O-^l^. ^1 = 4^ - 1, ^2 = 1.245. 

Table X. Eigenvalues £{fi\n,l) of Hamiltonian i4-113t with /3 = 0.5, for some sets («,/). First line: 
EnumCiS;",/) from numerical integration considered as the exact ones; second line: £^^'>{p;n,l) given by 
( l4TT31 l with e(/3) defined by l [4T22t and l l4423l l: third line: e(f)(j8;«, /) given by l ITOSl l with 2(j8) = gc- 



/ 


£(0.5;0,/) 


£(0.5;1,/) 


£(0.5;2,/) 


£(0.5;3,/) 





0.39711 


1.11714 


1.64558 


2.09628 




0.42779 


1.16223 


1.68099 


2.12205 




0.26827 


0.79105 


1.15440 


1.45987 


1 


0.90598 


1.45955 


1.92580 


2.34167 




0.88794 


1.46673 


1.93564 


2.34911 




0.79105 


1.15440 


1.45987 


1.73269 


2 


1.25749 


1.74247 


2.17133 


2.56288 




1.23307 


1.73892 


2.17323 


2.56506 




1.15440 


1.45987 


1.73269 


1.98358 


3 


1.55457 


1.99727 


2.39917 


2.77168 




1.52908 


1.98937 


2.39764 


2.77183 




1.45987 


1.73269 


1.98358 


2.21833 



Allowing a j3 -dependence for the coefficients b and c improves greatly the approximate 
eigenvalues since one gains a factor 10^^ to lO^-' on the chi-square values.l^S) The quality 
of the fits can be appraised for j3 = 0.5 by examining the values of approximate results 
compared with exact ones presented in Table |Xl The values Q = Quo and Q = Qc give 
respectively an upper and a lower bound on the exact energy. However, as can be seen from 
Table 1X1 for the case Q = Qc, these bounds are not very close to the exact value and, thus 
not very interesting (the error is of the order of 10% or more). The value Q = 2(j8) gives 
more accurate results (the error is less than 5% and is often of order of 2-3%), but we have 
no certainty concerning the (anti)variational character of the approximation. The global 
accuracy is less good than for power-law potentials but remains quite acceptable. 

§5. Two-body spinless Salpeter equation 

In this section, we will apply the AFM to the spinless Salpeter Hamiltonian for which 
the kinetic energy has a semirelativistic form 

H = ay'p^+m^ + V{r). (5-1) 

The term "relativistic" is generally not used because such a Hamiltonian is not manifestly 
covariant. The parameter a is equal to 1 for one-body systems (r is then the distance of the 
particle from the origin of the force field) and to 2 for two-body systems with equal masses 
(r is then the distance between the two particles). However, it is very convenient to choose 
it as an arbitrary positive real parameter. We will prefer the notation M instead of E for the 
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eigenvalue of this kind of Hamiltonian because it does include the rest mass of the particles. 

Obviously, there is much less material concerning the properties of eigenvalues for 
semirelativistic equations as compared to the huge bulk of results for nonrelativistic treat- 
ments. Nevertheless the quadratic potential has been studied taking the Fourier transform of 
of the potential.'^' Moreover, working in momentum space, Boukraa and Basdevant were 
able to obtain interesting properties concerning a class of potentials.'^' 

For two-body systems, the presence of two different masses cannot be dealt with a 
reduced mass, so that we need the introduction of two different auxiliary fields to treat the 
kinetic energy term. This situation makes the problem much more complex. It is briefly 
studied in Sect. 15.61 Thus, in the following, we focus on Hamiltonian (15-11) . 



5.1. Scaling laws 

Let us consider the spinless Salpeter equation 



a^/p^ + m^ + GV{ar)-M{m,G,a,a)) ^{r) = 0, (5-2) 



where G and a are two real numbers defining the potential, a being related to its extension 
and G to its intensity. The eigenvalues M depend a priori on the four parameters m,G,a, G. 
Scaling laws allow to obtain the general form of M{m,G,a, a) in terms of the energies of 
another system depending on less free parameters. Let us examine this point. 
Defining p =ar, Jt = p/a, one has 

G I \ 

n^ + im/aY + — Vip) Mim,G,a,o)] Wip/a)=0. (5-3) 

ao ao J 

This is just the spinless Salpeter equation 



'lt'^ + m'^ + G'V{p)-M'{m ,G\\,\)] (p{p)=0, (5-4) 
with 

' r^' ^ ^yr'^ 'n'^ y\ M{m,G,a,o) 

m=—, G= — , M(m,G, 1,1) = . (5-5) 

a ao ao 

We are thus led to the following scaling law for the mass spectrum 

Mim, G, a, o) = ao M { — ,— ,1,1] . (5-6) 
\a ao J 

It means that both one parameter of the potential (here the extension of the potential a) 
and the number of constituents in the system o can always be set equal to 1 in order to 
simplify the computations without loss of generality. The general energy formula will then 
be recovered thanks to (15-61) . 

5.2. Generalities 

To deal with the worrying square root, it is interesting to introduce one auxiliary field 
pL in order to build the AFM Hamiltonian corresponding to (15-11 ) 

H{n) = ?-^ + ^H + V{r). (5-7) 
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The semirelativistic operator is then replaced by a nonrelativistic counterpart with a param- 
eter n to be determined. This technique to get rid of the square root operator is not new 
and has been used previously to study hadronic physics.l^''^''^'l23'E3 Let us remark that, 
as H{lJ.) and H have not the same kinetic part, most of the results about lower and upper 
bounds developed in Sect. 12.41 are not applicable. At this stage, one can imagine two cases 
which are discussed below. 

5.2.1. Potential V{r) is solvable 

Let us assume that the analytical solution of the Schrodinger equation with the potential 
V{r) is known. In particular, this approach is interesting if V{r) = or V(r) = — 1/r. One 
can then write the eigenvalues of as 

M{n) = — + —^ + e{^), (5-8) 

where e{lJ.) is the eigenenergy of the Hamiltonian h{iJ.) = / + V(r). The final spectrum 
is given by M{}JLq), where /Iq is such that ^ai^(/^)|^=^o ~ ^ which implies that 

^ + e'(Ato) = ^. (5-9) 

The prime denotes a derivation with respect to /i. Provided that is analytically known, 
the only difficulty in this procedure is the analytical resolution of this last transcendental 
equation, especially when m / 0. Reporting the value of jUq given by (15-91 ) in the expression 
of the energy (15-81 ) allows to write 

M(Aio) = yMo+(M^'(m))'|^.^„- (5-10) 

From Ref.TT), we know that M{^) > M for any state, and that M{jXij) > M, since it is also 
true for the particular value jLto of the auxiliary field. The AFM yields in this case an upper 
bound. This can be also demonstrated using a simpler version of the AFM which can be 
used only when it is interesting to replace an operator by its square."^^ ' 

Using the Hellmann-Feynman theorem,l22) it has been shown that the value of /Iq is 
given bjISJED 

i"O = |0p' + m2), (5.11) 

where the mean value is computed with an eigenstate of ^(/Xq). An estimation of the differ- 
ence 5 between the exact and approximate eigenenergies is given b}0^ 



5 = aJ {p^ + m^) - a (y' + mA . (5-12) 



At the limit of low mass, (15-11 ) can be written 



2 

H^H"' + -^^ with = a + (5-13) 
2Vp2 
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The contribution A {m) of the mass m to an eigenenergy of the ultrarelativistic Hamiltonian 
//™ (m = 0) appears as a small contribution that can be computed as a perturbation. The 
mean value being taken with an eigenf unction of we hav^ 

^>».^ = f. (514) 

where is the auxiliary field obtained with the AFM applied to the Hamiltonian H™. 
5.2.2. Potential V{r) is not solvable 

If the eigenvalues of h = /}X + V{r) are not known, it is tempting to treat the kinetic 
energy term and the potential term on equal footing, in the spirit of what is done in Sect. 12.21 
introducing the auxiliary field pL to deal with T and the auxiliary field v to deal with V. The 
case of a general function P(r) cannot be solved analytically, but restricting to the important 
power-law potential P{r) = sgn(A)r'*', very interesting expressions can be obtained. 

The eigenenergies e(/i,v) of the nonrelativistic Hamiltonian h{}jL,v) = /pi + vP{r) 
are given by (13-21 ) and (13-351 ) and the AFM mass by (the K and J functions are defined as 
usual) 

2 2 

M(At, V) = ^'(At , V) + + V + ^(■^(^)) - sgn(A)v7(v)^ . (5-15) 

The extremization of M needs to solve two coupled equations with respect to pL and V which 
provide the values /Iq and Vq. Introducing the new variable 

f\X\ \ 2/(^+2) 

XQ={—pioVo\ , (5-16) 

it can be shown that both Hq and Vq can be expressed in terms of xq only. Namely 



Mxo) = ^^m^ + Qf^'^'ho, (5-17) 
Vo(^o)=^(!2f' +'VV^)- (5-18) 

Reporting these values into the definition (15-161 ) allows the determination of xq through the 
transcendental equation 



axi'-'^/' = |A| V-2 + g2A/(A.2)^^^^g2/(A.2)/^)_ ^^.^^^ 

One deduces easily the expression of the AFM mass 



It is possible to go further in the simplification using the expression of the K function and 
defining the new mean radius 

2/(A+2) 

\/xo 
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The transcendental equation (15- 191 ) and the value of the AFM mass (15-201) are given now by 
simpler expressions: 



CTe, = r2v'(ro)^l+(^) , (5-22) 

MAPM = ^^l+(^)%nro). (5-23) 

Despite the fact that two auxiliary fields were introduced, only one transcendental equation 
is necessary to solve the problem. This great simplification is due to the use of a power-law 
expression for the basic function P{r). 

It is also clear from (15-221 ) and (15-231 ) that, very much in the same way than the nonrel- 
ativistic case, the expression of the mean radius and of the mass depend on the A variable 
only through the value of the principal quantum number Q{X;n,l), which can thus be mod- 
ified at will in order to improve the results. So, in the following, this quantum number will 
be simply denoted Q. 

As in the nonrelativistic case, the AFM also yields approximations for the eigenstates. 
The approximant states, which are eigenstates of p^/ /Xo + Vosgn(A)r''', are characterized by 



a size parameter depending on the product /loVo (see (IE-31 ). (IF-31 ) and (IG-31 )). This quantity 
is given by 

\M4 



A lot of observables can then be computed in terms of ro only. 

Another procedure to solve the problem is to start from formula (15-81 ) and replace the 
exact value by a corresponding approximation given by an AFM solution obtained 
from the Hamiltonian/j(/i,Vo) = p^//i + VoP(r) with/'(r) = sgn(A)r''-. After this first step, 
the optimal value Vq depends not only on V{r) but also on /i. One can then show that the 
subsequent extremization on /i gives exactly the results (15-221 ) and (15-231 ) presented above. 

This procedure in two steps is nevertheless interesting to obtain informations about the 
bounds on exact energies. The replacement of the square root operator by a nonrelativistic 
counterpart without modifying the potential, H ^ H{ijl) (see (15-11 ) and (15-71 )). yields upper 
bounds. If the exact potential V{r) is then replaced by an AFM approximant V{r) >'V{r), 
the upper bounds will be further increased. On the contrary, if the new potential is such that 
V{r) <V{r), the upper bounds will be decreased. As a consequence, the energies given by 
(15-221 ) are upper bounds on the exact energies if V{r) > V(r), or equivalently if the function 
g{x) = V{P^'^{x)) is concave. When V{r) < V{r) or g{x) is convex, the antivariational 
character of the approximation cannot be guaranteed. 

The previous equations take a particularly simple form for massless particles. Indeed, 
for m = 0, they reduce to 

<yQ = rlv\ro), (5-25) 

Mafm = — + V(ro). (5-26) 
''0 
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Let us define the function D{x) as tiie inverse function of x^V'(x) 

D{x^V'{x))=x or {D{x)fv'{D{x))=x 
and tlie related function F (x) defined by 



(5-27) 



F(x) 



D{x) 



+ V{D{x)). 



(5-28) 



Once the potential V is given, these functions are universal in the sense that they do not 
depend on the mass m (except if the potential is mass-dependent) and can be computed once 
and for all. The mass for a system in the ultrarelativistic limit is thus written in a very simple 
form, namely 

MAFM=^'(CTe). (5-29) 
By rewriting (I5-22I )- (I5-23I ) and denoting the kinetic energy by r(|p|) with T{x) = 



o^/x^ + nfi, the AFM formula for the energies can be written into the form 

MAFM = r(po) + V(ro), 
Q 

Po = — , 
ro 

PqT'{pq) = roV'iro). 



(5-30) 
(5-31) 
(5-32) 



As in the nonrelativistic case, (15-311) defines the mean impulsion po from the mean radius tq. 
This radius is defined by (I5-25I ). and one can recognize again in (15-321 ) the general form of 
the virial theorem.ES Finally, (15-301 ) gives the energy as a sum of the kinetic energy evaluated 
at the mean impulsion po and the potential energy evaluated at the mean radius tq. These 
equations simplify greatly for m = and reduce to (I4-15I )- (I4-17I ). with Mafm = 2/m + ^afm. 
for m — )• oo and a = 2 (mis then the common mass of the two identical particles and not the 
reduced mass). 

The parameter jM) can be considered as an effective mass for the relativistic parti- 
ciem^JQ}^,32> Using (ISTTl) . (IFITl) and (IS^TT) . one finds 



Mo = -\Jpl + nfi 



(5-33) 



which is a kind of "total energy of a free particle". 

After some algebra, one can check that the virial theorem applied to Hamiltonian 
H{}io, Vo), whose solutions are M{pLQ, Vq) given by (15-151 ). implies that 



Aio,Vo 



Pi 



Ato,Vo 



Ato,Vo 



Mo,Vo 



1, 



(5-34) 



where |/io, Vq) is an eigenstate of H{iJ.o,Vo)- These equations, as well as the boundary char- 
acter of the solution, are only applicable when the exact form Q{X;n,l) is used. In practice, 
this occurs when A = —1 or 2, or A = 1 for S-states only. Nevertheless, a better accuracy 
can be obtained by an appropriate choice of Q as we will see below. 
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5.3. Power-law potentials 
5.3.1. General case 

The power-law potential was considered as a prototype for testing the AFM in the 
framework of the Schrodinger equation. It was shown that an analytical expression for 
the eigenenergy is obtained whatever the power A under consideration. We apply in this 
section the AFM to find approximate analytical energy formulae for the spinless Salpeter 
Hamiltonian (15-11) with the potential V{r) = sgn{X)ar^ where a is positive. We consider 
that A > —2, as needed in the nonrelativistic case to have bound states. In this cases, (15-221 ) 
and (15-231 ) become 



(jQ^ = a\X\r^+'jQ^ + m^rl (5-35) 



M(A) = cyW^ + m2 + sgn(A)ar^. (5-36) 



With the change of variable 



(15-351 ) can be recast into 



r2 



a\X\ 



2 

A+2 4 

G^^^o, (5-37) 



—x^+^-A^xo-m^ = with A^ = (-^j Qtti , (5-38) 

a form which will be very convenient to use. With this convention, the mass (15-361 ) is also 
written in a simple form 

MiX) = ^^!^^^pt±M^. (5-39) 

In order to obtain analytical energy formulae, one should be able to analytically solve 
(15-381 ). Let us now examine the cases for which this is possible. First, we set X + 2 = p/q 
and Xq = xH'^, where p and q are co-prime integers. Then, (15-381 ) becomes 

^XP-A^X^^-m^ = 0. (5-40) 

A polynomial possesses analytical roots if its order is less or equal to 4; therefore all the 
solvable potentials should verify the conditions 0<p<4, 1<^<4. An exhaustive research 
of all the solvable potentials leads to the following values for the power A : 

7 5 3 4 5 2 1 

Among these allowed values, three are of particular interest: A = — 1 and 2 because the 
Coulomb problem and the harmonic oscillator play a central role in theoretical physics, and 
A = 1 since a linearly rising potential is generally considered to be a relevant approximation 
of the confining potential in QCD.E3 These three cases are explicitly solved in the following 
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sections. Thus, in contrast to the nonrelativistic case where an analytical expression exists 
for an arbitrary value of A, the spinless Salpeter equation leads to a restricted list of favorable 
values (I5-41I) . Fortunately, the most interesting potentials belong to this list. 

Let us briefly discuss the existence of physical solutions, as function of the parameters. 
Starting from (15-381 ) and ( I5-39I ). it is possible to prove the following properties: 

• For A > 0, there exists always a physical solution with M(A) > 0. 

• For — 1 < A < there exists a physical solution with M(A) > as long as the parameter 
a is less than a critical value ^^.(A) given by 



Let us remark that ac(A) depends on m, except ac{—X) = oQ. 
• For — 2 < A < — 1 either there is no solution for (I5-38I ). or there exists two solutions 
but, in the latter case, these solutions aie not compatible with the nonrelativistic ex- 
pressions when m tends toward infinity. In both cases, the corresponding solutions are 
not physical so that these types of potential must be discarded from our study. 

In the following, we will denote Qx = Q{X;n,l) the best possible form of the principal 

quantum number for the power-law potential V{r) = sgn(A)<3r'*' in a Schrodinger equation. 

We have then Q2 = Qho^ Q \ = Qc and Qi given by (14-81 ). For other values of A, a formula 

of type (13-361 ) can be chosen for instance. 

5.3.2. Harmonic potential 

If A = 2, (15-381) becomes a quartic equation of the form 

—Xq -A2X0 -n? = with A2 = y/aQ. (5-43) 

Defining 



'0=(^^j X and , (5-44) 

(15-431 ) can be rewritten as 4X^ — SX — 3Y2 = 0, which is precisely of the form (IB-61) . Follow- 
ing (IB-71 ). the solution is given by X{Y2) = G-{Y2), and, after a rearrangement of (I5-39I ). the 
mass spectrum reads (under one of the equivalent forms) 

M(A = 2) = a.,/I^±i£=i£L (5-45) 

2am / 9 , , 1 \ 

With Q = Q2 these formulae yield an upper bound on the energy, since this setting gives the 
exact solution for the corresponding nonrelativistic Hamiltonian. 
In the nonrelativistic limit (7 — )• 00), the last equations reduce to 



2a a 

M-om^\ e, (5-48) 

m 
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as expected for a nonrelativistic harmonic oscillator. For large value of m, the choice Q = Q2 
is clearly optimal. For small values of m, another choice could give better results. 

5.3.3. Linear potential 

The resolution of the case A = 1 is rather similar to the one of the harmonic oscillator. 
Defining 



(15-381 ) simply becomes a cubic equation of the form (IB lb . that is — 3X — 2Y\ = 0. Notice 
that A 1 = {aQ/lf/^ so that 

r,=|^. (5-50) 
2 aQ 

Following (IB-21) . the solution is given by X(Y]) = F_(Fi). After a rearrangement of ( I5-39I) . 
the mass spectrum reads (under one of the equivalent forms) 

M(A = 1) = am.ll (5-51) 

= omJ I 3 + 1 (5-52) 

'3 + f2(F^)). (5-53) 



V2FiF_(Fi) 

In the nonrelativistic limit (F — oo), the last equations reduce to 



3 /aa222\ 1/3 
M^am + -[ —] , (5-54) 



m 

as expected for a nonrelativistic Hamiltonian with a linear potential. For large value of m, 
the choice 2 = 2i is a good one. For small values of m, another choice could give better 
results. 

5.3.4. Coulomb potential 

Equation (15-381) is considerably simplified when A = — 1 ; the value of xq can directly 
be extracted in this case and reads 

4m2 . . . ^ a y 



The energy spectrum ( 15-391 ) is finally given by 



MiX = -l) = am^l-^. (5-56) 

It is obvious from this last equation that bound states exist only when a < aQ, in agreement 
with (15-421 ). The most stringent upper bound for a is actually found for the ground state, 
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that is a < a 2|„^;^o- With Q = Q \, (15-561 ) yields an upper bound on the energy since this 
setting gives the exact solution for the corresponding nonrelativistic Hamiltonian. This is 
confirmed by the results of Refs. 17711.1781) in which formula (15-561 ) was already given. 

As the coupling constant a is dimensionless in this case, the only mass scale of the 
Hamiltonian is m. So the nonrelativistic limit cannot be obtained by setting m — > oo, as usual. 
It is well known that the mean speed of the particle is independent of m and proportional 
to a in the Coulomb case. So, the nonrelativistic limit is achieved for a — 0. When a <^\, 
(15-561 ) reduces to 

M^crm-^^, (5-57) 

which corresponds to the rest energy plus the Coulomb binding energy, as expected. For 
large values of m, the choice Q = Q \ is clearly optimal. For small values of m, another 
choice could give better results. 

Several works have been devoted to the spinless Salpeter equation with a Coulomb po- 
tential,'^''^''^® and it is interesting to compare our results with previously found ones. 
First of all, it has been shown in Ref. 1791) that the energy spectrum of the semirelativistic 
Coulomb problem is unbounded from below if 

2a 

a> = — . (5-58) 

71 



Moreover, a lower bound on the ground-state energy Mg is^ 

(5-59) 




A last result of interest is the analytical determination of the ground-state energy per- 
formed in Ref. 81). It is found that 

Mg\a=a,<CJmxOAS42564... (5-60) 

Our formula (15-561) leads in this case to 



M\a=a, = omJ\-——^ . (5-61) 

V Q\n^l=0 

This last expression is equal to 0.77 om if Q = but agrees with (15-601 ) if Q\„^i^q ~ 0.73 
is taken. 

5.3.5. Ultrarelativistic limit 

In the ultrarelativistic limit, that is to say m = 0, the various elimination equations for 
the auxiliary field and the energy formulae become simpler. Equation (15-381 ) is then written 

2 

(^4+'-Aa) =0, (5-62) 
and the final energy spectrum is given by (using the non trivial solution of (15-621 )) 

Ml\a,a,N) = ^ (a|A|)x^ (oQ)^ . (5-63) 
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Let us stress that, in this special case, an analytical expression is obtained whatever the 
power A of the potential, in contrast to the particular values (15-411 ) resulting from the general 
case. One can check that (15-451) and (15-511) reduce to (15-631 ) when m — 0. 

The quantity M™ is physically a mass. It appears that < when — 1 < A < and 
M'^y = 0. No bound state of massless constituent particles can be found in these cases. 
Although (15-631 ) is positive for —2 < A < — 1, this range of values for A has been proved to 
be unphysical. 

For the particular case A = 1, this last formula becomes M"^ = l^oaQ, and one has 
consequently (MJ"^)^ oc Q. Such a linear behavior between the squared masses and the quan- 
tum numbers of the different states (with Q = Q\ for instance) is a well-known property of 
the spinless Salpeter Hamiltonian with a linear potential and massless constituents. Such 
a Hamiltonian is one of the simplest ways to describe a light meson in a potential ap- 
proach of QCD, and it is experimentally checked that the square masses of the light mesons 
mainly grow linearly with their angular momentum (Regge trajectories). See for example 
Refs.l47l).l82l) for a discussion of that point. 

It is worth looking at the ultrarelativistic harmonic oscillator. The energy spectrum of 
the Hamiltonian 

}i^,o = l\/~p^^a? (5-64) 
can be analytically computed for / = and readd^ 

Mho = -(4a)i/3a„, (5-65) 

where a„ < are the zeros of the regular Airy function (see Appendix [G]). Using (I4-8I ). 
(15-651 ) reads 

Mho = 3(V^e(l;«,0))^/\ (5-66) 
This last expression can be compared to the AFM result (15-631) for A = a = 2, i.e. 

M^' = 3{V^Qf\ (5-67) 

With (15-641 ) as a starting point, a natural choice for 2 is Q2, but it is clear that both Mho 
and are identical provided that Q = Q{l;n,0), the principal quantum number giving the 
exact energy formula in the case of a nonrelativistic kinetic energy with linear potential. The 
explanation of this fact is that the Fourier transform of Hamiltonian (15-641) is a nonrelativistic 
Hamiltonian with linear potential. This is a supplementary indication that Q should be 
different in the semirelativistic and nonrelativistic cases for the same potential. 

5.3.6. Improved formulae for ultrarelativistic power-law potentials 

Exactly as it was the case in the nonrelativistic framework, it is possible to improve 
drastically the energy formulae of the spinless Salpeter simply by changing the value of the 
principal quantum number Q. To study this problem, we employ the scaling law properties 
and consider the dimensionless Hamiltonian valid for massless particles 

H = 2y^ + r^ (5-68) 

with A > 0. Very accurate numerical values for the eigensolutions of this Hamiltonian can 
be obtained with the Lagrange mesh method.ES 
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The approximate energy spectrum is given by (15-631) with a = I and a = 2. With the 
choice Q = Q2 = 2n + 1 + 3/2, upper bounds ai^e obtained. As mentioned before, another 
choice for the n- and /-dependences of Q can greatly improve the results. By using the form 

Q = b{X)n + d{X)l + c{X), (5-69) 

we find smooth variations for coefficients b, c and d for A €]0,2] (/ < 3 and n < 3). In 
particular, d{X) « 1 with relative variations less than 2%. Coefficients b{X) and c(A) can 
be fitted with various functions and similar agreement. Finally, we choose 

3.00 A +3.67 2.69A + 8.69 

A + 3.40 ' , ^W = l. (5-70) 

Agreement with exact results is very good but the variational character of the approximation 
is no longer guaranteed. With the choice (I5-70I) . the maximal relative error for / < 3 and 
n <3 and for A G [0.1,2] is located between 0.3 and 1.1%. With the choice Q = Q2, the 
corresponding error is located between 4.5 and 12.7%. Just to give a quantitative idea of the 
quality of the approximation, results for A = 1 are presented in Table IXII 



Table XI. Eigenvalues eln. I) of the Hamiltonian J5-68t with X = 1 , for some sets (n,/). First line: value from 
numerical integration (exact values); second line: approximate result l l5-63b with Q defined by i5-69t and 
l l5-70b : third line: approximate result J5 63t with Q = 2n + 1 + 3/2 (upper bounds). 



I 


£(0,Z) 


£(1,0 


£(2,/) 


£(3,0 





3.1577 


4.7109 


5.8913 


6.8742 




3.1338 


4.6849 


5.8374 


6.7973 




3.4641 


5.2915 


6.6333 


7.7460 


1 


4.2248 


5.4575 


6.4837 


7.3767 




4.2215 


5.4725 


6.4866 


7.3623 




4.4721 


6.0000 


7.2111 


8.2462 


2 


5.0789 


6.1304 


7.0470 


7.8671 




5.0814 


6.1602 


7.0764 


7.8869 




5.2915 


6.6333 


7.7460 


8.7178 


3 


5.8108 


6.7425 


7.5775 


8.3387 




5.8156 


6.7785 


7.6207 


8.3787 




6.0000 


7.2111 


8.2462 


9.1652 



It was shown that Q\i^q must be given by (14-81) for a quadratic potential. Formu- 
lae (15^ giveb{2) = 1.79 close to 7r/vt« 1.81 and c(2) = 1.37 close to 71^/3/4 « 1.36. 
This is also in agreement with formulae (71) in Ref. |38]) which predict, in the case of a 
nonrelativistic Hamiltonian with a linear potential, b{l) = 1.79 and c(l) = 1.38. 

5.3.7. Improved formulae for relativistic Coulomb potential 

In this section, we study the relativistic Coulomb potential and use the scaling properties 
to set m = 1. With dimensionless variables, the semirelativistic Coulomb Hamiltonian is 
written in this case 

h = 2^/p^ + l-- (5-71) 
r 
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with <a <ac. Again, very accurate numerical values for the eigensolutions of this Hamil- 
tonian can be obtained with the Lagrange mesh method. 



Table XII. Eigenvalues e(n, /) of the Hamiltonian l l5-71t with a = 1, for some sets {n, I). First line: value from 



numerical integration (exact values); second line: approximate result ( l5-56b with Q defined by l l5-72b and 
( l5-73b : third line: approximate result l l5-56t with Q — n + l + I (upper bounds). 



/ 


£(0,Z) 


£(1,0 


£(2,0 


£(3,0 





1.65817 


1.92184 


1.96739 


1.98231 




1.65982 


1.92356 


1.96680 


1.98151 




1.73205 


1.93649 


1.97203 


1.98431 


1 


1,93515 


1.97122 


1.98389 


1.98973 




1.93476 


1.97012 


1.98291 


1.98895 




1.93649 


1.97203 


1.98431 


1.98997 


2 


1.97187 


1.98416 


1.98987 


1.99297 




1.97296 


1.98416 


1.98961 


1.99266 




1.97203 


1.98431 


1.98997 


1.99304 


3 


1.98428 


1.98993 


1.99301 


1.99487 




1.98528 


1.99021 


1.99302 


1.99477 




1.98431 


1.98997 


1.99304 


1.99489 



The approximate energy spectrum is given by (15-561) with m = 1 and a = 2. With the 
choice Q = Q i =n + l + 1, upper bounds are obtained and the results are exact in the limit 
a — )• 0. As mentioned before, another choice for the n- and /-dependences of Q can greatly 
improve the results. By using the form 



Q = b{a) n + d{a) I + c{a) 



(5-72) 



and imposing b{0) = c(0) = d{Q) = 1, we find smooth variations for coefficients b, c and d 
for a G]0,ac] (/ < 3 and n < 3). These coefficients can be fitted with various functions and 
similar agreement. Finally, we choose 



bia) 



1.03a- 1.48 
a -1.48 



c[a 



1.07a -1.64 
a -1.64 ' 



d(a) 



0.96a- 1.56 
a -1.56 



(5-73) 



Agreement with exact results is very good for a < 1.2 but the variational character of the 
approximation is no longer guaranteed. With the choice (I5-73I ). the maximal relative error 
for / < 3 and « < 3 and for a G [0.2, 1.2] is located between 0.005 and 0.3%. With the 
choice Q = Q i, the corresponding error is located between 0.004 and 17.3%. To obtain 
a good accuracy in the domain a ?s a^. = 4/?! ?a 1.273, special method must be used as the 
one presented in Ref. ,81J . To recover the value obtained for the ground state in this paper, 
it is necessary to have c{ac) = 0.73. Our formula gives c(ac) = 0.77. Results for a = 1 are 
presented in Table IXIll 

One can see that Q = Q \ = ?i + / + lisa better choice for large values of n or /. This 
can be understood as a kind of nonrelativistic behavior since the limits a — )• and °° 
have similar effects. 



63 



5.4. Square root potential 

In Sect. 14.31 we introduced the square root potential V{r) = Va^r^ + and commented 
its important role in the theoretical description of hybrid mesons in a nonrelativistic frame- 
work. Nevertheless, such a type of potential allied with a relativistic kinetic energy operator 
is important to describe hybrid mesons with low mass quarks as well. The general formula 
with a non vanishing quark mass looks rather difficult and we prefer here to present the two 
extreme limits and show the close connection with the previous section through a Fourier 
transform. 

5.4.1. Nonrelativistic kinematics 

The nonrelativistic limit follows from the Hamiltonian H = ^ + Va^r^ + and as 



with 



been studied in Sect. 14.31 It is easy to show that the Fourier transform of this Hamiltonian, 
denoted as H^^, reads 

= (jy'p^+M^- + Kr^, (5-74) 
4a \ b ma 

— , M = -, K = a—. (5-75) 
m^y a 8 

This last Hamiltonian is nothing else than a spinless Salpeter one (15-11 ) with a harmonic 
potential. Provided that the proper substitution rules are taken into account, the energy 
spectrum in this case is the same than the one computed in Sect. 15.3.21 We have thus 

E- = b,[I^±^2=m= (5-76) 
V Y2 ^8G_(F2) + 3F2 

' ' ' (5-77) 



^G-iY2) + jA^], (5-78) 



with 



Y2 \G4Y2) 2Gl{Y2 
3^(^-(^^) + d^ 

b^ / 32/71 



a^Q^ 



(5-79) 



This potential was studied in more detail in Sect. 14.31 and we recovered here the results 
presented there (see (14-261 )). This remarkable result proves that AFM is an approximation 
consistent with the Fourier transform. 

5.4.2. Ultrarelativistic limit 

In the ultrarelativistic limit (m = 0), the corresponding Salpeter Hamiltonian is H = 
a \pp^ + \J d?-r^ ^h^. To get the eigenenergies, it is not necessary to start a new calculation. 
In fact, the Fourier transform of H is simply 



= OyJp'^+M'^ + ar, (5-80) 



with M = b/o. But this last Hamiltonian appears to be a spinless Salpeter one (15-11 ) with 
a linear potential. Provided that the proper substitution rules are taken into account, the 
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energy spectrum in this case is the same than the one computed in Sect. 15.3.31 We have thus 

M- = ^#^±^^i^ (5-81) 



^1 




' * (5-82) 

(5-83) 

with 

^3/2 a2 

20 0^/ 

When — 0, these last equations logically reduce to (15-631 ) with A = 1 . 

5.5. Funnel potential 

As we stressed in Sect. 14.6.71 the funnel potential Vif) = ar — b/r, with a and b both 
positive, is of crucial importance in hadronic physics. It was shown that an analytical ex- 
pression is available in a nonrelativistic framework. However the glueball resonances or the 
spectroscopy of very light mesons need a semirelativistic description and a spinless Salpeter 
equation governed with a funnel potential is the most economic way to consider this situa- 
tion. 

5.5.1. Massless particle 

The ultrarelativistic limit {m = 0), is particularly important for the description of a 
glueball composed of two massless gluons or a meson composed of two massless quarks. 
Then, the extremization condition (15-221 ) becomes ar^ + b = oQ whose solution is 



oQ-b 

?-o = y^^. (5-85) 
Inserting this value in the expression of the mass (15-231 ) leads to the mass formula 



Mf = 1 yJa{oQ - b). (5-86) 

It is astonishing that such a complicated potential used in conjunction with a semirelativistic 
kinetic energy admits analytical approximate eigenenergies of such simple form. For Z? = 0, 
(15-861 ) reduces to l^aoQ, that is the expected expression from (15-631) in the case A = 1. 
For £2 = 0, — )• as shown in Sect. 15.3.51 If the prescription 2 = 22 is chosen, then 
(15-861 ) is an upper bound on the exact result. As already mentioned, other choice for Q could 
improve the accuracy of the formula but without guarantee about the variational character 
of the result. This point is studied in the next section. 

5.5.2. Improved formulae for the ultrarelativistic limit 

Written with dimensionless variables, the ultrarelativistic Hamiltonian for two massless 
particles interacting via the funnel potential is given by 

2 ^ 



+ (5-87) 
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with j3 > 0. The approximate energy spectrum is given by (15-861) with a = 2, a = I and 
b replaced by j8 (to avoid confusion with the coefficient of n). This kind of Hamiltonian 
is often used in hadronic physics with typical values for j8 ^ 0.4. Very accurate numerical 
values for the eigensolutions of this Hamiltonian can be obtained with the Lagrange mesh 
method.E3 

With the choice Q = Q2 = 2n + 1 + 3 /2, upper bounds are obtained. As shown in 
previous cases, another choice for the n- and /-dependences of Q can greatly improve the 
results. By using the form 

Q = b{p)n + d{l5)l + c{l5), (5-88) 

we find smooth variations for coefficients b, c and for jS G [0, 1] (/ < 3 and n <3). These 
coefficients can be fitted with various functions and similar agreement. Finally, we choose 

1.88/3 -5.34 1.99/3 -4.40 0.76/3 -2.46 

^^^^= /3-3.51 ' '^^^= /3-3.49 ' ^^^^ = /3-2.54 " ^^'^^^ 

Agreement with exact results is very good but the variational character of the approximation 
is no longer guaranteed. With the choice (I5-89I ). the maximal relative error for / < 3 and 
n <3 and for j3 G [0, 1] is located between 0.6 and 4.9%. With the choice Q = Q2, the 
corresponding error is located between 12.7 and 42.2%. Results for j8 = 0.4 are presented 
in Table IXim 

For j8 = 0, one obtains b = 1.52, c = 1.26, and d = 1.09. From (15-701) with A = 1, one 
obtains b = 1.52, c = 1.23, and d = I. These values are close to each other as expected. 

Table XIII. Eigenvalues e{n,l) of the Hamiltonian ( l5-87t with j8 = 0.4, for some sets {n,l). First line: value 



from numerical integration (exact values); second line: approximate result ( l5-86t with Q defined by fS-l 
and ( l5-89l l: third line: approximate result l l5-86b with Q = 2n + 1 + 3/2 (upper bounds). 



I 


e(0,/) 


£(1,/) 


£(2,/) 


£(3,0 





2.7821 


4.3709 


5.5874 


6.5938 




2.7804 


4.4196 


5.5977 


6.5678 




3.2249 


5.1381 


6.5115 


7.6420 


1 


3.9944 


5.2365 


6.2744 


7.1772 




3.9737 


5.2529 


6.2765 


7.1552 




4.2895 


5.8652 


7.0993 


8.1486 


2 


4.8993 


5.9549 


6.8772 


7.7028 




4.8837 


5.9710 


6.8887 


7.6978 




5.1381 


6.5115 


7.6420 


8.6255 


3 


5.6588 


6.5927 


7.4311 


8.1957 




5.6489 


6.6115 


7.4508 


8.2046 




5.8652 


7.0993 


8.1486 


9.0774 



5.5.3. General case 

In order to treat the general case 



H = aVp^ + m^ + ar--, (5-90) 

r 
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we use the dimensionless Hamiltonian h = H / \/o a given by 



/2 = \/F+F + '---, (5-91) 



with % = m^Joja and j3 = b/o. The transcendental equation (l5-22b writes in this case 



Q = {ri + ^)\l^ + {Xro/Qf. (5-92) 
Defining a new variable x and a new parameter a by 

2 2 

^=|2^2 . a = |^^, (5-93) 

the previous equation can be recast as a third order polynomial equation: 

+ {\+2a)x^ + {a^ + 2a)x + a^{\-Q^/^'^) = 0. (5-94) 



In consequence, it is analytically solvable and can be put under the reduced form (IB lb . 
Reporting the obtained x value in the expression (15-231) for the mass leads to the final ex- 
pression. 

Explicitly, the procedure to get the analytical AFM expression for the spinless Salpeter 
equation based on the Hamiltonian ( I5-91I ) needs the following steps. Calculate the a value 
from (15-931 ). then introduce the parameter 

The X value is obtained from the expression 

x='^[\a-\\F^{Y)-{l+2a)], (5-96) 



while the AFM mass is given by 

1 



M. 



AFM 



X^\/l+x+{x-a)Q 



(5-97) 



Once again, one sees that the AFM is very powerful to get analytical expressions, even for 
quite sophisticated problems. We have checked that (15-861 ) [with a = 1, a = 1 and b = {5] 
is recovered in the limit ;t — )• 0. It can also be checked that (15-971 ) and (14-1091) [with a = 1, 
m = X, a = I and b = {5] tend towards the same limit when X ^ °°- 

5.5.4. Low mass expression 

By a "small mass", we mean the condition m <C ^/a. For m = 0, we saw already that 
ro = ^ (oQ — b^ja. The principle of the method is based on a limited expansion of the 
equations at first order in n?. Thus, we set 

rl=^-^{\^z). (5-98) 
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From the transcendental equation, it is easy to show that 

Inserting these values in the expression (15-231) for the mass, we obtain the final result 

^ oQ —b ( on? 



Here again, it is amazing that the approximate eigenvalues of a so sophisticated Hamiltonian 
take such a simple form. 

5.6. Unequal masses 

Some particular problems require to deal with a system of two particles with unequal 
masses. In this case, a general spinless Salpeter Hamiltonian is given in the rest frame by 



H=^p^+ml + ^p^ + ml + V{r). (5-101) 

A general result can be obtained for Hamiltonians with two different masses. Let us 
consider the following two-body Hamiltonians H = Ti + T2 + V , Hi = 2T\ + V , and H2 = 
2T2 + V whose ground-state energies are respectively M = {<^\H\<^) , My = {^i\H\\^i), and 
M2 = (^21^21^2)- Since H = {Hi+H2)/2, we can write 

= i ((0i//ii0> + ■ (5-102) 

The Ritz theorem implies that 

M >^{Mi+M2). (5-103) 

For particular cases, this approximation can be quite good. In Kef's?)), it is shown that M ^ 
{Ml +M2)/2 for a relativistic Hamiltonian of kind (15-1011 ) with V{r) = ar and ^ ^/a. 

The square roots appearing in the kinetic terms can still be treated by resorting to the 
AFM. But this time, two auxiliary fields, pLi and IJL2, have to be introduced. One is led to the 
Hamiltonian 

/I1+/I2 , ni\ ml ^ 



H{ixi ,1x2) = ^^^r^ + ^ + Tr^ + - 



2 2iii 2^2 2M {1x1,1x2) 
+sgn(A)vr^ +y(7(v)) - sgn(A)v/(v)^ , (5-104) 

with 

M(Mi,M2) = -^ (5-105) 
Ml +M2 

playing the role of a reduced mass. It has been shown in Ref . [85]) that the AFM leads to the 



system (I5-30I )- (I5-32I ). but with T{x) = yjx^ + m\ + _^ ^2 

It is then easy to obtain numerical upper bounds, but one can wonder if it is possible to 
obtain analytical solutions? Due to the presence of the square roots, it is only conceivable 
with very special conditions: one mass is vanishing {mi =0 and m2 = m) and the potential 
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is well chosen. Analytical solutions can be obtained for the funnel potential, often used in 
hadronic physics.!^ Let us note that an analytical upper bound for a logarithmic potential 
(also relevant in hadronic physics) can also be computed. 

The solution for the general funnel potential is very complicated. As an illustration, 
the case of the Coulomb interaction, —ajr, is only presented here. The solution for the 
associated system (I5-30I) - (I5-32I) gives 



Q ^a{lQ-a) 

ro = -- (5-106) 

m a — Q 



and 



This formula is an upper bound if Q = n + 1 + I, but, as we pointed out several times, it is 
justified to take a more sophisticated expression in order to get a better accuracy. It can be 
seen from (15-1071 ) that limm^o^ = 0, as expected. 

From equations above, a solution M exists only when the following condition is verified: 

^<Q<a. (5-108) 

The right part of this inequality corresponds to the cancellation of the binding (M ^ m 
when Q — a), and the left part corresponds to a collapse (M — when Q — a/2). Indeed, 
the mean radius ro — )• 0° when Q ^ a, and ro — )• when Q ^ a/2. A sufficiently strong 
interaction must exist to bind the system, but unphysical values of the mass could appear 
if it is too strong. With the Coulomb potential, the nonrelativistic limit is not defined by 
wi — )• 00 but by a — > 0. So, (15-1081 ) implies that this limit is irrelevant for the solution (15-1071 ). 

5.7. Semiclassical interpretation 

A semiclassical interpretation of the system (I5-30I )- (I5-32I ) is also possible. Let us as- 
sume a classical circular motion for the two particles, as illustrated in Fig. H] In this case, 
the force Fj acting on particle / is given by 

Fi = m^y{vi)^-^ with ^ = ^ . (5-109) 
at at ri 

Taking into account that both particles are characterized by the same momentum pQ, This 
equation becomes 

^- ^0 ^ (5-110) 



pl + m^' 

The rigid rotation constraint, v\/r\ =V2/r2, implies that 



lpl + m\ + Jpl + ml 

ro = r,+r2 = ri^ , ^ , (5-111) 

Ipl + m] 
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with / / j. If the force acting on / comes from the potential V{r) generated by j, then 
F\ = F2 = V(ro). (15- 1101) and (151111 ) can be recast onto the form (15-321) . and it is obvious 
than (15-301 ) gives the mass of the system. The total orbital angular momentum is r\ po + 
''2 7^0 = fopo- A semiclassical quantification gives thus ropo = L + 1/2, and we obtain a 
system very similar to (I5-30I )- (I5-32I ). Nevertheless, the AFM produces more general results. 



m 




Fig. 4. Classical circular motion of the two relativistic particles. 



§6. A'^-body problems 

Up to now, the AFM was applied essentially to one- and two-body systems with ei- 
ther a nonrelativistic or a relativistic kinetic energy operator. It was shown that, in many 
situations, the method was able to provide an accurate analytical expression for the corre- 
sponding eigenenergies. In this section we address the important point of whether AFM 
is also able to give informations concerning the eigenenergies of a system composed of N 
particles. Except the very special case of particles whose dynamics is governed by quadratic 
interactions, which is treated in detail in Appendix HI no analytical solution is known for 
A'^-body systems. Nevertheless a number of very interesting physical problems need to solve 
a A/^-body equation. As in the 2-body case, the AFM cannot pretend to give a very high 
accuracy in the results, but it exhibits the interesting virtue to clearly show the dependence 
of the results in terms of the physical parameters and the various quantum numbers. 

Indeed, the envelope theory was used in the A/^-body problem to find bounds for the non- 
relativistic binding energies if the two-body potential fulfills certain restrictive conditions.!^ 
In this section, we apply AFM to the A'^-body problem in a more general context 

• we consider the possibility of dealing with a semirelativistic kinetic energy potential; 

• we allow the presence of one-body terms in addition to two-body interactions. 
Moreover, we are not interested only by the existence of bounds which can sometimes 

give rather poor values, but the purpose of the present section is to show that the AFM can 
be successfully applied to find approximate analytical mass formulae for general relativistic 
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A/^-body Hamiltonians of the form 



H = Y^^Jp] + m] + Y^Ui{\ri-R\)+ £ V,-y(|r,- - r,-|), (6-1) 

!=i (=1 i<i=i 



with Y!i=\ Pi = P = 0, where r, and pi are respectively the position and momentum of par- 
ticle / with a mass m,, and where i? is a global variable defined in Appendix U As stated, 
we consider a relativistic kinematics and allow for both pairwise V and one-body U interac- 
tions. Particular systems exhibiting one-body potentials of the previous form will be studied 
below. Most of the material of this section can be found in Ref. 68 ). 

A quantity which appears very often in this section is the number of interacting parti- 
cles. In order to simplify the notations, we use a special symbol to designate this number 
and set 

C. = ^. (6.2) 

6.1. Main mass formulae 
6.1.1. General case 

As stated, the AFM can be straightforwardly generalized to the case of N-hoAy Hamil- 
tonians of the form (16- II ). The basic idea is to introduce auxiliary fields so that this Hamilto- 
nian is formally replaced by a Hamiltonian for which an analytical solution can eventually 
be found. The auxiliary fields are denoted as /x,, v,, and p,y, and are introduced as follows 



//(At.,v,,pH) = I ^ + ^ + 4-] +I[v,-P(r,) + t/,a(vO)-v,-^(/KvO)] 



N 



+ t [PijS{rij) + Vij{Jij{pij))-pijS{Jij{pij))], (6-3) 

i<j=\ 



where (notice the definition of r,- ^ |r,|) 

r,- = |r,- - /?| , nj = |r,- -rj\, (64) 



and where 



h{x)=K7\x), Ki{x) = ^, J,j{x)=LT.\xl L,Xx) = |^, (6-5) 

where the prime denotes the derivative with respect to the argument. It is sufficient to set 
i^k = "^/fe V/: to treat a nonrelativistic kinematics. We stress that, in order for the AFM to 
apply, it is assumed that the mass spectrum of (16-31 ) is analytically computable. In practice, 
this is only possible for the choice P{x) = S{x) = x^ under some conditions (see Appendix|Jl). 

The reason of such a definition is that Hamiltonians (1611 ) and (16-31) are equivalent if the 
auxiliary fields are properly eliminated. It can indeed be checked that 



8^,.H{Hk,Vk,pki)\^^^^^ = ^ = yP^ + fnj, (6-6) 

5y,H{^k, Vk,Pkl)\y_ = y_ = ^ V/ = ^/(r,'), (6-7) 

dp,jH{iXk,Vk,pki)\p^^^^^^ = ^ Pij = Lijinj), (6-8) 
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and finally that H{jlk, Vk,Pki) = H. 

The approximation underlying the AFM is now that, as we have done in previous sec- 
tions, the auxiliary fields will be seen as real variational parameters. Let us note 

H{Hkfl,Vk,o,Pki,o)\(Po) =Mo\(po), (6-9) 

where Hkfl, Vkfl, Pkifl are the real values of the auxiliary fields such that Mq is extremal. 
\(po) is an eigenstate with fixed quantum numbers and a given symmetry. Using the same 
mathematical techniques than in Sect. [2l one can show that 

^lkfi = \J{(po\pj + mj\(po) with P = 0, (6-10) 
{(po\P{ri)\(j\^) = Pinfl) with r,-o = /,(v,-o), (6-11) 
{(po\S{rij)\(j\)) = S{rijfi) with nj^ = Jij{pijfl). (6-12) 

Using the Hellmann-Feynman theorem, it can be shown that (<po|At(^|<Po) = M,^o> 
{(po\P{Ii{Vi))\(po) =P(/i(v,-,o)) and {(po\S{Jij{pij))\(po) = 5(/;y(pi;-o)). Let us rewrite the 
Hamiltonian (16-31 ) under the form 

N N 

H{^k,Vk,Pki) = T{^k) + Y,Oi{Vi,ri)+ £ Vij{pi,nj), (6-13) 

/=i i<j=i 

where T{iJ.k) stands for the kinetic part. One can also check that 

f/,-(v,-o,r,-o) = Uinfl), 0'i{Vifi,ri^o) = U'in^o), (6-14) 
Vij{pUfi,rij,o) = V(r;y,o), t>/.(p,-,-o,r,7.o) = V'(r,-y,o). (6-15) 

This means that the approximate potential Ui (Vij) and the corresponding genuine potential 
Ui (Vij) are tangent at, at least, one point. If the following conditions Ui{Vifi,ri) > U{ri) 
and Vij{pijfi,rij) > V{rij) are fulfilled V/ and Vj and for all values of the radial arguments, 
Mq is an upper bound on an eigenstate of H (16- lb . If the kinematics is non relativistic 
(T = T{Hk = nik)) and if the following conditions f7,(v;,o,r,) < U{ri) and Vij{p,-jfl,rij) < 
V(rjj) are fulfilled V/ and Vj and for all values of the radial argument, Mq is a lower bound 
on an eigenstate of H. This restriction about T comes from the fact that the replacement 
of the genuine relativistic kinetic operator by the form r(/iyt) yields upper bounds on the 
eigenvalues of the genuine Hamiltonian. In the other cases, it is not possible to obtain a 
relevant information about the position of Mq. 

At this stage, approximate numerical solutions of the A^-body problem can be easily 
computed. First, the choice P{x) = S{x) = %- allows a precise determination of the eigen- 
values of Hamiltonian (16-31 ) (see Appendix IJ]). From now on, this choice is adopted for the 
rest of this section. Second, the extremization of the eigenvalues with respect to the real 
auxiliary fields is a classical numerical problem which can be solved with a high accuracy. 

6.1.2. Identical particles 

To obtain analytical closed solutions of the eigenvalue and extremization problems as- 
sociated to Hamiltonian (16-31 ). it is necessary to simplify the system. First of all, we will only 
consider systems with identical particles, that is with m,- = m. In this case, it is reasonable to 
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consider identical interactions between them, namely Uk{x) = U [x) and Vki{x) =V{x). This 
means that Kk{x) = K{x), Ik{x) = I{x), Lj^i{x) = L{x), and Jki{x) = J{x). The global variable 
R (see (II- 11 1) is then the center of mass of the system, even for relativistic kinematics. 

Let us denote Pij the permutation operator exchanging particles / and j. We can write 
PijWo) = i|'Po) if this state |(po) is completely (anti)symmetrical. Then 

P{rjfl) = {(po\P{rj)\(po) = {(po\PijP{ri)Pij\(j\^) = {(po\P{ri)\(j\^) = P(r,-.o). (6-16) 

If P{x) is monotonic, which is always the case in practice, then rj^ = r,,o- Finally, if I{x) is 
invertible, which must be the case to solve the problem, we have V; o = K{rj o) = K{rj o) = 
Vjfi. So all optimal values V; o are the same. Using the same reasoning, we can draw the 
same conclusion for other auxiliary fields. 

Under these conditions, we can set /i,- = /i, v, = V, and p;y = p in the expression of the 
eigenenergies. We are thus led to replace the original Hamiltonian (16-31) by the following 
simpler Hamiltonian, which now depends on only 3 auxiliary fields (remember that P{x) = 
S{x)=x'^) 



H{^,v,p) = - 



21 

m 

At + — 
At 



+ [U{I{v)) - vl{vf] +Cn [V{J{p))-pJ{pf 



N 



i=l i=l i<j=\ 

which, by virtue of 0-161) . has the following mass spectrum {P = 0) 



(6-17) 



M{^,v,p] 




+ N [U{I{v)) - v/(v)2] +Cn [V{J{p))-pJ{pf 



{v+Np)Q 



with 2= £(2«; + /,) + - (A^-1 



(6-18) 
(6-19) 



(=1 



which is characterized by a high degeneracy due to the form of the principal quantum num- 
ber 

The last step needed to get the final mass formula is to find the optimal values Vq, 
and po from the extremization conditions 



a^M(/i,v,p)|^^^ = 0, dyM{n,v,p)\^^^^ 



0, dpM{^,v,p)\p^p^ 



0. 



Writing explicitly (16-201 ) and after some algebra, we arrive at the mass formula 

M{^,Vo,Po)=N^ + NUiIivo))+CNV{J{po)), 
where the auxiliary fields are solutions of 



m 

Mo = — + 



2Q^Vo+Npo 



1/2 



(6-20) 



(6-21) 



(6-22) 
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^(vo) 
J{po) 



2N^Hoivo+Npo) 
{N-iyiM){vo+Npo)^ 



1/4 



(6-23) 
(6-24) 



At this stage, functions / and / are not known. They depend on the specific forms of U and 
V. Optimal values of Vq, and po are not known anymore. Moreover, they also depend 
on the state considered through the variable Q. Nevertheless, formula (16-211) makes clearly 
appear the mean-field nature of the AFM, because Hq can be interpreted as the average 
kinetic energy of one particle, while U (/(vq)) and V{J{po)) can be respectively seen as the 
average potential energy of one particle in the potential U (r) and of a pair in the potential 
V{r). 

A general remark about the AFM can now be done: The analytical results obtained for 
a given potential U (r) approximated by a potential P{r) can be used as a starting point for 
finding analytical results for an other potential W{r), approximated this time by U (r). The 
same procedure can be applied for the potential V{r). This considerably enlarges the domain 
of applicability of this method. 

Notice that the idea of rewriting a A/^-body Hamiltonian with pairwise interactions of 
the form g{rjj) as a A/^-body harmonic oscillator has already been investigated in Ref. l87l ) 
within the framework of the envelope theory. That method shares many similarities with the 
auxiliary field method, as shown in Appendix |Al but the results presented hereafter have, to 
our knowledge, be presented for the first time in Ref. ,68,). 

6.1.3. Simplified formulations 

The four equations (I6-21I) - (I6-24I ) give the AFM approximation for the mass spectrum of 
a A/^-body system where the constituent particles are identical. They can be even simplified 
by defining 

Xo = ^2^{vo+Npo). (6-25) 



We have indeed 



M(Xo) =Nxlm^ + ^Xo + NU 



NXo 



+ CnV 



2Q 



{N-l)Xo 



(6-26) 



where the remaining auxiliary field is a solution of 



X' 



2,/m2 + |x, 



K\ 



NXo 



+ NL 



2Q 



{N-l)Xo 



(6-27) 



and where K{x) = U'{x)/{2x) and L{x) = V'{x)/(2x) (see (16-51 )). This condition ensures 
that M{Xo) is extremal. In so doing, we have simplified a lot the original formulation. 
The set (I6-22I )- (I6-24I ) is a complicated system of non linear coupled equations, while the 
equivalent formulation (16-271 ) is just a transcendental equation which is a priori easier to be 
solved. As we shall see in the next sections, this task can be achieved in several cases of 
interest. It is worth mentioning that the AFM results (16-261) and (16-271 ) are valid not only for 
the ground state but also for excited states. 
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For particles with null mass (ultrarelativistic limit), one obtains even simpler formulae 
by simply setting m = in formulae (16-261 ) and (I6-27I ). In the nonrelativistic limit, the 
auxiliary field /i tends towards m. In this case also the various formulae look simpler but 
cautions must be taken in the limit. Explicitly, (16-261) reduces to 



with tUt = Nm and where Xq is the solution of 



= 2m 



NXo \y{N-\)Xo 



(6-29) 



Let us introduce the distance ro = \/NQ/Xq and the momentum pQ = Q/r^, and let us 
define T by T{x) = ^/x^ + or by its nonrelativistic counterpart ni + It is a simple 
algebra exercise to show that formulae (I6-26I )- (I6-27I ) can be written aJ^ 

M,, = NT{p(,)+Nu(^-^)+Cmv{^-^^, (6-30) 
Po = -, (6-31) 

'"0 

With this formulation, an AFM eigenvalue given by (16-301 ) is simply the kinetic operator 
evaluated at the mean momentum po plus the potential energy computed at some mean ra- 
dius depending on ro. This generalizes the relations found previously for two-body systems. 
As one could expect, the kinetic energy and the one-body potential energy are proportional 
to the number of particles and the two-body potential energy is proportional to the number 
of pairs. Formula (16-301) looks like a semiclassical approximation but this is absolutely not 
the case. The AFM yields an approximate A'^-body wavefunction,'5SJ'Il2)Jl2} and the relation 
(16-311 ) between po and ro is a full quantum link, function of the quantum numbers of the 
system. Lastly, the value of ro (and thus of po) is the solution of a transcendental equation 
(16-321 ) which is the translation into the AFM variables of the generalized virial theorem^ 
which comes from very general properties of quantum mechanics. These considerations 
prove that the AFM really relies on very sound physical basis. Once the system (16-301) - 
(16-321 ) is written, it can appear finally quite natural to obtain such a result. The problem is to 
find a relevant link between the mean values ro and po. This is solved by the AFM. 

The eigenstates of H{iJ.o,Vo,Po) are built with harmonic oscillator states. A A/^-body 
state with a mass M{Xo) is written 

N-l 

¥=Y[Vnjiji^j,Xj), (6-33) 

where \if„.i.{Xj,Xj) is a three-dimensional harmonic oscillator wavefunction (see Ap- 
pendix 10), depending on the Jacobi coordinate Xj (see Appendix IT]) and decreasing asymp- 
totically like e^^i^il"^ (the magnetic quantum number is omitted). The quantum numbers 
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are such that LyL/ + Ij) = B = Q — ^{N — \) and the scale parameters are given by 




The quantum number B is generally called the band number. The state (16-331 ) has neither a 
defined total angular momentum nor a good symmetry, but its is characterized by a parity 
(— 1)^. By combining states (16-331 ) with the same value of B (or Q), it is generally possible 
to build a physical state with good quantum numbers and good symmetry properties, but the 
task can be technically very complicated.^' ' Some observables can be easily computed: 

--pI (6-35) 




Since these results only depend on the quantum numbers via Q, and since a physical state 
must be a combination of eigenstates with the same value of Q, formulae (16-351) and (16-361 ) 
are also valid for a physical state. This shows that ro can be considered as a mean radius for 
the system and po as a mean momentum per particle. Indeed, (16-351) and (16-361) imply that 



for arbitrary / 7^ j since the mean values are taken with completely symmetrized states. 
These results can also be obtained using the more general relations (I6-10I) - (I6-12I) relevant 
for P{x) and S{x) different from x^. Identifying (16-211) and (16-301 ). it appears that 



^o = ^m^ + p^^. (6-38) 

In the above expression, Q is the principal quantum number defined in ( 16-191 ). However, 
it has been shown many times for the two-body problem that a much better approximation 
of the exact energies can be obtained with a slight modification of the principal quantum 
number. A particularly simple form which seems to work quite well for A/^-body systems is 
given by 

e= + i8Z,-) + 7(A^- 1). (6-39) 

For instance, such a formula is tested numerically for baryons in Sect. 16.4. 3[ it improves 
substantially the results as compared to the expression (16-191 ). It is worth noting that there 
is no systematic procedure to determine the values of parameters a, p and 7, which depend 
on both the interaction and the kinematics. 

6.2. Critical constants 

Some interactions, as the exponential or the Yukawa potentials, admit only a finite 
number of bound states (see Sects. l4A2] and l4~5.2l ). Let us assume that such an interaction 
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can be written as W{x) = —Kw{x), where k: is a positive quantity which has the dimension 
of an energy and w{x) a "globally positive" dimensionless function such that limv^oow(x) = 
0. The critical constant K{{d}), where {6} stands for a set of quantum numbers, is such 
that, if k: > K{{d}), the potential admits a bound state with the quantum numbers {6}. 
The interaction energy for the state with quantum numbers {6} is then just vanishing for 

fC=JC({0})P 

Let us consider a nonrelativistic A/^-body system (no manageable calculation can be 
performed for a semirelativistic kinematics) with one-body potentials U{x) = —ku{x) and 
two-body potentials V{x) = —gv{x), both independent of the particle mass and both admit- 
ting only a finite number of bound states. The system (I6-30I )- (I6-32I ) for a vanishing energy 
gives 

N^^=NkNur-^)+CNgNv(^^, (6-40) 



N-^ = -kNrou'r^)-^/c^gNrov' (^], (6-41) 

where k^ and g^ are the critical constants for the system with N particles. The elimination 
of the ratio N /{mr^) from both equations yields the equality 

2NkNu(^^^ +2CNgNV (^--^=^ = -kNrQU (^^^ - y^gNrov' (^-^^^ . (6-42) 

When potentials u and v are both taken into account, nothing interesting can be said. So let 
us consider one type of potential at once. 

Assuming that only two-body forces are present, (16-421 ) reduces tc 



where the parameter g^ has disappeared. Introducing the new variable yo = ro/^/CN, we 
can rewrite (I6-40I ) and (16-411 ) as 

1 2 

2v(3;o)+3'ov'(jo)=0. (6-45) 

The variable y^, determined by (16-451) . is independent of T^i, Q and m, and depends only on 
the form of the function v{x). So, (16-441 ) gives precise information about the dependence of 
the many -body critical constant g;^ as a function of all the characteristics of the system. With 
the system (I6-44I) - (I6-45I) . it is easy to recover previous AFM results obtained in Sects. |4A2] 
and 143:21 and in Refs.|67l),|68]). 

Within the AFM approximation, the ground state (GS) of a boson-like system is char- 
acterized by 2 = |(A'^ — 1). We obtain in this case the following very general relation valid, 
at the AFM approximation, for all pairwise potentials with a finite number of bound states 

g^+i(GS) 
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This ratio has previously been obtained and numerically checked for several exponential- 
type potentialsP^-IS^ Similarly, in the same general situation, 

^;v(GS) = ^^2(GS), (6-47) 

indicating that in order to bind a A/^-body system, a coupling N /2 times smaller than the 
coupling for a two-body problem is sufficientP^''^ 

Assuming that only one-body forces are present, a similar calculation give^^ 

^^ = ^^^^ — ^ (6-48) 

2m(3;o)+3'om'(3'o)=0, (6-49) 

where the change of variable yo = ro/N has been used. Again, (16-481 ) gives precise infor- 
mation about the dependence of the one-body critical constant ki\i as a function of all the 
characteristics of the system. These results are strongly different from those for pairwise 
forces. 

If the AFM gives upper (lower) bounds for the exact eigenvalues, the critical constants 
predicted by formulae above are upper (lower) bounds for the exact critical constants. In the 
favorable situation where the AFM gives both upper and lower bounds for the eigenvalues, 
it is possible to approximate from above and below the exact critical constants. 

6.3. Connection with the perturbation theory 

It has been shown in Ref. I40l ) that, for one- and two-body nonrelativistic systems, the 
AFM and the perturbation theory give similar results when the potential is an exactly solv- 
able one plus a small perturbation. This result is extended here for the general Hamiltonian 

Let us first assume that each pairwise potential V(r,y) is supplemented by a term £v(r,y), 
with e <C 1 in order that ev{x) <C V{x) in the physical domain of interest. In the system 
(I6-30I )- (I6-32I ). the potential V{x) is replaced by V{x) + ev{x). In this case, new values r\ and 
pi for the mean radius and momentum will be the solution of the new system 



Ml =NT[pi)+Nu(^-^]+Cn 



(6-50) 



Pin = Q, (6-51) 



NpiT'{pi) = rrU'Q') + y^r, 



(6-52) 



Writing r\ = (1 + 5)ro, we can expect 5^1 since £ ^ 1. In this case, power ex- 
pansions at first order can be computed. We have pi ^ {I — 5)po from (I6-51I) . and 
we can write T{p:) « T{po) -dpoT'{po), T'{p,) « T'{po) -dpoT"{po), U{r,/N) « 
U{ro/N) + 5rQU'{ro/N)/N, etc. Equation (16-521) reduces to an expression of the form 
5 £h{ro) where /i is a quite complicated function of T', U', V' and their derivatives. It 
simply confirms that 5 ~ 0(e). It is then possible to perform an expansion of Mi which 
reduces, using (16-301 ) and (16-321 ). simply tc^^' 

Mi=Mo + Cn£v( +0{e^). (6-53) 
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This result could seem quite obvious, but it demonstrates that the knowledge of ro is suffi- 
cient to obtain the contribution of the perturbation at the first order. 

Let us now assume too that each [one-body potential U / kinetic operator 
is supplemented by a term [tj u{\si\) I Tf(|/?,|)], with [tj ^ 1 / T ^ 1] in order that [r\ u{x) <^ 
U (x) / Tt{x) <^ T{x)] in the physical domain of interest. With similar calculations, we finally 
find 

Ml=Mo+^^T^(po)+^^TJMQ)+C;v£v(^-^^+0(T^77^£2). (6.54) 
The parameter 5 is determined at the same order by the following relation 



NJ ^ " " VVQv 
2Npo T' ipo) +Npl T" ipo) + §U"Q)+ri V" (-^ 



(6-55) 



Perturbed observables and wavefunctions can then be computed at first order, since ri = 
(1 + 5)ro and p\ = {\ — 8)po at this order. 

The contribution of a perturbation at the first order can thus be very easily computed 
within the AFM once the unperturbed problem is solved. In order to check the quality of this 
approximation, let us consider a case in which the unperturbed Hamiltonian H can be solved 
exactly by the AFM, that is Mq is the exact solution. If the small perturbation potential is 
written £Li<y=i ^(''o)' quantum perturbation theory says that the solution is given 
by 

M, =Mo + Cw£(v(r,y)) + 0(£^), (6-56) 

for any pair {ij). The mean value is taken with a completely symmetrized eigenstate of 
the unperturbed Hamiltonian H. The comparison of (16-561) with (16-531 ) shows that {v{rij)) 
is replaced by v {tq/ \/Cn) within the AFM. This is to be compared with the exact relation 
{S{rij )) =S[ro/ \/Cn) for the auxiliary potential (see (16-121 )). So, the AFM does not give the 
same result as the perturbation theory. But the agreement can be very good, as shown with 
several examples calculated above. Similar- discussions can be made for small one-body 
perturbation potentials or small perturbations of the kinematics. 

6.4. Power-law potentials 
6.4.1. General results 

The first explicit example that will be considered below is the case of power-law poten- 
tials, i.e. the Hamiltonian 

N I N N 

^pi = IV^? + '^' + «sgn(A)£r<^+Z7Sgn(T]) £ r"^., (6-57) 

!=1 i=l i<j=i 

where A, T] > — 1 (in the nonrelativistic case, one can further consider A, tj > —2). When 
only a one-body or a two-body interaction is present, parameters a or b must be positive. 
If both types of potentials are present, there are less constraints on the sign provided that 
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a bound state can exist. Following the definitions (I6-5I ). it is readily computed that K{x) = 
a\X\x^-^/2 mdL{x) =b\ri\x'^-^/2. Then, by defining 

2-A 2-n 

A^=a\X\i^-j , B^=Z,|T]|A^(^^J , (6-58) 
Equations (16-261 ) and (16-271) can be recast under the form 



M{Xo) = N^m^ + ^Xo + Q {^^X,'/' + ^^o"'^^') , (6-59) 



X^ = \lm^ + ^Xo 



2-A 2-1) 



(6-60) 



The sufficient condition to get a closed analytical formula for M{Xq) is to solve (16-601 ) ana- 
lytically. This is possible if this last equation can be rewritten as a polynomial equation of 
the fourth degree at most. If A 7^ Tj , it can be computed that the following couples will lead 
to such an analytically solvable equation 

(A, 77) or (77, A) = (-^,-1), (0,-1), (1,-1), (-1,-2), (0,-2), (-1 -^) . (6-61) 

Note that the last three cases are only allowed with a nonrelativistic kinematics. 

The problem actually becomes particularly simple when A = 77 . In this case, the fol- 
lowing values lead to an analytical solution 

2 1 7 5 3 4 5 

■^ = -l,-2'~2'^'^'^'~^'~4'~3'~2'~3'~4" (6-62) 

The last six values are only allowed in the case of a nonrelativistic kinematics. The solution 
of (16-601) reads 



^0 + = i^^ + Ba )M "^^ + ^^0 ) , (6-63) 



and the mass formula ( 16-591 ) becomes 



MiXo) = ^^^!^^±m±2)^. (6-64) 
AJm2 + £Xo 



Among the values (I6-62I ). three cases are of obvious physical interest: A = — 1, 1, and 2. 
The case A = 2 corresponds to the harmonic oscillator and is solved in Appendix |Jl The 
Coulomb problem, i.e. X = —I, will be specifically considered in Sect. 16.4.41 while we will 
discuss the case A = 1 in Sect. 16.4.21 Notice that closed mass formulae for any A can be 
obtained in the nonrelativistic and ultrarelativistic limits, as we also show in the following. 
In the case m = 0, (16-631) allows us to extract the Xq value for arbitrary power 



1 



Q 

-( 



Xo={tM?.+Bi)^] . (6-65) 
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Inserting this value in (16-641 ) provides us with the expression of the energy in the ultrarela- 
tivistic limit 

A + 1 



M(Xo) 



X 



2(A + 1) 



(6-66) 



We already know from Sect. I5.3.4I that a bound state of massless particles exists only if 
A > 0. Setting N = 2 and a = 0, this equation reduces to (15-631 ) with a = 2 (and a renamed 
b). 

On the contrary, assuming a large mass, m, « M{Xo), it is possible to extract Xq from 
(1^631) 

Xo=[m{A^+B^)]^-. (6-67) 

Except the Coulomb case (A = —1) which will be treated subsequently, the term is 
dominant with respect to Xq. Expanding the expression of the energy (16-641) at lower order 
in 1/m leads to the nonrelativistic value of the energy 



M(Xo) =mt + 



A + 2 



2A 



-Q 



m" 



1 

A +2 



(6-68) 



Setting N = 2 and a = 0, this equation reduces to (14-71) with the reduced mass m replaced by 
m/2 (and a renamed b). 

6.4.2. The linear potential 

The case A = 1 corresponds to a linearly rising confining potential. It is of great interest 
in hadronic physics since a linearly rising potential appears to be the best way of modeling 
the QCD confining interaction within potential approaches, see for example Ref. |94| ) for 
more details. For the three-body problem, the one-body potential (term in a) corresponds 
to the so-called Y-j unction, while the two-body potential (term in b) corresponds to the 
so-called A approximation. For physical problems, one has the choice of retaining one 
approximation or the other or both. Here for sake of generality we present the general case 
including both. In the real QCD world, one must take N = 3. But a general N value is 
interesting to study alternative approaches of QCD.'^SJJSn) 

For A = 1 , one can write 



Ai+Bi = 

The solution of (16-631) when A = 1 is given by 



— c, with 
Q 



- a + bWCN- 



Xo 



where F_(F) is given by (IB -21 ) and where 



2Qc 



(6-69) 



(6-70) 



(6-71) 



Introducing this value in the expression of the energy (I6-64I ). a simple calculation allows us 
to obtain the mass under the form 



M(Xo) = Nm 



F-{Y) 
2Y 



F-{Y) 



(6-72) 
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In the ultrarelativistic limit, where m = 0, one obtains, from (I6-66I ). the very simple 
relation 

M{XQf =ANcQ. (6-73) 

Such a linear behavior of the square mass versus the principal quantum number is a well- 
known fact in hadronic physics, where light mesons and baryons are known from experi- 
mental data to exhibit Regge trajectories. 

In the nonrelativistic limit, the mass formula (16-681 ) can be recast under the form 

M{XQ)=mt + -[^—] . (6-74) 
2 \ m J 

The second term is the "binding energy", which is a positive quantity in case of a positive 
linear potential. 

In the special case of a linear potential, there is a peculiar relationship between the mass 
of the two-body problem and the mass of the A/^-body problem. More precisely, let us call 
M^^\o,m,b,Q2) the AFM eigenvalues of the Hamiltonian 



h'^^^ = ay/p^ + m^ + br (6-75) 

calculated with the natural quantum number Q2 = 2n + l + 3/2. Then, using (16-731) and its 
counterpart for A/^ = 2 with the value of the quantum number Q defined by (16-191 ). the mass 
M^^\N ,m,a,b,Q) of the original Hamiltonian is given by 

M^^\N,m,a,b,Q)=M^^^ (N,m,a + by^,Q) . (6-76) 



6.4.3. Baryonic case 

In this section, we want to discuss a special situation which can have an immediate 
application for hadronic systems, especially baryons or glueballs.^-'SS} Let us consider a 
one-body potential with a linear shape and a two-body potential of Coulomb type (A = 1 
and 77 = —1). We noticed that an analytical solution does exist, but one needs to solve a 
general polynomial of degree 3 and the corresponding solution is quite involved. We prefer 
to give here the solution in the ultrarelativistic limit (m = 0). In this special case the solution 
of (16-601) is quite simple 

and ( 16-591 ) can be simplified as follows 



M{XQ)=2^/^^QN-bc]l^. (6-78) 

Since the argument of the square root must be positive, it is readily checked that, either 
there exists a maximal allowed number of particles for a fixed value of b, either there exists 
a maximal allowed value of b for a fixed number of particles. This kind of mass formula has 
been applied to the computation of light baryon masses for various theories of QCD with a 
large number of colors .'^■'Sol 
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In order to test the relevance of our method, we apply it for a particular relativistic 
three-body system: the light baryon composed of three massless quarks. In the framework 
of constituent models, the Hamiltonian for such a system is given bjl^US 

^' = ty^ + ^.vII-'-^|-^«^ t TT^- (6-79) 

i=i i=i ^ i<j=i l'^' 01 

The dominant interaction is a confinement by three strings with density meeting at the 
center of mass. The short-range part is given by pairwise interactions of Coulomb nature 
with a strong coupling constant as- In order to allow a more realistic calculation of the 
correct baryon masses, this Hamiltonian must be completed with terms that can be computed 
in perturbation.^^ ' Since these contributions represent no interest for the AFM, they are not 
considered in this study. 

Let us denote the AFM counterpart of H^. Approximations of eigenmasses of 
the genuine Hamiltonian are given by the application of (16-781) with N = 3, a = Os and 
b = 2as/3, 

M^^^ = ^12a.(B + 3-^), (6-80) 

where B = 2{ni + «2) + ('1 + ^2) is the band number (see (1619b ). One can notice the strong 
degeneracy due to the particular form of B. This formula generalizes a result obtained for 
a small value of as in Ref. "95^. Following above considerations, these masses are upper 
bounds on the exact masses. 

In order to test the accuracy of this formula, it is necessary to compute accurately eigen- 
values of H^. It is possible, for instance, to use a variational method relying on the expansion 
of trial states with a harmonic oscillator basis.'^ We can write 

\¥)=t I (6-81) 

B=Qq{B) 

where B characterizes the band number of the basis state and where q summarizes all the 
quantum numbers of the state (which can depend on B). This procedure is specially inter- 
esting since the eigenstates of are expanded in terms of eigenstates of (up to a length 
scale factor). In practice, a relative accuracy better than 10^^ is reached with Bmax = 20. 
Such results are denoted "exact" in the following. 

Before comparing exact masses with formula (I6-80I) . it is significant to describe the 
baryon wavef unctions. Quarks are fermions of spin 1/2 with isospin and color degrees of 
freedom. The global color function is unique and completely antisymmetrical. For total spin 
S = 3/2 (isospin T = 3/2) the spin (isospin) wavefunction is completely symmetrical. The 
corresponding wave function are of mixed symmetry for 5 = l/2orr = 1/2. Using degen- 
erate eigenstates of H^, it is always possible to build three-quark states completely antisym- 
metrical,'^''^ which are characterized by the same mass. For instance, when S = T = 3/2 
or S = T = 1/2, the baryon states possess a spatial wavefunction completely symmetrical. In 
the real world, the degeneracies are removed by spin- or isospin-dependent operators which 
can be computed as perturbations.^^ ' When an eigenstate of is computed for given values 
of S, T, and the total angular momentum L, we determine the band number of its dominant 
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components. It can then be compared with the eigenstate of with the same spin-isospin 
quantum numbers and the same band number. 

In Table IXIVI exact masses are compared with the predictions of formula (16-801 ) for 
usual values of the parameters and as- The relative error can be large but the quality of 
the approximation is quite reasonable for a so simple formula. Due to the strong degeneracy 
of the harmonic oscillator, all AFM eigenstates with the same band number have the same 
energy. This approximation is better for high values of L. Let us note that the relative error 
can be reduced by a factor around 2 when as — )• 0. When as = 0, the relative error becomes 
independent of the energy scale factor 



Table XIV. Exact eigenmasses of the Hamiltonian ( l6-79b as a function of the band number B and the total 
angular momentum L for Os = 0.2 GeV^ and = 0.4. The number in brackets is the probability (%) of the 
component with the band number B in the harmonic oscillator expansion J6-81l >. These results are compared 
with the masses given by formulae J6 80b . ( l6-82t . and i6-83t . The number in parenthesis is the relative error 
(%) with respect to the exact value. All masses are given in GeV. 



B 


L 


Exact 


M^"-" (|6-80|l 




(I6-82D 


M4^^« (16-831) 








2.128 [92.9] 


2.468 (16.0) 


2.168 


( 1-9) 


2.168 (1.9) 


1 


1 


2.606 [95.9] 


2.914(11.8) 


2.596 


( 0.4) 


2.596 (0.4) 


2 





2.739 [89.4] 


3.300 (20.5) 


2.962 


( 8.1) 


2.811 (2.6) 




2 


2.959 [96.0] 


3.300(11.5) 


2.962 


( 0.1) 


2.962 (0.1) 


3 


1 


3.125 [91.7] 


3.646 (16.7) 


3.288 


( 5.2) 


3.152 (0.9) 




3 


3.299 [96.7] 


3.646 (10.5) 


3.288 


( 0.3) 


3.288 (0.3) 


4 





3.260 [80.8] 


3.961 (21.5) 


3.585 


(10.0) 


3.332 (2.2) 




2 


3.422 [92.3] 


3.961 (15.8) 


3.585 


( 4.7) 


3.460(1.1) 




4 


3.581 [96.8] 


3.961 (10.6) 


3.585 


( 0.1) 


3.585 (0.1) 


5 


1 


3.584 [86.3] 


4.253 (18.7) 


3.858 


( 7.7) 


3.625 (1.1) 




3 


3.716 [93.6] 


4.253 (14.5) 


3.858 


( 3.8) 


3.743 (0.7) 




5 


3.861 [97.0] 


4.253 (10.2) 


3.858 


( 0.1) 


3.858 (0.1) 


6 





3.721 [74.4] 


4.527 (21.7) 


4.114 


(10.6) 


3.782 (1.6) 




2 


3.838 [86.4] 


4.527 (17.9) 


4.114 


( 7.2) 


3.895 (1.5) 




4 


3.966 [93.6] 


4.527 (14.1) 


4.114 


( 3.7) 


4.006(1.0) 




6 


4.103 [96.9] 


4.527 (10.3) 


4.114 


( 0.3) 


4.114(0.3) 



We have shown many times that it is possible to improve two-body mass formulae by 
changing the structure of the principal quantum numbers Q. By fitting another form on 
exact eigenvalues, a very high accuracy can sometimes be reached. Here, we will proceed 
differently and will try to use analytical results to find the best shape for the mass formula. 
Using a good upper bound on the ground state given by the trial state of (16-811 ) at Bmax = 0, 
a modified AFM formula was proposed 

Mf™ = y^a,(B + 3- V3as). (6-82) 

One can see in Table IXIVl that the masses are greatly improved, but the variational character 
of the formula (16-821 ) cannot longer be guaranteed: some masses are now below the exact 
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ones. At last, using information coming from the WKB method,^^' the problem of the strong 
degeneracy has been improved by a new formula 



Details about these procedures are given in Ref. 68j. One can see in Table IXIVI that the 
relative error is now around 1% and generally below. Despite its simplicity and its non- 
variational character, (16-831 ) is then a very good mass formula for the eigenstates of the 
Hamiltonian (16-791 ). It is not sure that the procedure used here to improve the mass formula 
for baryons could work so well for other Hamiltonians. But, this shows that an improvement 
is possible, at least in some particular cases. 

6.4.4. Atomic systems 

We turn now our attention to atom-like systems, i.e. systems described by the following 
Hamiltonian 



For this section, we remain completely general, the only hypothesis being that the parame- 
ters a and a should be real positive numbers. The obvious application could correspond in a 
first approximation to identical electrons, with charge e, feeling the attraction of a central 
static source (the nucleus) and their own repulsion. In this particular situation, one must 
make the identification a = Ne^ and a = e^. At this stage, we point out that our formalism 
is spin-independent and thus that this model cannot be directly applied to a real physical 
atom. Moreover, one should be very careful in applying the formulae already given. Two 
cautions are in order: 

• In (16-841) . r; must be the distance between the particle / and the center of mass of the 
system of N electrons. In the atomic interpretation, this center of mass must coincide 
with the nucleus. This is not always the case but this prescription should be valid in 
case of a spherical atom. 

• Electrons are fermions and the total wavefunction (space and spin) must be completely 
antisymmetrical. Since the spatial wavefunction cannot be completely symmetrical, all 
possible values are not allowed for Q given by (IJ-171 ). For closed shell atoms, 2ags 
(IJ-191 ) should be used to estimate the ground-state energy only. 

Nevertheless, solving Hamiltonian (16-841 ) has an intrinsic interest since, to our knowl- 
edge, no corresponding analytical mass formula is known so far. A general solution can 
moreover be found with the relativistic kinematics, starting from (16-631 ) and (16-641 ). Apply- 
ing them to the case X = —\,a = a,b = —a and defining 





(6-84) 




(6-85) 



the solution of equation (16-631 ) is given by 



(6-86) 



^0 



2(1 
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The mass formula (16-641 ) then reads 

M{Xo)=mty/l-D^. (6-87) 

In the ultrarelativistic Umit, this expression gives the value M = 0. This property is well 
known in the two-particle system and prolongates to the A/^-body problem. It is quite obvious 
that a system of massless particles cannot be characterized by a (negative) binding energy 
since the resulting mass would be negative. Moreover, with a Coulomb-type potential, the 
only energy scale of the problem is the particle mass. So, if this mass vanishes, no bound 
state can exist. In the nonrelativistic regime, one must assume that the D quantity is small 
in (16-871 ). In this case 

M"' = mt-^m,D^. (6-88) 

It is clear from formula (16-871 ) that the value of the parameters cannot be arbitrary. One 
must satisfy the condition D < I. Explicitly this means 

aN^-aC^^^ <QN. (6-89) 

Let us emphasize that this condition is independent of the mass of the particle.'^ Let us 
notice that, for a = Ne^ and a = e^, a. maximal value of A'^ exists if Q increases less rapidly 
that A^2. 

6.5. Duality relations 
6.5.1. Generalities 

In this section, we will show that AFM is a very powerful method to give relationships 
between the energies of states (excited or ground states) for two different systems. We call 
the relations between the energies of both systems "duality relations". The case for which 
one of the two systems contains 2 particles is specially interesting since the corresponding 
eigenenergies are rather easy to obtain. These relations are exact if we consider the energies 
obtained in the AFM; they are no longer exact for the true eigenstates, but if AFM gives a 
good precision on the true states, we hope that the duality relations are good approximations 
of the physical situations and gives, at least, general trends for the exact levels. A rather 
complete treatment of duality relations can be found in Ref.|99b. We present here a simplified 
version where we restrict our discussion to systems of N identical particles interacting only 
through a two-body potential V{r). Readers interested in more specific details may have a 
glance to Ref. 99) for further informations. 

Our original Hamiltonian is thus 

i=l i<j=l 

In this expression r, is the position operator for particle /, /?, its conjugate momentum and m 
the common mass of all particles. The AFM mass of the system is given by (16-261 ) which is 
written now as 

M(X„)^Nf^„ + C,v(.r^], (6.91) 
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while the Xq quantity is solution of the transcendental equation (16-271 ) which turns out to be 
now 



and where L{x) = V'{x) / {2x) as usual (see (I6-5I )). 

In general, the potential V{r) depends on physical parameters {ti , T2, . . . , Tp} = {t} 
and should be noted more precisely V{{x};r). The other parameters of the problem are the 
number of particles N and the mass m of the particles. Lastly, we want to have a description 
of the whole spectrum, so that the principal quantum number Q also enters the game. In 
consequence, a complete notation for the eigenmasses would be M^^\{T};m, Q). 

Let us stress now a very important point: for realistic potentials, the parameters {t} 
could depend on N or/and m. We do not consider such a particular behavior here. Thus we 
assume that the parameters of the potentials are independent of and m. This means 
that the m dependence of H is only through the kinetic energy and its A'^ dependence through 
the numbers of terms in the summation. We suppose that the potentials are given once for all 
and consider that their form and parameters do not vary for all studied systems. Therefore 
we are finally interested only by the N, m, Q dependence of the eigenmasses and use the 
simplified notation M^^' {m, Q). What we call a "duality relation" is just a relation between 
M^^\m,Q) and M^P\m',Q'). 

6.5.2. General case 

Instead of the Xq quantity, let us define the new variable 



Then, the transcendental equation (16-921) and the mass expression (16-911) are replaced by the 
following ones 

slJl+CNimso/QfV'iso), (6-94) 



{N-1)VCn 
M^^\m,Q)=CN 



2Q 1 

{N-\)VC^so 



\+CNimso/Q)^ + Viso) 



(6-95) 



A little algebra on those equations allows to arrive at the following duality relation 
between the A^-body and the p-body systems with particles interacting via the same two- 
body interaction 



In this case, the spectrum of the A/^-body system is the same as the /j-body system (with the 
same two-body potential) provided we consider different particle masses in both situations 
and different excitation states. 
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6.5.3. Ultrarelativistic limit 

The case of ultrarelativistic systems, characterized by a vanishing mass m = 0, presents 
some very specific and interesting features. The case of systems composed of gluons or/and 
light quarks can be well represented in this scheme. The Hamiltonian of the system is then 

H = i^/pf+ I V{\ri-rj\). (6-97) 

i=i i<j=\ 

Indeed, the formulation is simpler for this particular situation. Putting the value m = in 
(16-911 ) and (16-921) . one gets a new set of equations. The transcendental equation looks like 



^xr- = 2NLLU0], (6.98, 



while the corresponding AFM eigenmass is given by 



MuiXo) = ./NQXo + CnV • (6-99) 

The eigenmass, depending now only on N and Q, will be noted M^^^ (Q) (the index u stands 
for "ultrarelativistic"). 

With the same definition ( 16-931 ) for sq, the transcendental equation ( 16-981 ) reduces to 

22 2,„ 



Defining the D{x) and F{x) function as in (15-271) and (15-281) calculated with the potential 
V{x), one has 



As we already mentioned, the D and F functions are universal. From expression (16-1021) . 
one deduces immediately the duality relation 

which is the special case of (16-961) with m = 0. 

6.5.4. Nonrelativistic limit 

Another interesting limit of the theory is the nonrelativistic one, valid when the mass of 
the particles is large compared to the mean potential. In this case, the considered Hamilto- 
nian is simply 

^ = l||+ £ yi\ri-rj\). (6-104) 

i=l i<j=l 
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Instead of dealing with the total mass M^^\ it is better to consider the binding energy 
obtained by removing the total rest mass: E^^'> = M^^^ —Nm. The AFM approximation 
E^'^\m, Q) of the binding energy is given by the following equations, 



2Q 



and 



.,«(x„).e^,c..y^j. (,,06, 

We will see that, in this nonrelativistic limit, one obtains additional interesting properties. 
Introducing by ( I6-93I ) as above, the transcendental equation (16- 1051 ) reduces to 



slV'iso). (6-107) 



Defining the ^(jc) and ^{x) function as in (14- 111) and (14-131 ) calculated with the potential 

V{x), one has 

^o = ^(^), (6-108) 
£(^)(m,e)=C^^(^). (6-109) 

As the case studied in the previous section, the ^ and ^ functions are universal. For a 
nonrelativistic system, an additional property appears. Remarking that E^^^ depends only 
on the ratio Q^/m, one has the general duality relation 

£W(m,e)=£W(j8Vi8e), (6-110) 

valid for any value of the real parameter j8. 

From expression ( 16-1091 ). one deduces immediately the duality relation 

(6.in, 

which is identical to the general case (16-961) . One can use the property (16-1101) to obtain 
many other possibilities. For example, choosing the value jS = a/ p{N — \)/{N{p — 1)), one 
has the alternative simpler duality relation 

£W(m,e) = ^£(^)(^/n,^e). (6-112) 
In this last relation, let us choose j3 = ^jWpp, one arrives at the relation 

E^^\m,Q) = L^J?-q\ . (6-113) 
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In this expression, we decide to keep the same mass for systems with N and p particles. The 
duality relation leads to a link between different excited states of both systems. Choosing 
the value j8 = Cn/Cp in equation (16- 1121) . one obtains the alternative expression 

In this expression we decide to maintain a one to one correspondence in the spectrum but 
for systems with different particle masses. 

6.5.5. Passing from nonrelativistic to ultrarelativistic limits 

We showed that both the ultrarelativistic limit and the nonrelativistic limit for the eigen- 
masses share the property of being expressed in terms of universal functions. The F function 
for the ultrarelativistic case and the ^ function for the nonrelativistic case are independent 
of the system and of the excitation quantum numbers, but depends only on the form of the 
potential under consideration. 

If the same form of potential is used in both situations, the expressions of the F function 
and of the ^ function are not identical so that the corresponding spectra are quite different. 
However, we showed in Ref. ,99) that if the potentials are different but Unked by a certain 
relationship, one can arrive at very interesting conclusions. 

We report here the main conclusions, skipping all the rigorous proofs that can be found 
in the previous reference. Let us consider a system described by a nonrelativistic treatment 
based on a two-body potential V{r) whose binding energy is E^^\m,Q). One knows that 

E'^^\m,Q)=CN^{R) with /?=^. (6-115) 

If one uses rather the potential W{r) = V{a^/r) , one has 

E^'^\m,Q)=CNF{S) with S=:^. (6-116) 

If we take the special value S = 2Q/{{N — 1)\/Cn), the value CnF{S) represents the 
ultrarelativistic mass Mi,^\q) of the same system but obtained with the potential W{r). The 
condition on S and the link between S and R leads to the condition defining the value of a, 
namely 

Q = 2m^/C^a'^. (6-117) 

In consequence, one can state the following theorem: 

If is the binding energy of a nonrelativistic system governed by the 

two-body potential V{r) and if Mf\Q) is the mass of the related ultrai^ela- 
tivistic system governed by the two-body potential W{r) defined by W{r) = 

V (^\/Gr/(2mVCW)) , then one has the general property Mf^ (Q) = (m, Q). 

From its definition, the potential W{r) depends on m so that the notation Mif^\Q) which we 
have used up to now, depends indirectly on m, as imposed by the theorem. 
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6.5.6. Applying duality relations to exact levels 

All the duality relations presented above are exact for the AFM solutions of quantum 
systems. One can wonder to what extent these constraints are satisfied for the corresponding 
exact solutions. We discuss briefly this point and, again, refer the reader to Ref. i99t for 
additional information. We are interested essentially to the case of nonrelativistic systems. 

Let us call £(^)(m; {«,■,/,}) the exact eigenenergy for state of the Hamiltonian (16-901 ) 
with {«,•,/,} = {ni,li,n2,l2, ■ ■ ■ ^'^w-i j^A^-i}. while e{m;n,l) = e^^\m;{n,l}) gives the ex- 
act spectrum of the corresponding two-body problem. The first thing to do is to connect 
the exact level £(^)(m; {«,,/,•}) to the corresponding AFM approximation E^'^\m,Q). Of 
course, one can always impose an equality between both 



This equality defines in practice the value of the principal quantum number Q that must be 
used in the AFM solution. The problem is that the value obtained in that way does depend 
explicitly on the mass m, in addition to the basic quantum numbers {«;,/,}. This is contrary 
to the philosophy of the AFM approach for which Q depends on {«, ,/,} only. To apply 
AFM formulas, one must give up the dependence on m for Q and choose a dependence on 
quantum numbers only, With such a constraint, the equality (161181 ) does not 

hold anymore and is replaced by an approximate value 



The cleverness of the physicist is to guess the form for Q^'^-' ({«,•,/,•}) which makes the pre- 
vious approximation as precise as possible. The choice ( 16191 ) is the most natural one since 
it follows directly from the construction of the AFM levels. However, the corresponding 
approximation is not always good and, moreover, it is plagued by degeneracies which are 
absent in the real physical spectrum. A choice like the one proposed above ( 16-391 ) could be 
much better in this respect, but very often the determination of the parameters entering it 
needs a preliminary study of the problem. Let us suppose that we have the technical ability 
to choose a practical and correct form for the principal quantum number. 

The idea of the method is to approximate the exact values by the AFM ones and to take 
the exact duality relations on the AFM results to extend them to the true states. This can be 
done in two steps 

1 . the link between the excited state of the A/^-body problem to the ground state of another 
A/^-body problem; 

2. the link between the ground state of the A^-body problem to the ground state of a two- 
body problem which can be solved easily. 

These two steps are the consequence of the duality relations given above and their proofs 
can be found in Ref. 99). Here we focus essentially on the practical result. 

The two-body problem can be solved exactly (at least numerically) for two particles of 
mass m interacting through the potential V{r). Thus, given the potential V, the ground sate 
energy is function of the mass only. Let us define the / function by 



£W(m;{n,,/,})=£W(m,(2(m;{n,-,/,})). 



(6-118) 



(6-119) 



/(m) = £(m;0,0). 



(6-120) 



This function is universal in the sense that it depends only on the form of the potential V (r) 
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and it can be computed once for all. An example of the function /(m) for the Hamiltonian 
H = p^/m + r is given in Fig.|5j 

12 - 
10 - 
8 - 

/(™ ) 



Fig. 5. Universal function /(m) as a function of the mass m for tlie ground state of the two-body problem with 
V{r) = r. 

The conclusion of the two steps indicated previously can be summarized by the follow- 
ing approximation 



£(^) (m; {nuk}) « CNf{M{m- {mJi}) 



(6-121) 



with the following expression for the mass M entering as the argument in the / function: 

2 

I . (6-122) 



The conclusion of the relation (16-1211 ) is very strong. It proves that the whole spectrum 
of all systems can be obtained approximatively by the calculation of a universal function, 
f{m), corresponding to the ground state of the 2-body system with the same potential, and 
for arguments M(m; {«;,/,}) given by (16-1221 ). Getting /(m) is a very easy task. Solving the 
2-body system can be performed with a great accuracy for any potential (let us recall that 
this potential must not depend on m and N); moreover, obtaining the ground state energy 
is free from possible numerical complications arising for excited states. Let us note also 
that the duality relation (16-1211 ). which is obviously an approximation, concerns the exact 
eigenvalues only; the AFM values, which were very convenient intermediate quantities in 
our demonstration, have completely disappeared. The accuracy of this duality relation for 
exact states has been tested with success in Ref.l99l) for some systems. 



57. Conclusions 



In this paper, we presented a detailed review of the auxiliary field method (AFM) to 
solve eigenequations in the context of quantum mechanics. We showed that this method 
is very powerful to obtain closed analytical approximate expressions for the properties of 
many physically interesting systems. 
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The idea of the method is to replace a Hamiltonian H for which analytical solutions 
are not known by another one H for which they are known. For instance, a potential V{r) 
not solvable is replaced by another one P{r) more familiar, or a semirelativistic kinetic part 
is replaced by an equivalent nonrelativistic one. The bridge between both Hamiltonians is 
ensured by a special function including one or more auxiliary fields. These fields can be 
determined by an extremization procedure and Hamiltonian H reduces to Hamiltonian H. 
In so doing the original problem is completely equivalent to the new one. The approximation 
comes from the replacement of the auxiliary fields by pure real constants. The approximant 
solutions for H are then obtained by the solutions of H in which the auxiliary parameters 
are eliminated by an extremization procedure for the eigenenergies. 

The great advantage of this method is that it is able to predict the behavior of the observ- 
ables in terms of the various parameters entering into the problem and also in terms of the 
quantum numbers. We presented the principles of the method, its characteristics, its major 
properties and we discussed the quality of the results for many relevant situations appearing 
in various domains of physics. The only restriction is the fact that, for relativistic systems or 
many-body problems, the computation is manageable only for systems composed of identi- 
cal particles. 

Although we focused essentially on the search for closed formulae of the eigenenergies, 
we showed that the method is also able to provide good results for the eigenstates and for 
different types of observables. In this context the harmonic oscillator wavefunctions and the 
hydrogen-Uke wavefunctions are particularly important. As an example, we proved that the 
AFM is very simple to implement for the Schrodinger equation with a linear potential and 
gives very satisfactory results. 

The main part of the paper deals with the approximate values coming from the AFM 
concerning the binding energies of the Schrodinger equation (nonrelativistic) and the eigen- 
masses of the spinless Salpeter equation (semirelativistic). We proved a number of very 
important and pleasant general features of the method: 

• The AFM is essentially a kind of mean field approximation. In particular, using a 
power-law form for P{r), one can define a "mean radius" vq and a "mean momentum" 
Pq which allows to express the final AFM eigenvalue in the very simple form T{pq) + 
V(ro), where r(|p|) and are respectively the kinetic and potential part of H. 
A very simple link exists between po and tq. This last parameter is the solution of a 
transcendental equation which looks like a semiclassical version of a generalized virial 
theorem. 

• Even in complicated situations (for example A'^-body systems with relativistic kine- 
matics, and presence of both one and two-body interactions), the AFM can always be 
brought to the resolution of a transcendental equation, a procedure quite easy from the 
numerical point of view. In many cases, this transcendental equation can be solved 
analytically and the corresponding results have been presented in this paper. 

• We discovered, and demonstrated, an amazing property of the AFM solutions that we 
called "universality of the form". This means that, provided a power-law potential is 
chosen for the function P{r), the AFM gives an universal expression of the resulting 
energies for any potential V (r) under consideration, whatever the exponent chosen for 
power-law potential. The only remnant of this exponent is entirely contained in the 
form of a principal quantum number Q{n,l). This property of universality is valid not 
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only for the eigenvalues but also for the parameter ro and is maintained passing from 
a nonrelativistic treatment to a semirelativistic one. 

• The AFM is fully equivalent to the envelope theory (ET) and, as such, it can benefit 
of all the properties already demonstrated in the framework of this method which was 
introduced a long time ago. In particular, in many occasions, one can deduce that the 
AFM energies are upper or lower bounds on the exact energy depending on the con- 
vexity of a certain function linked to the potential. As in the ET, the energy function 
(potential and/or kinetics) appearing in H is, in a sense, the best possible energy func- 
tion tangent to the exact one. But, in the AFM, we also pushed the theory in tracks that 
have not been explored by ET. 

• If we consider only the AFM approximation for the energy, we were able to prove 
several exact duality relations between different types of problems. For instance, a 
nonrelativistic problem with a potential V{r) can be reduced to the ultrarelativistic 
(m = 0) problem with a related potential W{r) = V{a-s/r), and the mass of a A^-body 
system with identical particles is related to the mass of a two-body system with a 
simple renormaUzation of the parameters. Moreover, these duality relations are shown 
approximately verified for exact solutions. 

The 2-body Schrodinger equation with a power-law potential was the prototype for the 

starting potential P{r). The resulting AFM energies are particularly simple. Using the "uni- 
versality of the form" we proposed expressions for the principal quantum numbers Q{n,l) 
which allow to improve drastically the quahty of the results. In some cases, a sophisticated 
expression for Q gives a relative accuracy of the order of 10^ over tens of the lowest sates 
of the spectrum, while even a very simple prescription is already able to give an accuracy of 
the order of 10^^. Since the AFM can be used recursively, the power-law potential is very 
often chosen as the basic potential P{r). 

In the framework of a nonrelativistic approach, we applied the AFM to a great number 
of 2-body problems with potentials of various forms: sum of powers, square root, exponen- 
tial potentials. In each case the AFM reproduces in a very simple way the basic properties of 
the solutions, even when the starting interaction P{r) is very far from the genuine potential 
V{r). We were able to give a closed form for the energies of two particles interacting via a 
funnel potential. This situation is quite realistic in hadronic systems and, to our knowledge, 
it is the first time that an analytical expression for the corresponding eigenenergies is given. 

The AFM gives a very general formula for the critical coupling constants of nonrela- 
tivistic Hamiltonians with a finite number of bound states. The dependence on the quantum 
numbers, the mass m of the particles, the number A'^ of particles, and the structure of the 
potential are predicted. Different A'^ behaviours are obtained depending on the one-body or 
pairwise character of the interaction. If the AFM gives upper (lower) bounds for the exact 
eigenvalues, the critical coupling constants predicted are upper (lower) bounds for the exact 
critical coupling constants. 

The 2-body Salpeter equation is also tested for a variety of potentials: power laws, 
square root, fuimel. We showed that, at the Umit of a large mass, the nonrelativistic ex- 
pression is recovered in each case. The ultrarelativistic limit leads to particularly simple 
expressions. The expression for the funnel potential is really astonishing of simplicity in the 
case of massless particles. 

Once a problem is solved within the AFM (quantities po and ro found), it is very easy 



94 



to compute the contribution of a small perturbation at the first order. It is given by the 
perturbation Hamiltonian evaluated at the mean momentum po for a kinetic energy or at a 
function of the mean radius ro for a potential. The result does not coincide with the one 
obtained by the quantum perturbation theory, but the agreement can be very good. 

Lastly the AFM is able to provide analytical results for the A/^-body problems even for 
quite sophisticated types of potential: Coulomb, linear, funnel interactions are presented. 
Although the accuracy would probably become worse and worse with increasing values of 
N, we showed that the spectrum of a realistic 3-body system can be reproduced with a 1 % 
accuracy if a good definition of the principal quantum number Q is adopted. Let us mention 
also that the AFM A/^-body results were recently applied to the computation of light baryon 
masses for various theories of QCD with a large number of colors .'^■'^Si 

In summary, the auxiliary field method is easy to use, can be applied to a great varieties 
of problems in quantum mechanics regardless the number of particles, exhibits remarkable 
properties and is able to provide closed analytical expressions of sometimes very impressive 
accuracy. 

Appendix A 

The envelope theory 

A. 1. Presentation of the envelope theory 

As the AFM, the envelope theory (EXj^SI is a method aiming to get approximate ana- 
lytical energy formulae from an arbitrary Hamiltonian. Initially, it has been introduced to 
obtain bounds on the eigenenergies of the A^-body problem (see Ref. 86 )). Later, it has been 
explored and refined to simpler systems and to other purposes. We only present here its key 
features and refer the reader to Refs. l45] |. II00l l. lI01l l. [T02i) for a detailed discussion about the 
basis and applications of EX. 

Let us set V{r) = v f{r) in Hamiltonian (12-11 ). Then the energy spectrum of this Hamil- 
tonian is formally given hy E = F{v), where the dependence on the usual quantum numbers 
n and / will be dropped for simplicity. The function F{v) is concave but not necessarily 
monotonic. This allows to define a so-called kinetic potential k{S) by using the Legendre 
transformation (here the prime denotes the derivative with respect to v) 

k{S)=F'{v) and S = F{v) -vF\v). (A-1) 

This transformation can be understood as follows. \W) being the eigenstates of Hamilto- 
nian (12-11 ). one can define S = {W\ T{p^) \ W) and rewrite formally the energy spectrum as 
F{v) = S + v {W\f{r) \ W) = S + vk{S). The transformation (lA-ll ) follows from these rela- 
tions. One is consequently led to the exact formula 

E =F{v)=rmn\S + vkiS)]. (A-2) 

5>0 

What can now be done to go a step further in FT is to assume that V{r) = g{P{r)), 
where P{r) is a potential for which the solution of the eigenequation 

[T{p^) + vP{r)\ \Wa) = EAiv) \Wa) (A-3) 

is analytically known. Then, 

s={Wa\T{p^)\^a) (A-4) 
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can be analytically computed. It can moreover be shown that the kinetic potential corre- 
sponding to V{r), namely K{s), is given approximately by 

K{s)^g{kA{s)), (A-5) 

where kA{s) is the kinetic potential associated to P{r). One then obtains an approximate 
form for the eigenenergies that reads^ 

E^^ = min[s + g{kA{s))]. (A-6) 

.v>0 

The variable s actually plays the role of a variational parameter. But, thanks to (IA-31 ). the 
following equalities hold 

eA{v)=s + vkA{s), 4 (v) = ) , (A-7) 
and another approximate energy formula coming from the rewriting of (IA-61 ) is 

^ = min[eA{v)-ve^{v)+g{e;,{v))] . (A-8) 

This last formula is called the principal envelope formula in Refs. ll0l1 ). ll02l l. 

It is possible to understand (IA-81 ) as follows. If V{r) = g{P{r)), with g{x) a smooth 
function of x, then we can define the "tangential potential" V'{r) at the point r = t as 

V'{r)=a{t)P{r)+g{P{t))-a{t)P{t) 

with «(0 = ^=/(/'(0). (A-9) 

Such a particular form is obtained by demanding that V'{r) and its derivative agree with 
V{r) and V'{r) at the point of contact r = If £ ^ 1, one has indeed 

V{t + e)-V'{t + e)= ^lp'{tfg"{P{t)) + 0{e'). (A-10) 

The eigenenergies of Hamiltonian = T(p^) +V*{r), denoted by (S'{t), are given by 

^(0 = SAiait)) +g{p{t))-a{t)P{t). (A-11) 

Let us now set 

t = a^\v). (A- 12) 

It can be computed from (IA-91 ) that ^"^(v) = P"^(A(v)) with A{v) = g'^^{v), and (lA llI) 
becomes 

<f(v)=£A(v)+g(A(v))-vA(v). (A-13) 
The final energy spectrum has to be extremized with respect to v, so we have 

dv<^{^)l=v, = ^ A(vo) = £;(vo) (A-14) 
and the physical energy reads 

^(vo) = eA(vo) +g(£;(vo)) - vo£;(vo), (A-15) 
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that is nothing else than the principal envelope formula (IA-81 ). 

We have just shown that ET can lead to analytical approximate energy formulae, namely 
(IA-61 ) and (IA-81 ). which are both equivalent. Moreover, it has been shown in Ref. |45] ). as 
suggested by (lA lOl ). that (f is a lower (upper) bound on the exact energy if the function g 
is convex (concave), that is if g" > {g" < 0). Let us note that H and H' have the same 
kinetic part. The tangential potential indeed always underestimates (overestimates) the exact 
potential in this case. A clear interest of ET is thus that it allows to know the variational or 
antivariational nature of the approximation that is performed. In practice, this information 
can be obtained only for a nom-elativistic kinematics, since it is necessary to know exactly 
EAiv) in dMll- 

A. 2. Equivalence between AFM and ET 

The similarity of the starting points of ET and the AFM is obvious: In both cases, a 
potential for which no analytical solution is known is "approximated" by an other potential 
for which analytical solutions exist. It suggests that a connection between both approaches 
should exist; and it will indeed be established in this section. Let us apply the AFM as 
described above with V{r) = g{P{r)). We find the following expression for the energy (12-111) 

E{v) = eA{v)+g{P{J{v))) - vP(7(v)), (A-16) 

the function J{x) = K^^{x) being computed from the relation (12-31 ). Remarkably, this AFM 
formula is equal to the ET one (IA-131) since J{x) = P^^ {g'^^ (x)). Consequently, the AFM 
and the ET lead to the same final energy formula (IA-61 ). The link between both approaches 
is given by 

v = a(0. (A-17) 

Moreover, with the point ro defined by the relation ro = J{vq), the potential V{r) takes the 
form 

V{r)=K{ro) (P(r)-P(ro))+y(ro). (A-18) 

It is then easy to see that V'(ro) = V{ro) and that V'{ro) = V'{ro). So, the potential V{r) is 
tangent to the potential V{r). An explicit example is presented in Ref. I46l) . 

The function 7(x) can be defined if the function K{x) can be inverted. In order to 
fulfill this condition, it is sufficient that K{x) is monotonic, that is to say that K'{x) has a 
constant sign. But, from the definitions above, we have K{x) = g'{P{x)), which implies that 
K'{x) = g" {P{x))P' {x). Since K{x) must be monotonic, the convexity of the function g is 
well defined if P{x) is also monotonic. This is the case if P{x) is a power-law potential, for 
instance. In these conditions, the convexity of the function g can also be used to determine 
the variational character of the AFM. 

Let us summarize our results. The auxiliary field v can be introduced as an operator 
in the Hamiltonian (12-11) . and leads to an equivalent formulation of this Hamiltonian. If one 
considers it as a variational parameter rather than an operator, as in the AFM, the results 
are approximate but can be analytical. We have shown in this section that the auxiliary 
field, when seen as a variational parameter, is nothing else than the function a{t) generating 
the tangential potential in ET. This shows that, although obtained in different ways, the 
AFM and the ET lead to the same results. In this way, some formulae about the power-law 
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potentials obtained in Ref . 'IDS') by the ET were rediscovered with the AFM in Ref. [38] ). but 
supplementary results are given in this last reference. 

Taking this equivalence into account, we can now better understand the meaning of the 
variational parameter v in the ET: Its optimal value can be seen as given by a mean field 
approximation since Z(vo) = (Z(g'(f (ro)))) (see ( I2-19I )). Moreover, the properties of the 
AFM that have been proven in Refs.l38]).l40l) also hold for ET. Finally, we can now have an a 
priori knowledge of the (anti)variational nature of the AFM energy formulae provided that 
we express V{r) as g{P{r)) and compute whether g is convex or concave. Equivalently, since 
the potential V{r, Vq) is tangent to the potential V{r) at r = rp, the approximation ^(Vo) is an 
upper (lower) bound on the exact energy if V{r, Vq) > V{r) (y{r, Vq) < V{r)) for all values 
of r.ESt AH these results holds provided that H and ^(Vq) have the same kinetic part. Several 
examples are presented above. 

Appendix B 

Reduced equations 

Finding analytical energy formulae for the potentials that we study in this work requires 
an analytical knowledge of the roots of particular cubic and quartic polynomial equations. 
In each case of interest, algebraic manipulations allow to transform the original equations 
to one of the following reduced equations. We sum up these equations in this appendix and 
put their roots in a form that is as convenient as possible to deal with. 

B . 1 . Third order equation 

We begin by the cubic equation (J > 0) 

x^±3x-2Y = 0, (B-1) 

for which there exists only one positive root given by 

F±{Y) = (Y + VY^y^\(Y + ^Y^y^^\ (B-2) 

Written in the above form, it seems that F^(Y) is not properly defined for F < 1. But, for 
this range of Y values, one can show that 

=2cos QarccosF^ . (B-3) 

So F (Y) is well defined for all positive values of its argument. It can be checked that the 
following approximate forms hold 

2Y Y 

, F4Y)^V3 + - if F<1, (B-4) 

F±(F) (27)'/^ if F>1. (B-5) 

B.2. Fourth order equation 

The quartic equation which gives the most pleasant form for the roots is (Y > 0) 

4x'^±Sx-3Y = 0. (B-6) 
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There exists only one positive root given by 



1 



with 



1/3 



-1/3 



V{Y) = (l + ^jA + Y^j -7(2 + ^4 + 73 
The following approximate expressions can also be useful 



Y 

if F»i. 

Appendix C 

Lambert function 



(B-7) 
(B-8) 

(B-9) 
(B-10) 




Fig. 6. Plot of the two branches of the Lambert function, namely Wq(z) (solid line) and VK_i (z) (dotted line). 



Let us 



briefly recall some points concerning the Lambert function (also called Omega 
function or product-log), that we will denote W{z)- Complements of information can be 
;f. '104). First of all, W{z) is defined as the inverse function of ze~. Consequently, 



found in Ref. 104). First of all, W{z) is defined as the inverse 
it has the following properties: 

W(z)e'^W=W(ze^)=z, 
W{z) 



aw(z) 

But, it is readily observed that the inverse 



(C-1) 
(C-2) 



z[l+W{z)]' 
function of ze^ is multivalued. Two branches of 
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the Lambert function thus exist El respectively denoted as Wo{z), defined for z > — 1/e, 
and W-i{z), defined for — l/e<z<0 They are plotted in figure [6l Obviously, both 
branches of the Lambert function share the same properties (IC ll ) and (IC-21 ) as they 
meet in 1/e) = Wo(— 1/e) = —1. It can moreover be checked that Wo(|:t| -til) ^ x, 

limv^oW^-i W = and lim^^oo Wo(x) = oo. 

For our purpose, it is worth mentioning that the following equation in which a, b, p are 
real numbers such that a z. + b is positive 



is analytically solvable with the Lambert function. One has 



b 

pW 

a 



ap 



(C-3) 



(C-4) 



where either Wq or W_i has to be chosen following the range of Q and z that the function 
z{B) has to cover. 



Appendix D 

Extended virial theorem 



In this section and in the following ones, we will denote {O) = {nj\0\n,l) and 
{n\0\n) = (?i,0|C?|«,0), \n,l) being an eigenstate of H with a radial quantum number n and 
an orbital quantum number I. If A is an arbitrary operator, it follows that 

([//,A])=0, (D-1) 

due to the hermiticity of H. Note that the mean value is calculated on the exact eigenstate of 
H. If pr is the radial momentum (with [r,Pr\ = i) and /(r) an arbitrary function depending 
on r = |r| only, the computation of (ID ll) with A = Prf{r) yields to the following relation, 
called the extended virial theorem 

(lEf'ir) - V'{r)f{r) - lV{r)f{r) + ^ " ^ (^) ') = 0, (D-2) 

One recovers the usual virial theorem (see Ref. 47)) for the special choice /(r) = r. But 
new interesting relations can be obtained for other choices. In particular, if /(r) = r*+^ the 
previous equation becomes 

2(^+ l)E{/)-2{s+\){/V{r)) - {r'+'V'ir)) + 4^ " 1 -4/(/ + 1)) {r'-^) = 0. (D-3) 
Finally, if, in addition, V{r) = sgn(A)ar'*' for A 7^ 0, (ID-31 ) reduces to 

2(5 + !)£(/) -«sgn(A)(2^ + A + 2)(/+-^-) + ^{s^- 1 -4Z(/+ 1)) {r'-^) = 0. (D-4) 

4wi 

This recurrence relation is particularly useful to compute {r^) mean values for potentials 
with integer power, choosing s as an integer. 

Notice that the two branches Wq{x) and W_i (x) of the Lambert function W{x) are known by the software 
Mathematica package as ProductLog [0 , x] and ProductLog [- 1 , x] respectively. 
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Appendix E 

Observables with harmonic oscillator functions 

The eigenenergies of the harmonic oscillator HamiltoniarP 



are given by 



H=^ + vr' (E-1) 

2m 



2v 3 
E = \—Qho with QH0 = 2n + l + - (E-2) 
V m 2 



and the normalized eigenvectors by 



with A = (2mv)^^^. La is a Laguerre polynomial and 7^ a spherical harmonic. At the origin, 
the S-states are such that 

Mean values (r*) can be computed by performing directly the integrals or by using (ID-41) . 
One obtains 

, ^ 1 r(^ + / + 3/2) " r(/+p+^+(^+3)/2) 

n\ ''"''T{p + l + 3/2)r{q + l + 3/2y ^^'^^ 

where is the usual binomial coefficient. One can also write, with L = / (/ + 1), 
1 4r(?i + 3/2) 2| n 

1 8(4n + 3)r(n + 3/2) „ 1 / . , 3 



{n\r^\n) = j^^ ^, (nj^nj) = [6Qj,o -2L + - j . (E-6) 

Using the Fourier transform of a harmonic oscillator, it is easy to show that 

{p')=X^'{r^). (E-7) 

Appendix F 

Observables with hydrogen-like functions 

The eigenenergies of a hydrogen-like HamiltoniarP 

H = ^ , (F-1) 

2m r 



are given by 

mv^ 

p J with Qc = n + l+\ (F-2) 



20 
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and the normalized eigenvectors by 



with = + / + 1) and rj = mv. At the origin, the S -states are such that 



Mean values (r^) can be computed by performing directly the integrals or by using (ID-4b . 
One obtains 

^ jn + l+lf-' {n + 2l+l)\ " {p + q + k + 2l + 2)\ 

^ ' 2{2riY n\ J^J ' ^«^« (p + 2/+ 1)!(^ + 2/+ 1)! ' ^ 

One can write also, with again L = / (/ + 1), 

1\ T] / 1 \ 2t]2 



3 



('-) = ^ (32^ - L) , (r^) = {5Ql - 3L + 1) , 

= (352^ + 5Ql{5 - 6L) + 3L(L - 2)) , 

= ^ (63ec + 352^(3 - 2L) + 5L(3L - 10) + 12) . (F-6) 
Using the virial theorem and the square of the Hamiltonian, it is easy to show that 



rf ^ Sn + 2l + 5 



Appendix G 

Observables with Airy functions 



The eigensolutions with Z = of the Hamiltonian 

H=i—+ar (G-1) 
2m 

are analytically known in terms of the Airy function Ai.l^ The eigenenergies can be written 
in terms of the (negative) zeros a„ of this function, namely 

£ = -(^j ««, (G-2) 
and the normalized eigenvectors Ynoi^) = (''I") are given by 



K 

Yno{r)= r—,].,, Ai(K:r + «„), (G-3) 
/4;r|Ai'(a„)|r 
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with K = {Imayl^ . An approximate form for a„ is given b 
«« = -i8„(l + ^i8-^-^Ar' + 0(i8-'^)) with A, 
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2/3 

, (G-4) 



the series converging very rapidly with n. At the origin, the square of the wavef unction 
reduces to 

IV^«,o(0)|^ = g. (G-5) 

The remarkable fact is that it does not depend on the radial quantum number. This property 
is specific to the linear potential. 

Mean values (r*) can be computed by performing directly the integrals or by using 
(ID-41) . One obtains 

/ I I V 2 1 I / I 2 1 \ 8 1 I ^ 

[n\r\n) = — — , [n\r \n) 



3 16|a„p + 15 , I 4, s ^ ,8|a„|'^ + 25|a„| 
{n\r\n) = ^^^, {n\r \n) = 16 ^^^^^ . (G-6) 

Using the virial theorem and the square of the Hamiltonian, it is easy to show that 



{n\p^\n) = K^^, {n\p^\n) = K^^^. (G-7) 

Appendix H 

Overlap with dilated functions 



2 



The scalar product of two radial functions Rnj{r) and Rn'.i{r) of a set of orthonormal 
states is simply given by 5„„'. When one of these functions is scaled by a positive factor a, 
the overlap 

poo 

Fn.n'.i{a) =a^^^ Rn,i{x)Rn',i{ax)x^dx (H-1) 

JO 

satisfies the following properties P^S 

limf„,„/,/(a) = 5nn', 

\Fn,n',l{a)\ ^ 1' 

,/(!/«) =Fn',n,l{^)^ 

limF„.„/,/(a) = lim F„_nii{a) = 0. (H-2) 



The first relation stems from the definition (IHTb . the second one from the Schwarz inequal- 
ity, and the others are due to scaling properties. 

Using the dilation properties of the Laguerre polynomials and the various existing recur- 
rence relations,' ^'^^^ it is possible to compute analytically the formula F„ „/ / (a) for hydrogen- 
like systems and harmonic oscillators. 

In the first case (hydrogen-like systems), one obtains 

Fn,l'M = (-1)"+'V^»! {N + l)\n'\ {N' + l)\ [AaNN'Y ^^J^^J'^v" i 
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k=0 



\4aNN' J k\{n-k)\{N-k + l)\{ri-n + k+\)\ 
X (2{N -k){n' -n + k+\) + {n-k){N -k + l)^^ 

+ {n' -n + k){f^ -n + k+\)^^ (H-3) 

withA^ = « + / + l,A^' = n' + / + l, Q{a)=aN-N', andS{a) =aN + N'. 
In the last case (harmonic oscillators), the formula is given by 

F"9M = Vnln'irin + l + 3/2)r{n' + 1 + 3/2) (2a)2«+'+3/2 " 



(1 _|.^2')«+"'+'+3/2 

xVf-n* ^ > CH4) 

;t=o (2a)2^it!(n-it)!(n'-n + yt)!r(n-A: + / + 3/2)' 

Appendix I 

Jacobi coordinates 



Usually the many-body problem in the nonrelativistic framework is treated starting with 
the shell-model or variants. In this method the degrees of freedom are simply the positions 
r,- of the various particles (/ = 1, . . . ,N). But, in this case, the motion of the center of mass 
is not separated correctly and this leads to spurious components which spoil the results. 
On the contrary, in the few-body problem, one introduces new degrees of freedom which 
allow to solve exactly this important drawback. Among the various possible new degrees 
of freedom, the Jacobi coordinates are of common use. First, we choose a reference mass 
m (it can be that of a particle or the total mass of the system for example) and define the 
dimensionless quantities a, = m,/m, ai2...( = ^^k and a = 0L\2...n- Then, the standard 
Jacobi coordinates, x/, can be expressed as 

JCi = ^^^-^-r,+i, / = 1,...,^-1, and x^ = /f = ^^a^r^ (I-l) 

For convenience, the center of mass (for a nonrelativistic system) coordinate R is relegated 

as the last Jacobi coordinate whereas jc,- represents the vector joining the center of mass 
(for a nonrelativistic system) of the first / particles to the particle i+\. Using the matrix 
notation jc, = Ylj^i Uijrj, one is led to the following definition of the ?7 -matrix 

Uij = -^ if j<i 

Uii+i = -l (1-2) 
Uij = if j>i+l. 

It is not difficult to calculate the inverse matrix B = whose matrix elements are given 
by 

Bki = ^'^^ if k<l<N 

«12.../+1 
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Bi+u = - if KN (1-3) 

«12.../+1 

Bki = if k>l + l 

Denoting and JCi the conjugate momenta associated respectively to r, and JC/, it is easy to 
prove that 

N N 

^i=llBjiPj and Pi = Y^UjiTtj. (1-4) 

jTiv = P = /»! + ■ • ■ + PAf is the total momentum of the system. 

Appendix J 

The N-body harmonic oscillator 

The case of a quadratic potential (the A^-body nonrelativistic harmonic oscillator) is 
practically the only one for which an exact solution is reachable in three dimensions, at 
least formally. This section is devoted to this problem, since it is the basic ingredient of the 
A^-body AFM. 

J.l. Nonrelativistic general case 

Let us start with the most general harmonic-oscillator-like Hamiltonian, corresponding 
to N particles of arbitrary masses, a one-body and a two-body quadratic potentials with 
arbitrary spring constants, so that the Hamiltonian looks like 

H^o = 't^ + 'tai{ri-Rf+ £ bij{ri-rjf. (J-l) 

1=1 i ;=1 i<j=l 

The kinetic energy operator, expressed in terms of Jacobi variables, allows the correct 

separation of the center of mass motion and appears decoupled in the various variables (this 
is in fact the justification of the form of Jacobi coordinates) 

N p2 p2 N-l yl ^ 

^ 2m,- 2m, 2m ^' ^ 

where m, = am = mi+m2-\ m^v is the total mass of the system. The kinematical quan- 
tities A, are calculated as 

_ocn^Y\ (J.3) 

For further convenience, it is judicious to switch from the standard Jacobi variables to renor- 
malized conjugate Jacobi coordinates defined as 

y, = |; p,. = A,-;r,-. (J4) 
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Working in the center of mass frame {P = 0) and using these new variables, the kinetic 
energy operator has a very simple form 

2 N-l 

With these new variables, the one-body operator is written 

N N-l 

Vi = £ aiin -Rf= £ Fimyi ■ ym. (J-6) 

where the symmetric definite positive matrix F is defined by 

N 

Flm = h^m 1 ^iBiiBim, (J-7) 
(=1 

the matrix B being given by (II-3I ). In the very same way, the two-body operator is written 

V2= t ^iMi - nf = L '^"nyi -ym (J-s) 

i<j=l l,m=l 

where the symmetric definite positive matrix G is defined by 

N 

Glm=^lKn 1 bij{Bii—Bji){Bim—Bjm). (J-9) 
i<i=l 

Introducing the matrix J = F + G, the total potential V = Vi + V2 is expressed as 

N-l 

v= 1 Jimyi-y„r (J- 10) 

l,m=l 

For arbitrary masses or/and spring constants there is no reason why the matrix / should 
be diagonal. However, this matrix being a symmetric definite positive matrix, it can be 
diagonalized with help of a unitary matrix (in fact an orthogonal one since all the quantities 
are real). Thus 

J = O^DO, with 0^=0. (Ill) 

The elements of the diagonal matrix D are all positive (due to definiteness) and are chosen 
under the form di = m(of /2. 

The last step is an ultimate change of conjugate variables 

N-l N-l 

^i=lLOijyf, Oi=Y^0ijPj. (J- 12) 

7=1 7=1 

Expressed with these new variables, the original Hamiltonian (IJ ll) appears to be the sum of 
N —\ decoupled harmonic oscillators 



N-l 
i=l 



1 2 2 



(J- 13) 
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and consequently the energy of the system is given by 

£ho= £ «/(2?i; + /, + 3/2), (J- 14) 

i=l 

where and are respectively the radial and orbital quantum numbers associated to the 
coordinate z,. The problem is now completely solved. Moreover, this result is valid for any 
excited state. Even if the expression of ft), is in general not analytical, the result (IJ14I) is 



exact and coul d be calculated with a high accuracy. Notice that formula (IJ14I ) extends a 
previous result,'^^''^^ where an equivalent mass formula is obtained in the case a,- = 0. An 
explicit result for the case = 3 is given in Ref . , 



J. 2. Case of identical particles 

It is of interest to rewrite the solution in the case where all the particles are identical, 
implying that they have the same mass m,- = m and the same spring constants a; = a, bij = b. 
Hamiltonian (IJ ll ) indeed reads in this case 

//ho = £|^+a£(r,-/?)' + ^ £ {ri-rjf, (J- 15) 

and it is easy to see that the J matrix is already of diagonal from the very beginning, so that its 
eigenvalues are analytically known (they all read a + Nb). Consequently, the eigenenergies 
of the system are also analytical. Explicitly they are given by 



^ho 



-{a+Nb)Q, (J- 16) 

m 



where Q is the total principal number 



Q=Y,{2ni + h) + l-{N-l). (J-17) 
1=1 ^ 

Wavefunction with good global quantum numbers and various symmetries can be built 
by appropriate linear combinations of states char acterized by the same value of the principal 
quantum number Q^^^ All symmetrized states have then the same energy than the non- 
symmetrized states. For spatial wave functions completely symmetrical, the ground state is 
obtained for the values = Z, = V/ so that the principal quantum number is simply 

esGS = ^(A^-l). (J-18) 

For mixed symmetry or completely antisymmetrical spatial wavefunctions, the situation is 
much more involved. An estimation of the ground-state energy for the completely antisym- 
metrical case can be computed by choosing different values for the quantum numbers and 
piling the states {d identical values of the same nj and Ij per state in order to take care of a 
possible degeneracy due to internal degrees of freedom) up to the Fermi level. By consider- 
ing only particle number insuring a saturated Fermi level (closed shell), one obtains 

QAGS = \{N-l){Bf + 2), (J-19) 



107 



where Bf, the band number of the Fermi level, is the real positive solution of 

N-l = ^{Bf + \){Bf + 2){Bf + 3). (J- 20) 

6 

Asymptotically, we have 

tl\ 1/3^,4/3 



J. 3. Relativistic case 

Let us now treat the relativistic generalization of (IJ15I ). namely 

N I N N 



Hho = T.\/Pi+'^^ + ^ll ('•'• - + ^ L i^i - rjf , (J-22) 

i=l (=1 !<;=l 



by the AFM. Extremization formula (16-631 ) can be written 4X^ - 8X - 37 = with 



{a-vm)Q^ 



(J-23) 



The only positive root of this last equation is given by X = G_(F) (see Sect. lB.2l ). Finally, 
mass formula (16-641) becomes 



^ho(Mo) = l/f^^^(4G_(y) + F) 



Irut 



Notice the simple ultrarelativistic limit 



(J-24) 



lim Mho(Mo) = \ \m[a^m)Q'-\ ^'^ . (J-25) 

The duality relations presented in Sect. I6.5l can be easily tested with relations (IJ-161 ). (I J- 241 ) 
and (IJ-251 ). 
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